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DESCRIPTION AND EVALUATION OF DIGITAL-COMPUTER PROGRAM FOR 
ANALYSIS OF STATIONARY OUTSIDE-COIL LUNDELL ALTERNATORS 

by Gary Bollenbacher 
Lewis Research Center 

SUMMARY 


A digital computer program for analyzing the electromagnetic design of stationary, 
outside coil Lundell alternators is presented. The program, which is written in 
FORTRAN IV programming language, calculates the open-circuit saturation curve, 
field-current requirements at rated voltage for various loads, losses, efficiency, re- 
actances, and weights. The method used to calculate these items is outlined briefly, or 
appropriate references are cited. 

The results of the computer calculations are compared with experimental data. The 
comparison shows that the maximum error for the calculated field excitation require- 
ment is 13 percent. It is shown that this error is due largely to the error in the rotor 
leakage permeance calculation. The comparison further shows that for the efficiency 
calculations the maximum error is 16 percent, this error being caused by the inaccuracy 
of the no-load pole-face loss calculation. Improvements in the method of estimating 
no-load pole-face losses and rotor leakage permeances, although very desirable, are 
beyond the scope of this report. 

The program is useful for parametric studies and alternator design optimization. 
Instructions for using the program and typical program input and output for a 14. 3- 
kilovolt-ampere alternator are given in the appendixes. Also included are an alphabeti- 
cal lists of most FORTRAN symbols and complete program listings with flow charts. 


INTRODUCTION 


The first dynamic energy conversion systems developed at the Lewis Research 
Center for use in space incorporated homopolar inductor alternators. A computer pro- 
gram for analyzing this type of alternator was developed as a result. This program is 



described in reference 1. As more advanced power systems were studied, it became 
apparent that higher rotor speeds and higher gas pressures in the alternator cavity were 
necessary to achieve low weight-to-power ratios. These requirements increased alter- 
nator windage loss to such an extent that it became a critical parameter in system anal- 
ysis. This prompted the selection of the Lundell alternator for advanced, high-speed 
systems. The Lundell alternator offers potentially lower windage loss while still re- 
taining the desirable features of the homopolar inductor alternator: solid-rotor con- 
struction and absence of brushes and of rotating windings. 

With the selection of the Lundell alternator came the need for a Lundell computer 
program, similar to the homopolar-inductor-alternator program, which would facilitate 
parametric system studies and design optimization. The quickest way to develop such 
a program was to modify the existing homopolar program as necessary. Reference 2 
was particularly helpful in this process. 

Most of the necessary modifications stem directly from the differences between the 
two machines. Primarily, these differences alter the magnetic calculations. Other 
modifications include the addition of a more general method for calculating windage loss 
(ref. 3) and a change in the method of accounting for the demagnetizing effect of arma- 
ture currents. While the resultant Lundell program differs substantially from the 
homopolar-inductor-alternator computer program, there still remain many similarities. 
The homopolar and the Lundell programs are each written in FORTRAN IV programming 
language. Most of the FORTRAN symbols used are the same, and the input require- 
ments and output formats are very similar. 

This report describes the program in some detail and evaluates program accuracy 
by comparing the calculated and experimental data for a Lundell alternator used with the 
36 000-rpm 1200-hertz Brayton cycle power conversion system (refs. 4 and 5). Instruc- 
tions for using the program, the program listing and flow charts, and a FORTRAN sym- 
bol list are included in the appendixes. 


COMPUTER PROGRAM DESCRIPTION 
GENERAL DESCRIPTION 

The stationary- coil, outside-coil Lundell alternator computer program is an analysis 
program. This means that the program accepts as input a complete electromagnetic al- 
ternator design; from this, it calculates the open-circuit saturation curve, field-current 
requirements at rated voltage for various loads, losses and efficiency, several react- 
ances, and weights of electromagnetic components. The results of the calculations, 
together with the input, are then printed out to provide a complete, self-explanatory 
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record. An explanation of how to use the program, including a list of input variables 
with definitions, is given in appendix A. A typical program output is shown in appen- 
dix B. 

The program may be used with any computer system which accepts FORTRAN IV. 
For program execution, approximately 13 200 storage locations are needed. At Lewis, 
the program has been used on the IBM 7044-7094 Mod n direct-couple system using a 
FORTRAN IV version 13 compiler. For this system, typical preexecution time is 
1. 25 minutes and typical execution time is 0. 08 minute per alternator design. 

Figure 1 shows a simplified flow chart of the complete Lundell program. As can be 
seen, the program consists of a main program LUNDEL and four subroutines labeled 
SUBLUN, OUTPUT, MAGNET, ANDWINDGE. The subroutines are necessary, in part, 
because one program is too large to compile with the available core storage locations. 
FORTRAN listings and detailed flow charts are given in appendix C. 

Communication between the main program and its subroutines is by means of 
labeled COMMON blocks, except for subroutine WINDGE which uses a call vector. 



Figure 1 . - Simplified flow chart of Lundell alternator computer program. 
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TABLE I. - CROSS-REFERENCE TABLE SHOWING 


USAGE OF COMMON BLOCKS 


COMMON block 

Program or subroutine name j 


LUNDEL 

SUBLUN 

OUTPUT 

MAGNET 

COM1 


X 

X 


COM2 

X 

X 

X 


COM3 

X 



X 

COM4 

X 

X 



COM 5 

X 

X 

X 

i _ x _ . 


Table I lists the COMMON blocks and shows in which surbroutines they occur. 

DESCRIPTION OF ALTERNATOR TO WHICH PROGRAM IS APPLICABLE 

The basic alternator for which the computer program was written is the stationary- 
coil outside-coil Lundell alternator, henceforth for brevity called a Lundell alternator. 
The basic configuration of a Lundell alternator with each major electromagnetic compo- 
nent identified, is illustrated in figure 2. Although the alternator shown has its field 


Yoke (cylindri- 
cal section)--^ 


Yoke (end section)- 





Auxiliary 
air gap 


-Main air-gap 


CD-9847-09 


Figure 2. - Stationary-coil outside-coil Lundell alternator (four pole, two coil version shown). 
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Figure 3. - Exploded view of magnetic rotor parts of four-pole Lundell alternator. 


p Tapered section 



Figure 4. - Alternator rotor. 


068-2241 


winding divided into two identical coils, one coil positioned on each end of the alternator, 
it is also permissible to have the entire field winding located in the center of the alter- 
nator between the stator and the cylindrical section of the yoke. 

The rotor construction for a four-pole Lundell alternator is illustrated in figure 3. 
As shown, the rotor consists of two identical magnetic parts. In the finished rotor the 
two magnetic parts are positioned relative to each other as shown by the broken lines. 
Inserted in the space between them is a nonmagnetic separator. The three pieces are 
welded or brazed together. A finished rotor is shown in figure 4. 

To describe the flux path, a cross-sectional view of a Lundell alternator is shown in 
figure 5. For clarity, the figure shows a two-pole machine. The useful flux (shown by 
solid arrows) starts at a rotor north pole, crosses the main air gap and proceeds radi- 
ally through the stator teeth into the stator back- iron or core. It then goes circumfer- 
entially through the stator core for a distance of one pole pitch, enters the teeth and con- 
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Yoke 





Figure 5. - Cross-sectional view of Lundell alternator showing useful and leakage flux paths. 


tinues radially across the main air gap into a south pole. From there the flux goes axi- 
ally down the rotor and crosses the first auxiliary air gap into the yoke. It leaves the 
other end of the yoke, crosses the second auxiliary air gap, and completes its path back 
to the rotor north pole from which it started. In addition to the useful flux there are 
several leakage fluxes as shown by the broken arrows. 

A number of assumptions, in addition to those implicit in the geometric configura- 
tion, are made regarding the alternator. These assumptions are 

(1) Shaft and poles are made of the same magnetic material. 

(2) Alternator armature winding is three-phase and Y-connected. 

(3) There are no windings or circuits of any kind on the rotor. 

(4) The rotor diameter under the auxiliary air gap is constant. 

(5) The alternator can have only one field winding. The winding may be split into 
two identical coils and positioned as shown in figure 5. Both coils are assumed con- 
nected in series. 

(6) There is no leakage flux across the space between the yoke and the stator 
back iron. 

In contrast to the restrictions imposed on the alternator by the preceding assumptions, 
there are several options that are available to the program user. These options, which 
increase the applicability of the program, are 

(1) Armature conductors and field conductors may be round or rectangular. 

(2) Armature conductors may consist of any number of strands. 

(3) Yoke, rotor, and stator may each be made of a different magnetic material. 
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(4) There are six allowable field coil locations or configurations 

(5) Five different stator slot configurations may be used. 


METHOD OF CALCULATION 

This section of the report will outline in general terms the method of calculation 
used in the computer program. However, due to the length of the program and the large 


TABLE n. - FORTRAN VARIABLES USED IN MAGNETIC CALCULATIONS 



Areas, 

Lengths, 

Flux densities, 

Ampere turns 

Total 

Permeances, 


in. 2 

in. 

q 

kilolines/in. 



flux, b 

lines / amp- turn 




--- 


Sub- 

Nonsub- 

kilolines 





Sub- 

Non sub- 

scripted a 

scripted a 






scripted a 

scripted a 





Main air gap 

GA 

GE 


BG 

FGLL1 

FGL1 

C PPL 

P8 

Auxiliary air gap 

AG2 

G2 


BAG2 

FGLL2 

FGL2 

PPL3 

P9 

Teeth 

ATOOTH 

HS 

BTLL 

BTL 

FTLL 

FTL 

C PPL 

— 

Core 

ACORE 

DCORE 

BCLL 

BCL 

FCLL 

FCL 

C PPL 

— 

Pole 

PA and 

CL 

d BPLL 

BPL 

FPLL 

FPL 

FLUX3 

— 


PAA 








Rotor tapered 

(e) 

DSHMID 

d BROTL 

BROT 

FROTL 

FROT 

PPL2 


section 









Rotor end 

ASHFT 

LG2 

d BSHFTL 

BSHFT 

FSHFTL 

FSHFT 

(e) 

--- 

section 









Pole tip 

PA1 






PPLll 

— 

Pole base 

PA 2 






PPL1 

— 

Yoke end 

(e) 

10. *DYK2 

d BYCLL2 

BYCL2 

FYCLL2 

FYCL2 

f PPL3 or 

' 

section 







PPL4 


Yoke cylindrical 

AYK1 

DYK1 

BYCLL1 

BYCL1 

FYCLL1 

FYCL1 

f PPL3 or 

--- 

section 







PPL4 


Total demagnetizing 






FGML 


— - 

ampere-turns at 









rated load 









Demagnetizing 






FGXL 


— 

ampere turns 









(direct axis) 









Total field 





FFLL 

FFL 



excitation 










a Nonsubscripted variables are used in MAGNET. These variables are then stored in subscripted arrays 
after return to LUNDEL. 
b Total flux = Useful flux + Leakage flux. 

C PPL is the total useful flux per pole. 

d Flux density is not constant. The variable given is the maximum flux density in the element. 
e Area or flux are not constant. See figs. 10 and 11 for assumed variation. 

^Whether PPL3 or PPL4 depends on field coil location. 
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number of equations involved, specific equations will not, except in a few instances, be 
given. These equations can be found in volume 2 of reference 2. In addition, detailed 
information and specific equations may be found in the program listings and flow charts 
in appendix C. To assist in locating specific information in the listing COMMENT cards 
are used freely to identify the major calculations. Of further value are appendix D, 
which is a FORTRAN symbol list, and table II, which gives the most important variables 
used in MAGNET in a convenient format. 

Many variables, both in the report and in the computer program, are expressed in 
the per unit system. Per unit quantities are defined as follows: 

1 per unit voltage = Rated voltage 

1 per unit current = Rated current 

1 per unit volt-amperes = Rated volt-amperes 

1 per unit impedance = (Line-to neutral rated voltage)/(Rated current) 

Some variables are given in percent; these are per unit values multiplied by 100. Equa- 
tions given throughout this report will use the same FORTRAN symbols used in the pro- 
gram and as defined in appendix D. Where a FORTRAN symbol does not exist, an ordi- 
nary algebraic symbol will be used. 


Leakage Permeance Calculations 

The leakage fluxes, shown in figure 5 by the broken arrows, have a pronounced ef- 
fect on the saturation curves and, to a lesser degree, on alternator efficiency and 
weight. Because the leakage paths are through air (or some other gas exhibiting linear 
magnetic properties) they can be characterized by a permeance value. The permeance 
of a leakage path is defined as the ratio of flux through the path to the ampere-turn drop 
across the path. 

Three primary leakage permeances are considered by the program: 

(1) Field coil leakage permeance - This leakage permeance is denoted by P5. Leak- 
age flux through P5 affects only the flux density in the yoke. 

(2) Stator-to-rotor leakage permeance - This permeance, labeled P7, accounts for 
the leakage flux from the stator to the tapered section of the rotor. 

(3) Rotor leakage permeance - This leakage permeance has the greatest effect on 
alternator performance. The flux through this permeance is assumed confined entirely 
to the volume occupied by the nonmagnetic separator in the rotor. The rotor leakage 
permeance is made up of a number of components: PI, P2, P3, P4, and P6 (fig. 6). 

The component PI is the pole tip leakage permeance. The flux through PI is assumed 
to flow through the entire pole from pole tip to pole base while the flux through P6 flows 
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Section A-A 



from rotor end section to rotor end section without ever entering the poles. The com- 
ponents P2, P3, and P4 are distributed over the entire pole. For convenience they are 
combined into one permeance P10 = 2(P2 + P3 + P4). The total rotor leakage permeance 
from one magnetic pole piece to the other is 

^RL = ( px )( pl + P2 + P3 + P4) + P6 
= (PX)^Pl + ~^j + p6 

Because of the many possible geometries, and because all geometries are complex, it is 
difficult to program leakage permeance calculations. For this reason allowance has 
been made for the programmer to calculate leakage permeances as desired and to use 
the values so computed as input to the program. However, as a convenience, approxi- 
mate formulas for the leakage permeances PI, P2, P3, P4, P5, and P6 have been in- 
corporated in the program. The formulas used were taken from reference 6 and are 
applied very loosely to the machine geometry. No attempt was made to calculate P7 
even approximately. Thus, P7 must be read into the program, or it will be assumed to 
be zero. 
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Magnetic Calculations 


The purpose of the magnetic calculations is to compute the flux densities throughout 
the alternator, the ampere-turn drop across the various parts of the magnetic circuit, 
and the total field excitation. In the program, the computations, which are repeated for 
several combinations of load, voltage, and power factor are divided into three cate- 
gories: 

(1) No-load rated-voltage magnetization characteristics. 

(2) Alternator load characteristics. This category comprises magnetic calculations 
at rated voltage for several loads at a fixed power factor. These calculations may then 
be repeated for up to five different power factors. 

(3) No-load saturation data. This category gives up to 10 points on the open-circuit 
saturation curve starting at voltage VMIN. 

The magnetic calculations are carried out, for the most part, in subroutine 
MAGNET which is called by the main program LUNDEL as required (see fig. 1). How- 
ever, prior to calling MAGNET for the first time all areas and lenghts needed by 
MAGNET are calculated in LUNDEL. 

Just prior to each transfer to the subroutine two additional quantities needed by 
MAGNET are calculated in LUNDEL. The two quantities are the actual flux per pole in 
the main air gap (PPL) and the demagnetizing ampere-turns due to armature current 
(FGXL). Both PPL and FGXL are functions of voltage, load, and load power factor. 

At no-load, PPL is proportional to the output voltage, and, since the armature cur- 
rent is zero, the demagnetizing ampere-turns FGXL are zero. Expressed in equation 
form, 


PPL = (EDD)(FQ) 
FGXL = 0 


where 

FQ air-gap flux per pole at no-load rated voltage (calculated in SUBLUN) 

EDD constant of proportionality; here, equal to alternator output voltage normalized 
to rated voltage 

Under load (at rated voltage) the calculation of PPL and FGXL requires reference 
to the alternator phasor diagram shown in figure 7. Construction, of the diagram starts 
with the terminal voltage and line current phasors which, along with the angle between 
them are determined by the particular load for which PPL and FGXL are to be calcu- 
lated. The phasors representing resistance and reactance drops are then added to the 
diagram using values of armature resistance and synchronous reactance values pre- 
viously calculated in SUBLUN. This locates the direct axis. 
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Figure 7. - Phasor diagram. This phasor diagram is in per unit system and is for rated 
voltage and load G at power factor cos (AN). Hence the length of the terminal 
voltage phasor is one, and the length of the armature current phasor is G. Values of 
reactances XL, XB, and XA are in percent. 


The angle between the direct axis and the current phasor is of interest for calculat- 
ing FGXL. This angle, labeled XI and called the internal power factor angle, can readily 
be calculated from the geometry of the phasor diagram. The length OP is also readily 
determined from the phasor diagram. Both PPL and FGXL can now be determined from 
the following equation: 


FGXL = 2 • FGML • G • sin(XI) 
For power factors less than 0. 95, 

PPL = FQ • (length of OP) 


for power factors >0. 95 

PPL = FQ • (length of OP) ' 1. 10 


where 

FQ flux per pole at no-load rated voltage 

FGML demagnitizing ampere-turns at rated load 

G volt-ampere output of alternator normalized to rated volt-ampere 

The PPL and FGXL calculations are summarized in table in. 

With PPL and FGXL known, the calculations are turned over to subroutine MAGNET. 
The approach taken in MAGNET is to represent the actual magnetic circuit of the alter- 
nator by an equivalent electrical circuit having lumped resistive elements. The resistors 
representing air gaps or leakage flux paths are assumed linear. Those which represent 
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TABLE m. - SUMMARY OF CALCULATIONS FOR PPL AND FGXL 



No-load rated- 
voltage magnetization 
characteristics 

Alternator 
load characteristics 
(rated voltage) 

No-load 

saturation 

data 

Flux per pole 
(PPL) 

PPL = FQ 

For power factor <0. 95 
PPL = FQ • (Length 
of OP a ) 

For power factor SO. 95 
PPL = FQ • (Length of 
OP) • 1. 10 

PPL = EDD-FQ 

Demagnetizing 

FGXL = 0 

FGXL = 2 • FGML • G • sin(XI) a 

FGXL = 0 

ampere-turns 




(direct axis) 




(FGXL) 





a OP and XI are obtained from vector diagram shown in fig. 7. 


a part of the flux path through iron must be considered nonlinear. The exact equivalent 
circuit used depends on the value of TYPY, which, in code form, gives the field coil 
location. The permissible field coil locations, the corresponding equivalent circuit, and 
the code value of TYPY are shown in figure 8. The only difference between the six cir- 
cuits is the location of the leakage permeance P5. 

The equivalent circuit of the flux path between the base of one pole and the base of 
the adjacent pole is independent of the field coil location and is given in figure 9. The 
figure shows the outline of two adjacent poles. Each pole is arbitrarily broken into 
10 segments. Flux is assumed to leave each pole segment and to split into two compo- 
nents, the useful flux FLUX1, and the rotor leakage flux FLUX2. The useful flux flows 
across the main air gap, teeth, and core to the opposite segment of the adjacent pole. 

The leakage flux flows in an electrically parallel path through permeance P10. Within 
each pole the flux is assumed to flow axially; once the flux leaves the pole, the flow is 
assumed radial and circumferential, but not axial. These assumptions, when translated 
into the equivalent electrical circuit of figure 9 results in the 10 vertical branches rep- 
resenting the main air gap, teeth, and stator. Similarly, the leakage flux path is as- 
sumed to consist of 10 vertical branches, each having a permeance P10/10. 

The FORTRAN symbols used in figure 9 are as follows: 

FLUXl(I) useful flux leaving 1^ pole segment; from symmetry, FLUX1(I) = FLUXl(ll-I) 

FLUX2(I) rotor leakage flux leaving I th pole segment; from symmetry 

FLUX 2 (I) = FLUX2(11-I) 

XL 

FLUX 3 (I) total flux passing through center of I pole segment 
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Alternator configuration 


Equivalent magnetic circuit 


Yoke (cylindrical section) 


Yoke (end/ 
section)- 


a 


Core 


Teeth 


Field 
cofl— \ 




Pole 


Rotor 


-Rotor centerline - 


Yoke 

(end section)- 
Held coil ^ 


-Main air gap 


I W W, } 

i Yoke (cylindrical section) i 




-Auxiliary / 
air gap-/ 7 


End 


Rotor end 
section 


section 
V 


Tapered! 
, section j 




^-Rotor tapered 
section 

(a) Field coil location code TYPY : 


Rotor 





Main permeance 
Main flux path 
Leakage flux path 



Each leakage permeance is identified 
by its FORTRAN symbol 

Figure 8L - Alternator configurations and equivalent magnetic circuits. 
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(d) Field coil location codeTYPY = 4 


Yoke (cylindri- 
cal section)-* 



AA 




— — Leakage permeance 
— vvv, — Main permeance 

Main flux path 

Leakage flux path 


Each leakage permeance is identified 
by its FORTRAN symbol 


Figures, - Concluded. 





PPL1 



Main air-gap 


Teeth 


Core 


Demagnetizing 

ampere-turns 


Core 


Teeth 


Main air gap 


PPL11 


ATK 10) « FPL - 

Figure 9. - Equivalent circuit of stator, teeth, and poles. 
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ampere-turn drop from one pole segment to opposite pole segment of 
adjacent pole due to flux FLUX1(I) (AT2(I) » AT3(I) after iteration 
has converged) 

ampere-turn drop from pole tip to end of r n pole segment 
rotor leakage permeances 

FPL = AT 1(10) total ampere-turn drop along pole 

FBB ampere-turn drop from base of one pole (point a, figs. 8 and 9) to base 

of adjacent pole (point d’, figs. 8 and 9) 

The equivalent circuit is analyzed by MAGNET to obtain the flux, flux densities, and 
ampere-turn drops for each element. Since the only flux known is the air-gap flux per 
pole, PPL, the logical place to start the analysis is with the equivalent circuit of fig- 
ure 9. The solution of the circuit must satisfy the condition that the magnetic potential 
drop FBB must be the same regardless of the path along which it is computed. For ex- 
ample, paths a-b-b’-d', a-c-c’-d', and a-d-d’ must all give the same result. To meet 

this condition it is necessary to assume a nonuniform air-gap flux distribution, that is, 

10 

FLUX1(I) * FLUXl(J) I * J. However, the total air-gap flux ^ FLUX 1(1) must remain 

1=1 

equal to PPL. 

Determination of the flux distribution requires an iteration process. The steps in 
the iteration are outlined in table IV. Complete details of the iteration are in appendix C 
(program listing and flow chart of subroutine MAGNET). Very briefly the iteration in- 
volves the following steps: 

(1) Initially, assume a uniform air-gap flux distribution. 

(2) Calculate FBB for each of the 10 possible paths. 

(3) Check that 10 values of FBB are within a given tolerance of each other. If 
"yes, ” the iteration is complete; if "no, " proceed to step (4). 

(4) Adjust flux distribution to more nearly equalize all 10 values of FBB, while 
keeping the total air-gap flux equal to PPL. 

(5) Repeat steps (2) to (5) as necessary. 

With the iteration complete, not only is the air-gap flux distribution and ampere- 
turn drop known, but also the flux densities and ampere-turn drops for the teeth, the 
core, and the poles. In addition, the iteration gives the pole-tip leakage flux (PPL11), 
the flux leaving the base of the pole (PPL1), and the ampere-turn drop from pole base 
to pole base (FBB). 

It may happen that the iteration for the air-gap flux distribution fails to converge. 

In that case a uniform flux distribution is assumed and the calculations are continued 
on that basis. A message is then printed on the output record informing the program 


AT 2 (I) 1 
AT 3 (I) J 

AT 1(1) 

PI and! 
P10 J 
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TABLE IV. - DETAILED STEPS IN THE ITERATION FOR AIR-GAP FLUX DISTRIBUTION 


Step 

Explanation of steps in iteration 

Mathematical expression for steps in iteration 

1 

Assume initially that air-gap flux distribution 
is uniform and that flux in each segment 
is 1/10 of total air-gap flux 

FLUXl (1)=FLUX1(2)=. . . FLUX1(10)=PPL/10. 

2 

Calculate ampere-turn drop due to flux 
FLUX 1(1) across main air gap (FGL1), 
teeth (FTL), and core (FCL) for all 10 
paths; the total drop is called AT 3 (I) 

AT3(I)=(FGLl4FTL+FCL+FGXL/2. )*2. 

3 

Calculate pole-tip leakage flux, PPL11 

PPL11=AT3(1)*0. 001*P1 

4 

Calculate rotor leakage flux FLUX2(I) in 
each of 10 paths; the rotor leakage per- 
meance P10 is assumed to be divided 
into 10 equal, parallel paths 

FLUX2(I)=AT3(I)*P10*0. 0001 

5 

With FLUXl, FLUX2, and the pole-tip 
leakage flux known, flux at center 
of each pole segment (FLUX3) is 
calculated 

for 1=1 

FLUX3(I)=PPLll+(FLUXl(I)+FLUX2(I))/2. 
for I > 1 

FLUX3(I)= FLUX 3(1- l)+(FLUXl(l-l)+FLUX2(I-l))/2. 
+ (F LUX 1 (I) + FLUX2 (I)) /2 

6 

Calculate ampere-turn drop from pole- 
tip to end of each pole segment; these 
ampere-turn drops are called AT1(I); 

1 < I< 10; AT 1(10) is also labeled FPL 

See appendix C for details 

7 

Calculate FBB, the ampere-turn drop 
from pole base to pole base, along 
path a~b-b T ~d ? 

FBB= FPL+AT 3(l)-0. 5*(AT I (IK AT 1 (9) -AT 1(10)) 

8 

Calculate what AT3 should be along every 
path (such as a-c-c f -d f ) such that FBB is 
independent of path (i. e. , same value as 
in step 7); call value that AT3 should be 
AT2 

AT2(I)=FBB-2. *FPLh-0. 5*(AT 1(1)4 AT 1(H)) 

IF (I.NE. 10) AT2(I)=AT2(I)4(AT1(II-1))*0. 5 
IF (I. NE. 1) AT2(I)=AT2(I)4(AT1(I-1))*0. 5 
where II = 11 - I 

9 

Compare AT 2 and AT 3 and compute (a) error 
in each path, ERR OR 1(1) (b) cumulative 
error, ERROR 

ERROR1 (I)=AT2 (I) -AT3 (I) 

10 

ERROR = ABS(AT2(I)-AT3(I)) 

1=1 

10 

Recompute FLUX 1(1) such that AT2(I) and 
AT 3(1) will be more nearly equal 

FLUXl(I)=(AT2(I)/AT3(I))**(ERROR3(I)/AA)*FLUXl(I) 
where ERROR3(I) is the value of ERRORl(I) on the 
previous pass and AA=ERROR3(I)- ERROR 1(1) 

11 

Calculate total air-gap flux resulting 
from flux distribution calculated in 
step 10 

10 

FLUX = £ FLUX 1(1) 
1=1 

12 

13 

14 

Recalculate FLUXl (I) such that FLUX=PPL 

If ERROR < 0. 5x10“ 3 , iteration is com- 
plete; proceed to step 15 

If ERROR s 0. 5x10' 3 , return to step 2 
and repeat each step in the iteration pro- 
cedure; but, if, after 25 iterations, 
ERROR is still >0. 5 E-3, go back to 
step 1, execute steps 1 to 7 and then skip 
to step 15 

FLUXl (I)=FLUX1(I)*(PPL/FLUX) 

15 

Calculate total flux through base of pole, 
PPL1 

PPLl=FLUX3(10)4(FLUXl(10) + FLUX2(10))/2. 
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user that a uniform air-gap flux distribution was assumed. 

With the circuit of figure 9 completely analyzed, calculations proceed with the ap- 
propriate circuit of figure 8 to determine fluxes, flux densities, and ampere-turn drops 
for the auxiliary air gap, the yoke, and the end and tapered sections of the rotor. 

In the analysis of the equivalent circuit, a calculation that repeatedly occurs is the 
determination of the ampere-turn drop across an element. This calculation is most 
easily carried out for those elements in which the flux density is constant throughout. 
After calculating the flux density, the ampere-turn drop per inch of element length cor- 
responding to that flux density is found by interpolation between points on the appropriate 
material magnetization curve. Points on the magnetization curve for rotor, stator, and 
yoke materials are input to the program. The ampere-turn drop per inch of element 
length is then multiplied by the length of the element to obtain the total ampere-turn drop 
across that element. 

For the rotor end section, the tapered rotor section, and the yoke end section the 
flux density is not constant throughout. These three elements are each arbitrarily 
broken into 10 segments. The flux density is calculated at the center of each segment 
and assumed constant throughout that segment. The ampere-turn drops for each of the 



Figure 10. - Assumed flux and area variation 
along rotor end section and tapered rotor 
section. 
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10 segments are then calculated and added to obtain the total ampere-turn drop across 
the element. 

To calculate the flux density at the center of each segment requires that the area 
and flux at these points be known. In general, however, area and flux are known only at 
both ends of these elements. Thus, the assumption is made that the area and flux vary 
linearly between the values at both ends. This assumption allows the calculation of flux 
and area, and therefore flux density, at any point along the entire length of the element. 
Figures 10 and 11 illustrate these assumptions of linear variation graphically. At the 
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Figure 11. - Assumed thickness and area variation of yoke end section. 


bottom of figure 10, for example, is the outline of the tapered and end sections of the 
rotor. Directly above the outline drawing are two graphs, one gives the cross-sectional 
area of the rotor and the other gives the flux passing axially through each segment. Note, 
for instance, that for the rotor end section, the flux is assumed to vary linearly from 
zero to the value PPL3 while for the tapered rotor section the flux is constant at PPL2. 
The difference between PPL2 and PPL3 is the leakage flux through permeance P6, all of 
which is assumed to enter the rotor at the interface between the tapered section and the 
end section. Figure 11, which is for the yoke end section, is similar to figure 10. 


Efficiency and Loss Calculations 

Individual losses and efficiency are calculated at rated voltage for several loads of 
increasing magnitude, continuing until the alternator saturates or until calculations have 
been completed for five loads. While the first load at which loss calculations are made 
must always be zero, the program user has the option of specifying any or all of the re- 
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maining four loads. These loads are designated by G within the program (G is the 
volt-ampere alternator output expressed in per unit or percent). 

These calculations are first carried out for rated power factor. They are then re- 
peated up to four times for any other power factor specified in the program input. The 
individual losses that are calculated by the program, along with the method of calculation 
or references, are presented in -the following sections. 

Field conductor losses (PR) and arma ture cond uctor loss es (PS). - These losses 
are given by the expression I^R where I is the direct- current (dc) or root-mean- 
square (rms) current in the winding, as appropriate, and R is the direct current wind- 
ing resistance corrected for the winding temperature. Correcting the winding resistance 
for temperature involves several assumptions: 

(1) The average no-load temperature T NL is known or can be estimated. 

(2) The average rated-load winding temperature T^^ is known or can be esti- 
mated. 

(3) The average winding temperature is a parabolic function of the current in the 
winding. 

With these assumptions, the winding temperature T^ at any load G is 


T G = 


( t rl ' t nl- 


) t g + t nl i rl " t rl i nl 


a RL ~ NL 


current in winding at load G 
current at rated load 

current at no-load, equal to zero for armature winding 

If the programmer, by exercising one of the input options, elects not to have the 
program compute load characteristics for rated load, then the field current at rated load 
I_, T will not be calculated. In that case, when the preceding equation is applied to the 
field temperature calculation, the field current corresponding to the load closest to rated 
load will be used for I RL . 

Armature conductor eddy loss (EX). - References 2 and 7 present a discussion of 
armature conductor eddy losses. In the program these losses are assumed to be zero 
for round conductors. For rectangular conductors they are given by 

EX = (EZ - 1)(PS)(— \ 


where 


*G 


a RL 


a nl 
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where 


PS armature conductor loss 

CL length of armature conductor within stator slot 
HM total armature conductor length 

EZ eddy factor, which is function of conductor geometry, armature winding param- 
eters, frequency, and armature conductor resistivity 

Pole face losses (WN and PP). - The no-load pole-face loss (see refs. 2 and 8) is 
calculated from the equation 


WN = (D1)(D2)(D3)(D4)(D5)(D6)(GA) 


where 



GA 

main air gap area 

. 2 
, in. 

D1 

empirical factor, 

function of 

D2 

empirical factor, 

function of 

D3 

empirical factor, 

function of 

D4 

empirical factor. 

function of 

D5 

empirical factor, 

function of 

D6 

empirical factor, 

function of 


As given in reference 2, the empirical curve which defines the factor D3 is limited to 
a maximum slot frequency of 5800 hertz. In the program, for all slot frequencies 
greater than 5800 hertz, values of D3 are obtained by extrapolation. 

Under load (ref. 9, eq. (22)) the pole face loss is given by 


PP = (k) (WN) k > 1 


where 

k function of alternator current, ratio of slot opening to air gap, number of conduc- 
tors per slot, number of parallel circuits, and air-gap ampere-turns at no- 
load and rated voltage 

Stator core loss (WQL and WQ) and stator tooth loss (ST and WT). - The respective 
equations used to calculate these losses are 
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Stator core loss = 3(Stator core weight) tooth f lux density^ 

\ BK / 


Stator tooth loss = 3 (Stator tooth weight) (WL)('^ :ato r to °^ * lux density^ 

V BK / 

WL core loss at flux density BK and at rated alternator frequency, W/lb 

BK flux density at which WL is measured 

and where weights are given in pounds. 

Miscellaneous loa d loss (WMIS ) . - In an alternator there are generally additional 
electromagnetic losses not accounted for by the losses enumerated so far. These losses, 
labeled miscellaneous load loss, are assumed to be 1 percent of the kilovolt-ampere out- 
put of the alternator at load point G. 

Windage loss (WF|. - If an accurate value for windage loss is known, it may be read 
into the program for use in the efficiency calculation. If the windage loss is not read into 
the program, it will be assumed to be zero. The program user may also elect to have 
the program calculate the windage loss according to reference 3. For this purpose sub- 
routine WINDGE has been included. This subroutine is called twice by subroutine 
OUTPUT: once to calculate windage loss in the main air gap and once to calculate wind- 
age loss in the auxiliary air gap. 

Efficiency ( E). - At each load, efficiency is calculated from 


Efficiency - 


Alter nator power output 

(Alternator power output) + (Losses) 


x 100 


EVALUATION OF C OMPUTER PROGRAM 

Accuracy of the program is evaluated by comparing results of computer calculations 
with experimental data for the following calculations: the rotor leakage permeance, the 
open- circuit saturation curve, the field current at several loads, various losses, and 
efficiencies. The experimental data, taken from reference 5, is for a 1200-hertz 
36 000-rpm alternator. The analytic results for the same alternator come from the 
computer output shown in appendix B. 
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DESCRIPTION OF ALTERNATOR USED FOR PROGRAM EVALUATION 


The alternator used for the evaluation of the computer program is a 1200-hertz 
four-pole Lundell alternator (ref. 5). The alternator was designed for space application 
as part of a Brayton cycle energy- conversion system. The rated alternator output is 
14. 3 kilovolt-amperes at 120/208-volt, 0.75 power factor. 

The alternator rotor (shown in fig. 4) is made of SAE 4340 steel with a stainless 
steel, nonmagnetic separator. No effort was made to reduce pole-face losses by lami- 
nating or grooving the pole faces. Figure 4 also shows that the pole width at the tip is 
less than at the base. The effect of this is that the pole embrace is not constant; for 
computational purposes the average pole embrace is used. The stator is made of 0. 004- 
inch AL 4750 laminations. SAE 1010 steel is used for the yoke. 

For voltage regulation, the alternator has two separate field windings. Each wind- 
ing is split into two coils and one coil of each winding is located on each side of the 
stator. For the test results reported in reference 5 both field windings were connected 
in series. 

When used with the Brayton cycle system the alternator cavity will be filled with 
helium-xenon mixture having a molecular weight of 83. 8. The pressure in the cavity 
will be 44. 4 psia at a temperature of 177° C. These conditions, coupled with the 
36 000 rpm shaft speed, illustrate the importance of having a smooth rotor to minimize 
windage losses. 

Other alternator parameters and dimensions may be obtained from appendix B. 


COMPARISON OF EXPERIMENTAL AND CALCULATED RESULTS 
Rotor Leakage Permeance 

The values for the rotor leakage permeances PI, P2, P3, P4, P6, and P10 as cal- 
culated by the program are given in the first column of table V. When these values are 
used to compute the total rotor leakage permeance the result is 

^RL = ( px H pl + P3 + P2 + P4) + P6 

= (4) (1 . 04 + 2. 43 + 6. 22 + 3. 22) + 0. 41 

= 52. 95 lines/ampere-turn 
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While it is not possible to measure the individual components of it is possible to 

measure approximately. This was done as described in appendix E. The experi- 
mental value obtained was = 108 lines per ampere-turn. This shows that the true 

value of is considerably greater than the calculated value of 52 lines per ampere- 

turn, although, for reasons discussed in appendix E, it is also slightly less than the ex- 
perimental value of 108 lines per ampere-thrn. 

For the purpose of evaluating the computer program it was decided to obtain two 
sets of computer results: one set (shown in appendix B) for which the rotor leakage 
permeances are as calculated by the program and a second set for which the values of 
rotor leakage permeance shown in the second column of table V were used as input to the 


TABLE V. - TWO SETS OF ROTOR LEAKAGE PER- 
MEANCE VALUES FOR WHICH COMPUTER 
OUTPUT WAS OBTAINED 


Permeance, 
lines/am p- turn 

Value calculated 

by computer 
a 

program 

Values* 3 based on 
experimental data 

PI 

1. 04 

2. 16 

P2 

6. 22 

12. 94 

P3 

2. 43 

5. 05 

P4 

3. 22 

6. 70 

P6 

. 412 

. 86 

P10 

[ = 2- (P2+P3+P4) | 

23. 74 

49. 38 

9 

RL 

[=PX- (P1+P10/2. )+P6] 

52. 05 

C 108 


a From appendix B. 

^Except for > ^ ese values are obtained by multiplying 
the previous column by 108/52 = 2. 08. 
c Experimentally determined value (see appendix E). 


program. The second column of table V was obtained by multiplying each permeance 
value in the first column by the ratio 108/52 = 2. 08. 


Open-Circuit Saturation Curve 

The experimental and two calculated open- circuit saturation curves are shown in 
figure 12. The two sets of calculated points correspond to the two sets of rotor leakage 
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Experimental (ref. 5) 
Calculated using ex- 
perimental values 
for rotor leakage 
permeances 
Calculated using 
values for rotor 
leakage permeance 
calculated by pro- 
gram 


500 1000 1500 2000 

Field excitation, ampere-turns 


2500 


Figure 12. - Comparison of experimental and 
calculated open-circuit saturation curve. 


permeances given in table V. In general, the calculated values of field excitation are 
too high at low terminal voltages. The discrepancies, 115 and 150 ampere-turns, re- 
spectively, for the two sets of calculated values, are independent of terminal voltage 
from 20 to 90 percent of the rated voltage. This has the effect of making the calculated 
air-gap lines parallel to, but shifted to the right of, the experimental air-gap line. At 
high terminal voltages, as the alternator becomes saturated, the calculated values of 
field excitation are lower than experimentally observed. 

The effect of changing the rotor leakage permeance by a factor of 2. 08 is very 
minor at voltages less than rated voltage. At higher voltages the effect becomes more 
pronounced with the experimental value of rotor leakage permeance definitely giving 
better agreement with the experimental data. For example, at 135 percent of rated 
voltage, the calculated field excitation, using the values of rotor leakage permeance as 
calculated by the program, is low by 18 percent. This error is reduced to less than 
8 percent by using the experimental value of rotor leakage permeance. 

At rated voltage, the respective errors for field excitation, corresponding to the 
calculated and experimental values of rotor leakage permeance, are 2. 3 and 7 percent. 


Field Excitation Under Load 

Experimentally determined field excitation requirements as a function of load for 
two power factors are plotted in figure 13. Superimposed on each curve are two sets of 
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(a) Power factor, 1.0. (b) Power factor, 0.75 lagging. 

Figure 13. - Comparison of experimental and calculated field excitation requirements. 


calculated points which correspond to the two sets of rotor leakage permeances in 
table V. As can be seen from the figure, the higher value of rotor leakage permeance 
gives calculated values of field excitation higher than the experimental data while the re- 
verse is true for the lower value of rotor leakage permeance. As discussed earlier, the 
true value of rotor leakage permeance lies between those values given in the two columns 
of table V; thus, the corresponding calculation for field excitation should be intermediate 
in value between those plotted in figure 13. In other words, using a more accurate value 
of rotor leakage permeance should improve agreement between calculated and measured 
field excitation requirement under load. 

While the inaccuracy of the rotor leakage permeance calculation does introduce an 
error into the field excitation calculation, the result nevertheless shows fair agreement 
with experimental data. The maximum error, using the values of rotor leakage per- 
meance calculated by the program, is less than 13 percent. 


Losses and Efficiency 

The losses that will be compared for the purpose of evaluating program accuracy are 
the armature conductor loss, the field conductor loss, the open-circuit core loss, and 
the stray-load loss. The open-circuit core loss and the stray-load loss, though com- 
monly measured, are not specifically calculated by the computer program. To permit 
comparison of the open-circuit core loss and the stray-load loss reported in reference 5 
with those losses calculated by the program, the following equivalences are assumed. 
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Open-circuit core loss = (No-load pole-face loss) + (No-load stator tooth loss) 


+ (No-load stator core loss) 


= PP(1) + ST(1) + WQL(l) 

Stray-load loss = (Armature conductor eddy loss) + (Miscellaneous load loss) 

+ (Pole face, stator tooth, and stator core losses due to load) 

= [EX(J)j + [WMIS(J)j + |[PP(J) - PP(1)] + [ST(J) - ST (1)] + [WQL(J) - WQL(1)]J 

= [EX(J) + WMIS(J) + PP(J) + ST(J) + WQL(J)] - [Open-circuit core loss] J > 1 

While the stray-load loss equivalence is not completely valid, it is nevertheless useful 
for the purpose of evaluating program accuracy. 

Field conductor and armature c onductor losses . - The accuracy of the ohmic loss 
calculations for the armature and field windings is related directly to how precisely the 
winding resistance, the winding temperature, and the winding current are known. The 
accuracy of field-current calculations was discussed in the previous section. A discus- 
sion of armature current accuracy has no relevance since (in the program) armature 
current is essentially an independent variable, which is chosen rather than calculated. 
Thus, the comparison of calculated and measured winding losses will be preceded by an 
examination of winding resistances, and of winding temperature variation with current. 

Table VI gives both experimental and calculated values for the field and armature 
winding resistances at 25° C. The calculated results agree quite well with the measured 


TABLE VI. - COMPARISON OF CALCULATED 


AND MEASURED WINDING RESISTANCES 


Winding 

Measured 21 resistance 

Calculated resistance 


at 25° C, 

at 25° C, 


ohm 

ohm 

Field 

4. 882 

5. 118 

Armature 

.0307 

.0309 


a Obtained with a Kelvin bridge. 
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values; the error of field resistance calculation is less than 5 percent; the armature re- 
sistance calculation error is less than 1 percent. 

Winding temperature is estimated in the program, as described in the section 
METHOD OF CALCULATION, by assuming that the winding temperature is a parabolic 
function of winding current. The validity of this assumption is illustrated in figure 14. 

Of course, the accuracy of estimating temperatures in this manner depends largely on 
the accuracy of the no-load and rated- load temperatures read into the program. 


Power factor 



(a) Field winding. (b) Armature winding. 

Figure 14. - Temperature variation with current for field and armature windings. 
Program assumes temperature variation obtained by parabolic interpolation 
between points marked X. Temperature at points marked X are input to 
program. 



Figure 15. - Comparison of experimental and calculated field and armature losses. 
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A comparison of field and armature losses is made in figure 15. The maximum 
error in the field loss calculation is less than 13 percent. Some of the discrepancy be- 
tween calculated and measured field losses is due to the error in the field-current cal- 
culation. If this error were eliminated, the field- loss calculation would be more accu- 
rate at high loads and less accurate at low loads. Armature loss calculation, on the 
other hand, are within 2^ percent of experimental data. 

Open-circu it core loss. - The calculated open-circuit core loss, as defined earlier, 
is given by 


Open-circuit core loss = PP(1) + ST(1) + WQL(l) 

= 959 + 64 + 140 
= 1163 W 

The measured value given in reference 5 is 340 watts. Since the calculated value ex- 
ceeds the measured value by 823 watts, it is clear, from examining the magnitudes of 
the three losses, that the only component which can be in error by that amount is PP(1), 
the no-load pole-face loss. More specifically, the calculated value of no-load pole-face 
loss is too high by at least 619 watts. This error predominates any other possible error 
in the open- circuit core loss calculation to such an extent that future improvement in the 
open-circuit core loss calculation must begin with the no-load pole-face loss. 

One contributing factor to this large error may be the high slot frequency (21 600 Hz) 
of the alternator used for the program evaluation. Extrapolation to that frequency of the 
empirical equations used to calculate the no-load pole-face loss may not be valid. 

The large error in the no-load pole-face loss calculation will also affect the accu- 
racy of the stray-load loss, the total loss, and the efficiency calculations discussed in 
the next sections. 

St ray-load loss. - With the method of calculation used by the computer program the 
pole-face loss at any given load is directly proportional to the no-load pole-face loss. 
Thus, the error in the no-load pole-face loss calculation discussed in the previous sec- 
tion causes a large error in all pole-face loss calculations. This error, in turn, affects 
the stray-load loss calculations. To make a meaningful comparison of measured and 
computed stray- load loss, the calculated values of pole-face loss shown in appendix B 
were first adjusted downward to the value they would have been, had the no-load face loss 
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60 


Figure 16. - Comparison of experimental and calculated stray-load 
loss. 


been 136 watts. a These adjusted values were then used to compute the stray-load loss. 
The result is shown in figure 16. It is evident from the figure that the calculated values 
of stray-load loss are a function of load power factor, while the experimental values, 
which are measured with the alternator short-circuited, are not. This is not unexpected 
since the assumed equivalence of stray-load loss and the sum of several calcluated losses 
is not completely valid. Other sources of error in the stray-load loss calculation are 
errors in the individual losses that add to give the calculated stray-load loss. At rated 
current the error ranges from 37 percent at 0. 75 power factor to 54 percent at unity 
power factor. 

Electromagnetic efficiency . - Experimental and calculated efficiencies are com- 
pared in figure 17. Two sets of calculated points are shown. One set shows the effi- 
ciency points exactly as calculated by the program, that is, using the obviously incorrect 
value of no-load pole-face loss and the correspondingly high values of stray-load loss. 

The other set of points was computed using the adjusted no-load pole-face loss and the 
corresponding values of stray-load loss. The computed efficiencies, using the adjusted 
pole-face loss, differ from the measured efficiencies by no more than 2. 4 percent out of 
92 percent for an error of only 2. 7 percent for the range of 50 to 100 percent of rated 
kilovolt-ampere. 

a Assuming that the no-load core and tooth loss calculations are correct, 136 watts 
must be the value of the no-load pole face loss to give agreement between the experi- 
mental and calculated values of open-circuit core loss. 
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Figure 17. - Comparison of experimental and calculated efficiencies. 


CONCLUDING REMARKS 


This report describes a digital computer program for the analysis of Lundell alter- 
nators. The accuracy of the calculations is evaluated by comparing the computer results 
with experimental data for one alternator. The evaluation begins by comparing the cal- 
culated and measured values of rotor leakage permeance. The comparison shows, not 
unexpectedly, that the calculated value is very approximate, being low by a factor of two 
for the particular alternator used in the evaluation. This error affects the accuracy of 
the field- excitation calculations. It is shown that elimination of the error would improve 
the agreement between experimental and calculated open circuit saturation curves and 
between experimental and calculated field excitation requirements under load. However, 
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even without this correction, the maximum error in field excitation at rated voltage is 
13 percent. 

Field and armature loss calculations are accurate to within 13 and 3 percent, re- 
spectively. Other loss calculations could not be evaluated precisely because of a lack of 
detailed experimental data. However, it is clear that the no-load pole-face loss calcu- 
lation is greatly in error for the alternator used in this evaluation. One possible expla- 
nation for this large error is that the empirical method of calculation used in the pro- 
gram is not valid for the high slot frequency encountered in this alternator. 

The error in the no-load pole-face loss calculation in turn affects the accuracy of 
the stray load loss and efficiency calculations. Thus, caution must be used in inter- 
preting the efficiency results, especially if the calculated no-load pole-face loss is high. 
When allowance is made for the inaccuracy in the no-load pole-face loss and its effect on 
stray- load loss, efficiency calculations are close to the experimentally observed value. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 14, 1969, 

120-27. 
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APPENDIX A 

HOWTO USE COMPUTER PROGRAM 
Input Data Requirements 

To use this computer program for the analysis of Lundell alternators, the complete 
electromagnetic design of the alternator must be known. This includes physical dimen- 
sions, armature and field winding parameters, and the magnetic characteristics of the 
materials to be used in the stator, rotor, and yoke. The design information must then 
be transferred onto data cards for use with the program. A typical set of data cards is 
shown in figure 18. It consists of three material decks followed by any number of alter- 
nator design decks. The material decks must be in the order shown in the figure, that 
is, stator material, rotor material, and yoke material. There must be exactly three 
material decks in each data deck even if two or all three materials are identical. The 
last card of the alternator design deck ($WIND) is needed only if windage loss is to be 
calculated by the program. 

If more than one alternator design deck is included in the data deck, the program 
will treat each design deck independently. Each will result in a separate alternator anal- 
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ysis complete with an individual output record. However, the same three material decks 
will be assumed to apply to each alternator design deck. 


Preparation of Material Decks 

A material deck consists of five cards. The first card contains the material name. 
This serves two functions: it identifies the material deck, and it is read by the computer 
and stored for later printout on the output record. The remaining four cards contain 
the coordinate values of 14 arbitrary data points located on the magnetization curve of 
the material specified on the first card. 

During program execution, the original magnetization curve is approximately re- 

TABLE VII. - DEFINITION OF INPUT VARIABLES 


USED WITH MATERIAL DECKS 


Card 

FORTRAN symbol 

Format 

Description of data 

1 

SMAT 

6A6 

Stator material name 

2 - 5 

AI{J), lfijs29 

8F10. 1 

AI(1) is the saturation flux density of the stator 
material; AI(2) to AI(29) are, alternately, 
values of flux density (kilolines/in. ) and 
magnetizing force (A-turn/in. ) for 14 points 
on the stator material magnetization curve; 
the 14 points are arranged in order of in- 
creasing flux density; note that AI(1) = AI(28) 

6 

RMAT 

6A6 

Rotor material name 

7-10 

AI(J), 31sJ<59 

8F10. 1 

AI(31) is the saturation flux density of the rotor 

material; AI(32) to AI(59) are, alternately, 

2 

values of flux density (kilolines/in. ) and 
magnetizing force (A-turn/in.) for 14 points 
on the material magnetization curve; the 
14 points are arranged in order of increasing 
flux density; note that AI(31) = AI(58) 

11 

YMAT 

6A6 

Yoke material name 

12 - 15 

Al(J),61sJfi89 

8F10. 1 

AI(61) is the saturation flux density of the yoke 

material; AI(62) to Al(89) are, alternately, 

2 

values of flux density (kilolines/in. ) and 
magnetizing force (A-turn/in. ) for 14 points 
on the yoke material magnetization curve; 
the 14 points are arranged in order of in- 
creasing flux density; note that 
AI(6 1) = AI(88) 
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constructed by interpolation between points. The interpolation assumes a straight line 
on semilog paper between data points. 

Table VII lists those variables associated with the three material decks. It also 
states the FORMAT used on each card and gives a brief description of the data appearing 
on each data card. 

To illustrate preparation of a material deck, SAE 4340 steel will be used as an ex- 
ample. The first card of this material deck will appear as shown in figure 19. The 
material name should start in column 1 and may extend up to column 36. 



To prepare the remaining four cards of the material deck, the magnetization curve 
of the material is needed. The magnetization curve for SAE 4340 steel is shown in fig- 
ure 20. The units must be kilolines per square inch for the magnetic flux density and 
ampere-turns per inch for the magnetizing force. Fourteen points on the curve must 
then be chosen. In the figure, 13 points are indicated by data symbols; the 14th point is 
off the graph. These points are listed in the table insert. Careful attention must be 
paid to the sequence in which the numbers are punched onto data cards. The first 
number must be the maximum flux density of the points chosen. In the example, this 
value is 122. 5 kilolines per square inch. This is followed in ascending order, by alter- 
nate values of magnetic flux density and magnetizing force. Again, in the example, with 
reference to the table insert, the values appear in the following sequence on the data 
cards: 122.5, 0, 0.01 6.45, 12, 12.9, 24, . . . , 116, 280, 122.5, 454. The complete 
material deck for SAE 4340 steel is shown in figure 19. 
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Figure 20. - Average magnetization curve for SAE 4340 steel. 


Preparation of Alternator Design Deck 


The alternator design deck contains all the dimensions, the geometric configuration 
(in numerical code), and the winding parameters needed for an electromagnetic analysis 
of the alternator design. Unlike the material decks, which are read according to a 
FORMAT statement, the alternator design decks are read with a READ statement ref- 
erencing a NAMELIST name. For each NAMELIST name one or more data cards are 
required to numerically define the variables included in that NAMELIST name. In all 
there are 10 NAMELIST names. Each name is suggestive of the type of variables in- 
cluded in its list. Detailed information about each NAMELIST name is provided in 
table vm which lists all variables used with the alternator design deck. All variables 
belonging to the same NAMELIST name are grouped together. The NAMELIST names 
are arranged in the order in which the data cards must appear in the data deck. Units 
are given, where applicable, and each variable is classified as mandatory (M), condi- 
tional (C), or optional (O). A mandatory classification indicates that the variable must 
be read in. The conditional classification indicates that, for some alternator designs, 
the variable is required and that, for others, it may be omitted from the program input. 
Variables identified as optional are read in at the discretion of the program user. In 
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TABLE VIE. - DEFINITION OF INPUT VARIABLES USED WITH ALTERNATOR DESIGN DECK 


NAMELIST 

FORTRAN 

Definition 

Classification 


Remarks 

name 

symbol 


(a) 



RATING 

VA 

Kilovolt-ampere rating of alternator, 

M 





kVA 






EE 

Line-to-line design voltage, rms V 



Either one must be read in, or both may 


EP 

Line-to- neutral design, voltage, 


zi 

read in 



rms V 






F 

Frequency, Hz 






RPM 

Shaft rotational speed, rpm 


4 

Any two must be read in, or all 


IPX 

Number of poles 


4 

three may be read in 


G 

Volt-ampere load (normalized to 


D 

G is a subscripted variable (array size 



rated VA) at which load char- 



is 5); if not read in, program assumes 



acteristics are calculated (see 



values, 0, 0.75, 1.0, 1.25, and 1. 50; 



appendix B, pp. 51 to 53), 



any one or all (except 0) may be 



percent or per unit 



changed by reading in different values; 






program automatically arranges 






values in increasing order; any num- 






ber >9. 0 is assumed to be in percent, 






<9 

0 in per unit. If, at some value 






of 

G, the alternator saturates, then 






G 

will be decreased in steps of 0. 1 






until saturation no longer occurs 


PFC 

Power factors at which load char- 

M 

PFC 

is a subscripted variable (array 



acteristics are calculated (see 



size is 5); one or more values 



appendix B, pp. 51 to 53) 



must be read in; the first value will 






be assumed to be the design power 






factor 


VMIN 

Lowest voltage for which a point 


3 

If not read in, VMIN = 0.70 will be 



on the open-circuit saturation 



assumed 



curve will be calculated (see 







appendix B, p. 54), per unit 





STATOR 

DI 

Stator inside diameter, in. 

M 




DU 

Stator lamination outside 

M 





diameter, in. 






| CL 

Length of stator stack, in. 

M 




; HV 

Number of cooling ducts 


If there are no cooling ducts, these 


BV 

Width of cooling duct, in. 


'i 

need not be read in 


SF 

; Stacking factor (stator) 


Either one or both may be read in; if 


LTS 

Stator lamination thickness, in. 

CJ 

neither is read in, program assumes 






that stator is not laminated 






(SF 

= 1.0) 


WL 

Core loss at flux density BK 

M 





W/lb 






BK 

Flux density at which core loss 

M 





WL is given, kilolines/in. ^ 





SLOTS 

zz 

Slot type 

M 

See fig. 21 1 


BO 

Slot dimension (use with slot 

C 





type 2, 3, 4), in. 






B3 

Slot dimension (use with slot 







type 3), in. 






BS 

Slot dimension (use with slot 







type 1, 2, 4, 5), in. 






HO 

Slot dimension (use with slot 







type 2, 3, 4), in. 






HX 

Slot dimension (use with slot 







type 1, 2, 3, 5), in. 






HY 

Slot dimension (use with slot 







type 1, 2, 3, 5), in. 






HS 

Slot dimension (use with slot 

M 





type 1, 2, 3, 4, 5), in. 






HT 

Slot dimension (use with slot 

c 





type 2, 3), in. 





| IQQ 

Number of slots 

M 




a M, mandatory; C, conditional; O, optional. 



TABLE VIE. - Continued. DEFINITION OF INPUT VARIABLES USED WITH ALTERNATOR DESIGN DECK 


NAMELIST 

FORTRAN 

Definition 

Classification 

Remarks 

name 

symbol 


(a) 


WINDNG 

RF 

Type of coil 

M 

RF =1.0 for form wound coil; RF = 






for random wound coil 


SC 

YY 

Number of conductors per slot 
Slots spanned per coil (number of slots 



In fig. 24, SC = 4 






between coil sides plus one) 

Number of parallel circuits per phase 







DW 

Strand diameter or width, in. 



See fig. 24 


SN 

Strand per conductor in depth (radial 

C 


Read for rectangular wire only (in 



direction) 



fig. 24, SN = 2) 


SN1 

Total strands per conductor 

M 

In fig. 24, SN = 4 


DW1 

Uninsulated stator strand thickness 

C 


Read for rectangular wire only; see 



(radial direction), in. 



fig. 24 


CE 

Straight portion of coil extension, in 

M 

See fig. 24 


SD 

Distance between centerline of strands 

M 

See fig. 24 



in depth, in. 





PBA 

Phase belt angle, deg 

C 


If not read in, program assumes 






PBA = 60° 


SK 

Stator slot skew at stator inside 

c 


If not read in, program assumes 



diameter, in. 



SK = 0 


T1 

Rated-load armature winding tern- 

M 

Used for loss and efficiency 



perature, °C 



calculations 


RS 

Armature conductor resistivity 

C 


If not read in, program assumes 



at 20° C, (^ohm)(in, ) 



copper resistivity (0. 694) 


ALPHAS 

Armature conductor temperature 

c 


If not read in, program assumes 



coefficient of resistivity at 20° C, 



copper temperature coefficient 



(“c)' 1 



(0. 00393) 


Til 

No-load armature winding tempera- 

M 

Used for loss and efficiency cal- 



ture, °C 



culations 


TST 

Armature winding temperature, °C 

C 


; Program calculates and prints out 






armature resistance at this tem- 






perature; if not read in, program 
assumes TST = 25° C 


EL 

End turn length, in. 

c 


Read in if exact value is known; if 






not, program will calculate ap- 
proximate value 

AIRGAP 

GC 

Main air gap (radial dimension), in. 

M 

See fig. 23 


G2 

Auxiliary air gap (radial dimension), 
in. 

Auxiliary air-gap length (axial 

M 

See fig. 23 


LG 2 

M 

See fig. 23 



dimension), in. 





WF 

Windage losses, W 

O 


Read in actual value; if not read in, 

! 





program assumes windage to be 
zero; to have program calculate 
approximate windage loss, set 
WF = 1. 0 and add data card for 
NAMELIST name WIND to end of 
alternator design deck 


a M, mandatory; C, conditional; O, optional. 



TABLE vm. - Continued. 


DEFINITION OF INPUT VARIABLES USED WITH ALTERNATOR DESIGN DECK 


NAMELIST 

name 

ROTOR 


FORTRAN 

symbol 

PL 

PE 

WROTOR 

D1 


Cl 


CM 

CQ 


PW1 

PW2 

PTl 

PT2 

PA 


PA 1 
PA 2 

GAMMAR 

DISHl 


Definition 


Pole body length (axial direction), in. 
Pole embrace 
Rotor weight, lb 
Pole face loss factor 


Ratio of fundamental maximum to ac 
tual maximum value of field form 
(fielu form is air-gap flux density 
distribution due to field only) 

Ratio of average to maximum value 
of field form 

Demagnetizing factor (direct axis) 

Cross magnetizing factor (quadra- 
ture axis) 

Inside shaft diameter (for hollow 
shaft), in. 

Pole width at narrow end, in. 

Pole width at wide end, in. 

Pole thickness at narrow end, in. 

Pole thickness at wide end, in. 

o 

Pole cross-sectional areas, in. 


Pole cross-sectional area at 
pole-tip, in. ^ 

Pole cross-sectional area at 
2 

base of pole, in. 

Density of nonmagnetic rotor 
material, lb/in. ® 

Shaft diameter under auxiliary 
air gap, in. 

Rotor length at diameter DR, in. 


Rotor length between rotor end 
sections, in. 


Classification! 

(a) 


C 

O 

o 

o 


c 

M 


c 

c 

o 

M 

o 


See fig. 23; omit if rotor leakage per- 
meances rre input to program 
If not read in, program calculates PE for 
pole widt s (FW1+PW2)/? 

If not read in, program will calculate an 
approximate rotor weight 
If not read in, D1 = 7. 0; factor is used for 
no-load pole-face loss calculations; for 
a solid rotor factor is normally equal to 7 
(see ref. 2). Setting D1 = 0 reduces 
pole-face losses to zero 


Identical to those defined for conventional 
salient pole alternator (ref. 10); if not 
an input, values are calculated from 
formulas given in refs. 2 and 10 


Read in only for hollow shaft (see fig. 23) 


PA is a subscripted variable (array size 
is 11); PA(1) is the pole lip area: 

PA (11) is the area at the base of the pole; 
PA(2) through PA (10) are defined at equally 
spaced intervals between PA(1) and PA(ll); 
PA(1) must be read in unless PA1 is read; 
PA (11) must be read in unless PA2 is read 
in; any or all of the remaining PA's may 
be omitted. In that case the computer will 
fill in the missing values by linear inter- 
polation 

Need not be read in if PA(1) is read in 

Need not be read in if PA(ll) is read in 

If not read in, GAMMAR = 0. 283 will be as- 
sumed; if WROTOR is read in, then 
GAMMAR is not needed 
See fig. 23 

If not read in SHL1 = CL will be assumed 
(used for weight calculation only); see 
fig. 23 
See fig. 23 


Mj 


mandatory; C, conditional; O, optional. 



TABLE vm. - Concluded. DEFINITION OF INPUT VARIABLES USED WITH ALTERNATOR DESIGN DECK 


NAMELIST 

name 


FORTRAN 

symbol 


TY 

Yoke dimension, in. 

M 

TYL 

Yoke dimension, in. 


DYC 

Yoke outside diameter, in. 


TYE 

Yoke dimension, in. 


TYR 




TY2 



C 

TY3 



M 

WYOKE 

Yoke weight, lb 

O 


PCOIL 

DCOIL 

PT 

RD 

RT 

BCOIL 

T2 

T22 

RR 

ALPHAR 

TF 


P2 

P3 

P4 

P5 

P6 

P7 

P10 


RO 


VIS 


(Classification 
(a) 


Field coil inside diameter, 
in. 

Field coil outside diam- 
eter, in. 

Number of field turns per 
coil 

Field conductor diameter 
or width, in. 

Field conductor thickness, 
in. 

Field coil width, in. 

Rated -load field tempera- 
ture, °C 

No-load field temperature, 

°C 

Field-coil resistivity at 
20° C, (pohm)(in. ) 

Temperature coefficient of re- 
sistivity at 20° C, °C * 

Field-coil temperature, 


Field coil location code 


Leakage permeance, 
lines/am-turn 


Density of gas in alternator 
cavity at pressure PRESS 
and temperature TEMP, 
lb/ft 3 

Viscosity of gas in alternator 
cavity at temperature TEMP 
and pressure PRESS, 
Ib/sec-ft 

Pressure of gas in alternator 
cavity, lb/in. ^ 

Molecular weight of gas in al- 
ternator cavity 

Temperature of gas in alter- 
nator cavity, °C 


O 

C 

O 


Remarks 


See fig. 22 


If not read in, TY2 = TY3 will be assumed 
(see fig. 22) 

See fig. 22 

If not read in, program will calculate an 
approximate value 


Do not read in for round conductors 


Used in loss and efficiency calculations 


If not read in, 0. 694 is assumed 

If not read in, 0. 00393 is assumed 

Program calculates and prints out field - 
coil resistance and open-circuit time 
constant at this temperature; if not 
read in, program assumes TF - 25° C 
See fig. 8 

See figs. 6 and 8; if these permeances 
are omitted from program input, the 
program will calculate approximate 
values except for P7; if P7 is omitted 
it will be assumed zero; if P10 is read 
in, P2, P3, and P4 may be omitted 


If RO is not read in, then M must be 
read in; program then calculates RO 


If VIS is not read in, program assumes 
viscosity of air at 25° C (= 1. 24x1 0~ 3 ) 


If not read in, program assumes 
PRESS = 14.7 

This must be read in if RO is not read 
in; otherwise it may be omitted 
If not read in, program assumed 
TEMP = 25 


a M, mandatory; C, conditional; O, optional. 

^This card is needed only if WF of NAMELIST name AIRGAP is unity. 



i 


r 

HX 


□ 

.□ 

- BS 
■ BO* 


1 

HS 


{BS = BO) 


. r ? j g 


HXT 


o 

□ 

BS 


7 

HS 


Types 1 and 5: Open 
slot, constant slot 
width. Type 5 slot is 
same as type 1, but 
it contains only one 
coil side. 


Type 2: Partly closed slot, 
constant slot width. 



Type 3: Partly closed slot, Type 4: Round slot, 

constant tooth width. 


Figure 21. - Stator slot dimensions. (Starred variables are not 
input, they are shown for reference only. ) 



Figure 22. - Yoke dimensions. 



PT1, PW1, PA1 


PT2, PW2, 
PA2 and 




Figure 24. - Definition of variables used with NAMELIST name WINDNG. 


each case where an optional variable is omitted from the program input, an assumption 
regarding that variable is made internal to the program. The remarks column of the 
table supplies specific explanations in regard to all conditional or optional variables. 

To further clarify the definition of some variables, figures 21 to 24 are given. 
These figures are referenced in the table where applicable. 

Preparation of an alternator design deck will be illustrated with the construction of 
a typical data card for the NAMELIST name FIELD. The data that will be used is for 
the 1200-hertz 14. 3-kilovolt-ampere 120/208-volt Brayton cycle alternator (refs. 4 
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and 5). Figure 25(a) gives all pertinent design data for the Brayton cycle alternator field 
circuit. Figure 25(b) shows how the design data are related to the variables of 
NAMELIST name FIELD and how these data are transferred to the data card $ FIELD. 

Data cards for the remaining NAMELIST names are prepared in a similar manner. 
To illustrate the results, a complete data deck listing for the Brayton cycle alternator 
is shown in figure 26. The output, which resulted from the data deck shown in figure 26 
is presented in appendix B. 



Number of field coils 2 

Number of turns per coil 376 

Rated -load field temperature, °C 48 
No-load field temperature, °C 30 

Material Copper 

Field coil location code 3 


(a) Field coil design parameters. 



DCOIL 5.584 

BCOIL 1.55 

PT 376 

RD 0.0453 

RT Not applicable to round conductors 

12 48 

T22 30 

RR Need not be read in; program will automatically assume value for copper 

ALPHAR Need not be read in ; program will automatically assume value for copper 

TYPY 3 (see fig. 8(c)) 

TF Field resistance will be printed out at this temperature. Since it is not 

read in, program will assume TF * 25° C 

(b) Correspondence of design parameters to FORTRAN symbols and data card $FIELD. 


Figure 25. - Procedure for preparing data card for NAMELIST name FIELD. 



50 AT A 

AL 4750 STEEL 


1L 3. 

C. 

C.01 

7.74 

0.2 

22.6 

0.4 

51.6 

l.C 

64.5 

2.0 

72. 

3. 

74.5 

4. 

76. 

5. 

79.5 

7. 

90. 

15. 

93.5 

20. 

97. 

3 C . 

ICC. 

5C. 

103. 

200. 




SAE 4 3 4C 

STEEL (RCCP 

TEPP) 






122.5 

C. 

C.01 

6.45 

12. 

12.50 

24. 

19.35 

29. 

58. 

42. 

64.5 

44. 

77.4 

50. 

83.8 

56. 

5C.3 

7C. 

96.8 

92. 

103. 

126. 

109.5 

184. 

116. 

280. 

122.5 

454. 




SAE 1C10 

STEEL ( NCT 

ANNEALED) 






I 2 e. 

C. 

.Cl 

1.0 

1.0 

7.0 

3.0 

11.0 

4 . C 

17.5 

5.C 

46.0 

10. 

72.0 

20. 

83.5 

30. 

51. C 

4 C . 

55.5 

50. 

107. C 

100. 

116. 

200. 

121- 

3CC. 

128. 

500. 




SRAT ING 

VA = 1 4 . 3 , EP 

= 1 2 C » F= 12CC 

, 1 PX = 4 » 

PFC=0 . 75 , 

0.85, 1.0, 

V*IN=0.25, 



G=0. , 0.5, C.75, l.C, 1.5C $ 

SSTATOR C I =3 . 3C , CL=5.364, CL=l-65, LTS=0.C04, SF = 0.9, V*L=22, BK=77.4 S 
SSLOTS Z 7 = 3 , 8C = C.C5» B3=C.23, FY=0.05, H0=0.02, HT=0.C1, HX=0.423, 


HS = 0 .483 , ICC=36 $ 

SWINCNG RF = C * SC=1P, YY = 6, C = 4, OW=O.C201, SN=1, SN1=5, CE=0.25, SD=0.0215, 
PPA=60, T 1 = 1 5 3 , T 1 1 = 6C * 

SAIRGAP GC=C.C2, G2=C.C2, LG2=l-3, WF=1.0 S 

tRGTCft PWl=l.fc, Ph 2= 1 « 756 , PT1=C.215, PT2=L.27, PA2*1.74, 01 = 7.0, 

PA { 1 ) =0. 24 , PA(2)=C.34, PA(3)=0.47, PAt4)=C.65, PA(5)=0.82, PA(6)=0.98, 
PA (111 = 1. 73, CiSH = 2.19, P L = 2 • 4 , SHL1 = 1.65, SHL2=3.65, PAl=0-245, 

C I SH = C • 5 $ 

iYOKE TY=0. 26, 0YC=6.ie, TYL=5.C5, TYR=6.25, TYE=0.25, TY3=0.3, TY2=0.5 * 
5FIEID PCO I L =4 .284 , CCCIL=5.5e4, 8C0IL=1.55, PT = 376, TYPY=3, RD=0.0453, 
T2=48, T 22= 3C * 

SPERPCE P7=12.C $ 

SWINC V I S = 2 . 3E-5 , PRESS=44.4, TEPP=177, P=83.8 $ 

Figure 26. - Listing of data deck for 1200-hertz 120/208 -volt 14.3-kilovolt-ampere 
Brayton alternator. 



APPENDIX B 


TYPICAL COMPUTER PROGRAM OUTPUT LISTING 


The following is the output listing of the Lundell computer program resulting from 
the input data shown in figure 26. While this output listing is typical, the output format 
may vary somewhat, depending for example, on the type of stator slot configuration used 
or on whether or not windage loss was computed by the program. 

Note that on page 54 under ’’No-Load Saturation Data, " the last column contains all 
zeros. In general, when one or more columns contain all zeros, it indicates that the 
flux density calculated for some section of the magnetic circuit was greater than the 
maximum flux density given on the appropriate material data deck, that is, the alternator 
has saturated thereby causing the calculations to be terminated. When this occurs a 
message is printed on the affected page of the output identifying which section of the 

'Q alternator saturated. 

I 

f + ** STATIONARY OUTSIOE-COIL LUNDELL ALTERNATOR *** 

(TWO FIELD COILS) 


ALTERNATOR RATING 


ALTERNATOR KVA 

14.3 

LINE-LINE VOLTAGE 

208. 

LINE-NEUT. VOLTAGE 

120. 

PHASE CURRENT 

39.72 

POWER FACTOR 

0.75 

PHASES 

3 

FREOUENC Y 

12 00. 

POLES 

4 

RPM 

36000.0 

;lqts 


•CONSTANT TOOTH WIDTH 


BO 

0.050 INCHES 

81 

0-151 

BP 

0.154 

B3 

0.230 

BS = (B2+83I/2 

0.192 

HO 

0.020 

HX 

0.423 

HT 

0.010 

HW 

0.020 

HS 

0.483 

NO. OF SLDTS 

36 

SLOT PITCH 

0.288 INCHES 

SLOT PITCH 


AT 1/3 0 I ST. 

0.316 INCHES 


L 1 


1 B 1 1 

i I 

l -BO- 1 

! * * 1 

1 * * 1 H3 

1 * * 1 

1 * * 1 HT 

* * 

HS * * Hrf 

4c ************ * 

*1 * * 1* 

* 1 * * t * HX 

* 1 ************ 1 * 

* l l * 

****************************** 

11 11 

I 1 B2 L l 

II 11 

1 B3 l 

1 1 


45 



AIR GAP 


MAIN AIR GAP 

auxiliary AIR gap 

LENGTH Of AUXILIARY AIR GAP 
EFFECTIVE AIR GAP 


0*020 INCHES 
0.020 
l* 300 
0.021 


CARTER COEFFICIENT 1.074 


ARMATURE WINDING I Y-CQNNEC TED • RANDOM WOUND) 


STRANO DIAMETER 0.0201 INCHES 

DISTANCE ftTWN CL OF STRANDS (RADIAL) 0.0215 


STRANOS/CONOUCTQR IN RADIAL DIR. 
TOTAL ST RA NOS/CO NO UC TOR 
CONDUCTOR AREA 

CURRENT DENSITY AT FULL LOAD 


1 . 

5. 

0.0015 SO- I N. 
6262.41 AMP/ S3- IN. 


COIL EXTENSION BEYOND CORE 
MEAN LENGTH OF 1/2 TURN 
END TURN LENGTH 
StATOR SLOT SKEW 


0.250 INCHES 

5.121 

3.471 

0 . 


RESISTIVITY AT 20 UEG. C 0.6940 MICRO OHM INCHES 

StATOR RESISTANCE AT 25. DEG. C 0.0309 OHMS 


NO. OF EFFECTIVE SERIES TURNS 

TOTAL EFFECTIVE CONDUCTORS 

SLOTS SPANNED 

SLOTS PER POLE PER PHASE 

CONDUCTORS /SLOT 

NO. OF PARALLEL CIRCUITS 

PHASE BELT ANGLE 


22.44 

140.31 

6 . 

3.00 

18. 

4. 

60. DEGREES 


SKEW FACTOR 

DISTR IBUTION FACTOR 

PITCH FACTOR 


1.000 

0.963 

0.866 


FIELD WINDING 


CONDUCTOR DIAMETER 
CONDUCTOR AREA 


0.0453 INCHES 
0.0015 SO- IN. 


NO. OF FIELD COILS 
NO. OF TURNS (PER COIL) 
MEAN LENGTH OF TURN 


2 

376. 

15.501 INCHES 


RESISTIVITY AT 20 DEG. C 
FIELD RESISTANCE AT 25. OEG. C 
(BOTH COILS IN SERIES) 
COIL INSIDE DIAMETER 
COIL OUTSIOE DIAMETER 
COIL WIDTH 


0.6940 MICRO OHM INCHES 
5.1179 OHMS 

4.284 INCHES 

5.584 

1.550 


LOCAT ION CODE (TYPY) 


3 


CONSTANTS 

Cl* FUNOAMENTAL/MAX. OF FIELD FLUX 1.117 


CP. POLE CONSTANT 0.693 
CM* DEMAGNETIZATION FACTOR 0.848 
CO* CROSS MAGNETIZATION FACTOR 0.487 
Dl* POLE FACE LOSS FACTOR 7.000 
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STATOR 


STATOR INSIDE DIAMETER 

3.30 

INCHES 

STATOR OUTSIDE OIAMETER 

5.36 


OVERALL CORE LENGTH 

1.65 


EFFECTIVE CORE LENGTH 

1.48 


DEPTH BELOW SLOT 

0.55 


STACK ING FACTOR 

0.90 


NO. OF COOLING DUCTS 

0. 


WIDTH OF DUCTS 

0. 

INCHES 

CTRE LOSS AT 77.4 KI LOL I NE S/ SO. 1 N. 

22.0 

WATTS^LB 

LAMINATION THICKNESS 

0.004 

IN. 


YOKF 


l 


TYL 


1 


i F Y3 
1 


♦ ***************************** 

* * 

* * 

* ************************ — *- 
* * * * 

* * * * 


* 

* 

* 

* 

* 

* 

* 

*** 

♦ 


******** 

1 1 

y LG2 1 

1 1 


TY 


* * l 

* * 1 

* * 

* * 

* * 

1 * * 

1 * * 

*** * 

T YE * * 

******** 

1 l 1 

1 l TY2 

1 


1 


TY 


0.26 INCHES 


TY2 

TY3 

TYL 

TYF 

TYR 

LG2 


0.50 
0.30 
5.05 
0.25 
6.2 5 
1.30 


YOKE OUTSIDE DIAMETER 


6.18 



ROTOR 


I 



1 



1— 

— 

3.65 

— 1 



1 



1 



1 



1 



1 



I 



I 

************ 




1 

** 

** 

I 


1 

1 

* * 

* * 

1 


1 

* 

* 

* 

* 


1 

* 

* 

* 

* 


1 

****** 

* 

* 

****** - 


l 

* 

* 

* 

* 

L 

1 

* 

* 

* 

* 

l 


* 

* 

* 

* 

l 

3. 2< 

* 

* 

* 

* 



* 

* 

* 

* 

2. 19 

1 

* 

* 

* 

* 


1 

* 

* 

* 

* 

1 

1 

****** 

* 

* 

****** - 


1 

* 

* 

* 

* 


1 

* 

* 

♦ 

* 


1 


* * 

* * 



1 


** 

** 



l 


************ 

— 




1 1 

i 1 

1 1.65 1 

1 i 


POLE DIMENSIONS (INCHES) 


NARROW END 

WIDTH 1.60 

THICKNESS 0.21 

AREA 0.24 

WIDE END 

WIDTH 1.76 

THICKNESS 1.27 

AREA 1.74 

POLE LENGTH 2.40 

INSIDE SHAFT DIAMETER 

(FOR HOLLOW SHAFT) 0.500 

POLE EMBRACE 0.69 


PERIPHERAL SPEED 30748.32 FEE T/MI N» 

(AT MAIN AIR GAP) 

SPEC. TANGENTIAL FORCE 0.99 LBS/SQ.IN. 

DENSITY OF NQN-MAGNETIC 

ROTOR MATERIAL 0.283 LBS/CU. IN. 



WEIGHTS 


STATOR COND. 

F I ELD CGND . 

STATOR IRON 

ROTOR 

YOKE 

TOTAL 

( ELECTROMAGNETIC! 


WINDAGE 


FLUID PROPERTIES 
VISCOSITY 
DENSITY 
PRESSURE 
TEMPERATURE 
MOLECULAR WEIGHT 

REYNOLDS NUMBER 
MA IN AIR GAP 
AUX- AIR GAP 
WINDAGE LOSS ( WATTS! 
MAIN AIR GAP 
AUX. AIR GAP 
TOTAL 


1.689 POUNDS 
6. 031 
4. 584 
9.662 
11.287 


33.252 


2.300E-05 LBS/SEC-FT 
4.2664 LBS/CU.FT. 
44.40 LBS/SQ. IN 
177.00 DEG. C 
83.8 


15831.50 

10635.27 

652.43 

114.39 

881.22 


PERMEANCES 

WINDING LEAKAGE 

(PER INCH OF CORE LENGTH) 

STATOR SLOT 5.903 LINES/AMPERE TURN 


STATOR 

END 

8.009 

SPECIAL (SPS) 

0. 

PI (POLE 

TIP! 

1.042 

P 2 


6.218 

P3 


2.431 

P4 


3.224 

P 5 


20.218 

P6 


0.412 

P 7 


12.000 

Pft (MAIN 

AIR GAP! 

2540.897 

P9 (AUX. 

AIR GAP) 

1439.618 

PIO 


23. 748 


REACTANCES 


AMPERE CONDUCTORS/INCH 
REACTANCE FACTOR 

STATOR WINDING LEAKAGE 
ARM. REACTION (DIRECT) 

ARM. REACTION (QUAD.) 
SYNCHRONOUS (DIRECT) 

S YNC HRONOUS (QUAD. ) 

FIELD LEAKAGE 
TRANSIENT 

FIELD SELF INDUCTANCE 

OPEN CIRCUIT TIME CONSTANT 
(FIELD ONLY, AT 25. DEG. 

SHORT CIRCUIT AMPERE-TURNS 
SHORT CIRCUIT RATIO 


537.518 

0.993 

13.775 PERCENT 
142.623 
73.318 
155.398 
87.392 
<►51.382 
122.182 

2.355 HENRIES 


C) 0.46042 SECONDS 

1180.916 

0.922 


STATOR MATERIAL - AL 4750 STEEL 

ROTOR MATERIAL — SAE 4340 STEEL (ROOM TEMP) 

YOKE MATERIAL SAE 1010 STEEL (NOT ANNEALED) 


MAGNETIZATION CHAR AC Tt R I ST I C S (NO LOAD, RATED VOLTAGE) 

TG* TOTAL USEFUL FLUX 564.33 KILOLINES 

FQ, MAIN AIR-GAP FLUX/POLE 97.32 

PP Li, FLUX LEAVING BASE OF POLE 109.21 

PPL 2, FLUX IN TAPE RE 3 ROTOR SECT. 222.15 

PPL3, FLUX INTO ROTOR END SECT. 222. 4L 

PPL4, (LEAKAGE THROUGH P5)*PPL3 237.37 

FLUX DENSITIES ( M LOL I NE S/SQ IN) 


MAIN AIR GAP 

32.990 

AUXILIARY AIR GAP 

24.857 

STATOR CORE 

60.245 

STATOR TEETH 

75.339 

YOKE - CYLINDRICAL SECTION 

49.193 

END SECTION 

53.720 

ROTOR - POLE 

62.690 

END SECTION 

59.177 

TAPERED SECTION 

53.755 

PER E-TURNS 

MAIN AIR GAP 

223.342 

AUXILIARY AIR GAP 

154.493 

STATOR CORE 

5.789 

STATOR TEETH 

2.393 

YOKE - CYLINDRICAL SECTION 

29.127 

END SECTION 

15.300 

RJTUR - POLE 

65.638 

END SECTION 

38.375 

TAPERED SECTION 

43.537 

TO TAL 

L039.1L5 
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ALTERNATOR 

LOAD CHARACTERISTICS 

(RATED VOLTAGE, 

0.75 ? 3HET FACT3TI 



PERCENT KVA 

0 . 


53. 

75. 

100 . 

153. 

FLUX DENSITIES (KL/SQ-IN) 

TEETH 

75.04 


83.45 

83.18 

86. 22 

3 3. 4 V 

CORE 

60.25 


64.59 

66.78 

69.22 

75. n 

YOKE - CYLINDRICAL SECTION 

49.19 


59.68 

65.48 

71.50 

rr 

m 

CO 

ENO SECTION 

53.72 


65.18 

71.50 

78.19 

92.73 

ROTOR - POLE 

62.69 


73.24 

78.98 

85.32 

97. 95 

END SECTION 

59. 18 


59.98 

75.89 

82.11 

95.42 

TAPERED SEC. 

63.76 


75.35 

81.69 

88. 37 

132. 54 

AMPERE-TURNS 

MAIN AIR GAP 

223.04 


239.14 

247.25 

256.28 

277.7? 

AUX. AIR GAP 

154.49 


182.70 

198.13 

214.37 

2 49. L 3 

TEETH 

2.09 


3.62 

4.42 

5.51 

9.51 

CORE 

5.79 


7.31 

8.23 

9.3 9 

15.71 

YOKE - CYLINDRICAL SECTION 

29.13 


33.53 

44.96 

52.94 

3 4. 53 

END SECTION 

15.00 


19.55 

22.67 

26.74 

+ 3.75 

ROTOR - POLE 

65.64 


71.26 

74.89 

80. 52 

113.33 

END SECTION 

38.07 


41.19 

43.02 

45.22 

53. ?3 

TAPERED SEC. 

43.54 


48.58 

53.06 

63.28 

i L S . 54 

TOTAL AMPERE TURNS 

1089.11 


1653.88 

1982.35 

2342. 84 

323?. 85 

FIELD CURRENT IAMPS) 

1.45 


2.20 

2.64 

3. 12 

4. 33 

CURRENT DENS. (FIELD) 

898.60 


1364.58 

1635.59 

1933.02 

2557.35 

FIELD VOLTS 

7.56 


11.75 

14.35 

17.36 

25. 71 

TEMPERATURES IDEG.CI 

FIELD 

30.00 


36.48 

41.43 

48.30 

59.75 

ARMATURE 

60.00 


83.25 

112.31 

153.00 

2 >9. ?5 

RESISTANCES (OHMS) 

FIELD 

5.22 


5.34 

5.44 

5.57 

5.93 

ARMATURE 

0.0350 


D.D378 

0.0412 

3.0461 

0. 0599 

ALTERNATOR LOSSES (WATTS) 

FIELD 

10.94 


25.85 

37.82 

54. 08 

113.54 

WINDAGE 

881.22 


881.22 

881.22 

881.22 

931.2? 

STATOR TOOTH 

64.50 


74.14 

79.25 

85.15 

1 33. 33 

STATOR CORE 

139.70 


163.59 

171. 66 

184.43 

215.59 

POLE FACE 

958.70 


1224.29 

1556.27 

2021.35 

3343.99 

STATOR COPPER 

0 . 


44.71 

109.79 

218.39 

537.92 

EDDY 

0 . 


3. 

0 . 

0 . 

3. 

MISC. LOAD 

0 . 


71.50 

107.25 

143.30 

214.53 

TOTAL 

2055.06 


2432.29 

2943.27 

3587. 02 

5539. 75 

ALTERNATOR OUTPUT (KVA) 

0 . 


7.15 

10.72 

14. 30 

21.45 

ALTERNATOR OUTPUT (KW) 

0 . 


5.36 

8.04 

10.72 

15.39 

EFFICIENCY (PER CENT) 

ELECTRO-MAGNETIC 

0 . 


77.01 

79.60 

79.85 

77.53 

OVER-ALL 

0 . 


63.36 

73.21 

74.94 

7 4. 49 
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ALTERNATOR 

LOAD CHARACTERISTICS 

(RATED VOLTAGE, 0.85 

ER FACTOR) 



PERCENT KVA 

0 . 

53. 

75. 

100. 

153. 

FLUX OENSITIES (KL/SQ-IN) 

TEETH 

75.0* 

73.93 

80.8* 

83. 39 

37.93 

CORE 

60.25 

63.37 

6*. 93 

66. 71 

7 3.53 

YOKE - CYLINDRICAL SECTION 

*9.19 

58.33 

63. 19 

68.76 

33.93 

END SECTION 

53.72 

63.37 

69.03 

75.39 

3 3 . V+ 

ROTOR - POLE 

62.69 

71. *3 

76.37 

81.77 

93. 67 

ENO SECTION 

59.18 

53. 17 

73.3* 

78.9* 

91.19 

TAPERED SEC. 

63.76 

73. *0 

78.95 

8*. 95 

99.39 

AMPERE-TURNS 

MAIN AIR GAP 

223.0* 

23*. 61 

2*0.29 

2*6.96 

261.37 

AUK. AIR GAP 

15*. *9 

177.96 

191.46 

206.09 

233.37 

TEETH 

2.09 

3.20 

3.73 

*.39 

5. 24 

CORE 

5.79 

5.85 

7.** 

8.20 

13.12 

YOKE - CYLINDRICAL SECTION 

29. 13 

35.87 

*2.30 

*9.38 

7 3.43 

END SECTION 

15.00 

18.75 

21.33 

2*. 73 

35. 46 

ROTOR - POLE 

65.6* 

73.22 

73.15 

77.2* 

96. 94 

END SECTION 

38.07 

*3.5* 

*2.2* 

**.3* 

49.95 

TAPERED SEC. 

*3.5* 

*7.56 

50.58 

56.52 

93.32 

TOTAL AMPERE TURNS 

1089.11 

1585.72 

1891.78 

2226.48 

3315.73 

FIELD CURRENT (AMPS) 

1 . *5 

2.11 

2.52 

2.96 

4.31 

CURRENT DENS. (FIELD) 

898.60 

1333.3* 

1560.87 

1837.32 

2433. 25 

FIELD VOLTS 

7.56 

11.23 

13.62 

16.37 

23.54 

TEMPERATURES ( DEG.C ) 

FIELD 

30.00 

35.56 

*0.01 

*5.78 

63. 39 

ARMATURE 

60.00 

33.25 

112.31 

153.30 

259.25 

RESISTANCES (OHMS ) 

FIELD 

5.22 

5.33 

5. *1 

5.53 

5. B 7 

ARMATURE 

0.0350 

3.3378 

0.0*12 

0.0*61 

0. 05 99 

ALTERNATOR LOSSES (WATTS) 

FIELD 

10.9* 

23.68 

34.26 

*8. *6 

9 4.39 

WINDAGE 

881.22 

881.22 

881.22 

881.22 

331.22 

STATOR TOOTH 

6*. 50 

71.36 

7*. 86 

79.37 

33.57 

STATOR CORE 

139.70 

15*. 56 

162. 1* 

171.27 

191.94 

POLE FACE 

958.70 

122*. 29 

1556.27 

2021. 35 

3343.99 

STATOR COPPER 

0 . 

**•71 

109.79 

218.39 

637. 92 

EDDY 

0 . 

3. 

0 . 

0 . 

0 . 

MISC. LOAD 

0 . 

71.50 

107.25 

1*3. 00 

2L4.53 

TOTAL 

2055.06 

2*71.32 

2925.79 

3562.15 

5457.44 

ALTERNATOR OUTPUT (KVA) 

0 . 

7.15 

10. 72 

14.30 

21.45 

ALTERNATOR OUTPUT (KW) 

0 . 

5.38 

9. 12 

12.15 

13.23 

EFFICIENCY (PER CENT) 

ELECTRO-MAGNE TIC 

0 . 

79.26 

81.68 

81.93 

7 9.94 

OVER-ALL 0. 71.09 

< STARTING WITH COLUMN 5 ALL CALCS ASSUME UNIFORM AIR-GAP 

75.70 

F„ JX DISTRIBUTION) 

77.3* 

76. 95 


ALTERNATOR 

LOAD CHARACTERISTICS 

(RATED VOLTAGE, 

1. 00 ? 3WE* FACTO*) 



PERCENT KVA 

0 . 

53. 

75. 

100. 

1 3 0. 

FLUX DENSITIES (KL/SQ-INI 

TEETH 

75.0* 

83.1* 

78.0* 

76. *3 

75.3b 

CORE 

60.25 

6*. 3* 

62.66 

61.36 

60.90 

YOKE - CYLINDRICAL SECTION 

*9.19 

55.91 

58.00 

61.01 

69. ** 

END SECTION 

53.72 

61.06 

63.3* 

66.52 

75.33 

ROTOR - POLE 

62.69 

69.89 

71. 13 

73.29 

30.43 

END SECTION 

59.18 

66.39 

67.98 

70.51 

73.30 

TAPERED SEC. 

63.76 

71.50 

73.20 

75.90 

3*. 2* 

AMPERE-TURNS 

MAIN AIR GAP 

223.0* 

233.20 

231.97 

227. 18 

2 2 5. 4 3 

AUX. AIR GAP 

15*. *9 

173.32 

177. *8 

18*. 09 

?3V. 4? 

TEETH 

2.09 

3.5* 

2.9* 

2.52 

2.36 

CORE 

5.79 

7.21 

6.59 

6.15 

5. 30 

YOKE - CYLINDRICAL SECTION 

29.13 

3*. 8* 

36.8* 

39.91 

49.96 

END SECTION 

15.00 

17.76 

18.73 

20.22 

? 5 . 1 3 

ROTOR - POLE 

65.6* 

69.31 

70.08 

71. *0 

76.27 

END SECTION 

38.07 

*3.11 

*0.58 

*1.3 5 

43.31 

TAPERED SEC. 

*3.5* 

*6.78 

*7.56 

*8.8* 

55.54 

TOTAL AMPERE TURNS 

1089. 11 

1396.77 

1612.30 

1877.58 

2505.29 

F1ELO CURRENT (AMPS! 

1 . *5 

1.86 

2.1* 

2.50 

3.33 

CURRFVT DENS. ( FIELD J 

898.60 

l L52 • ** 

1330.27 

15*9.15 

2067. 36 

FIELD VOLTS 

7.56 

9.81 

It. *3 

13.51 

18.73 

TEMPERATURES (DEG.C) 

FIELD 

30.00 

33.20 

35.91 

39.79 

51.29 

ARMATURE 

60.00 

83.25 

112.31 

153.00 

269. 25 

RESISTANCES (OHMS) 

FIELD 

5.22 

5.28 

5.33 

5. * 1 

5.64 

ARMATURE 

0.0350 

3.3378 

0.0*12 

3.0*61 

3. 3599 

ALTERNATOR LOSSES (WATTS) 

FIELD 

10.9* 

19.21 

24.52 

33.72 

62.56 

WINDAGE 

881.22 

831.22 

881.22 

881.22 

381.2? 

STATOR TOOTH 

6*. 50 

73.56 

69.76 

66.91 

55.91 

STATOR CORE 

139.70 

159.33 

151.11 

14*. 93 

142. 77 

POLE FACE 

958.70 

122*. 29 

1556.27 

2021.05 

3343.99 

STATOR COPPER 

0 . 

**• 71 

109.79 

218.09 

637.92 

EDDY 

0 . 

3. 

0 . 

0 . 

3. 

MISC. LOAD 

0 . 

71.50 

107.25 

1*3.00 

2L4. 50 

TOTAL 

2055.06 

2*72.81 

2899.92 

3508.93 

5353. 33 

ALT FRNATOR OUTPUT (KVA) 

0 * 

7.15 

10.72 

14.30 

21.45 

ALTERNATOR OUTPUT (KW) 

0 . 

7.15 

10.72 

14.30 

21.45 

EFFICIENCY (PER CENT) 

ELECTRO-MAGNETIC 

0 . 

81.79 

8*. 15 

84.48 

32.75 

OVER-ALL 

0 . 

7*. 30 

78.72 

80.30 

30. 33 



cn 


NU-LOAD SATURATION OAT A 


VOL TAGE 

PERCENT 

25.00 

40. 00 

55.00 

70.00 

85.00 

L IN E-NEUTRAL 

30. UO 

48.00 

66.00 

84.00 

102.00 

L RE-LIME 

5L.96 

83. L4 

1 14.32 

145.49 

176.67 

FIELD CURRENT 

0.47 

0.67 

0.86 

1.05 

1.25 

FLUX PER PULb 

24.45 

39.13 

53.80 

68.47 

83.14 

FLUX DEMSITIES 

CO*E 

15. 35 

24.21 

33.19 

42.23 

51.24 

TEETH 

YOKE 

L9.L2 

30.16 

4L.34 

52.56 

63.82 

END SECTION 

13.62 

21.65 

29.63 

37.63 

45.64 

CYL. SECT. 
ROTOR 

12.47 

19.82 

27.14 

34.46 

41.83 

POLE 

15.64 

25.03 

3 4.42 

43.82 

53.23 

END SECTION 

14. 78 

23.64 

32.49 

41.36 

50. 24 

TAPERED SEC. 

15.92 

25.46 

35.09 

44.56 

54.12 

AMPERE-TURNS 

MAIM AIR-GAP 

56.82 

89.64 

122.83 

156.24 

189.70 

AUX AIR-GAP 

38.59 

61. 71 

34.33 

107.97 

131.16 

CORE 

1.04 

1.53 

2.03 

2.70 

3.60 

TEETH 

YOKE 

0.16 

0.25 

0.35 

0.51 

0.93 

END SECTIUN 

6. 03 

7.57 

8.95 

10.61 

12.53 

CYL. SECT. 
ROTOR 

11.25 

L 4 . 15 

15.91 

20.20 

24.15 

POLE 

36.45 

46.41 

51.26 

56.06 

60.79 

END SECTION 

15.74 

23.97 

28.52 

32.57 

35.57 

TAPERED SEC. 

26.08 

30.66 

33.55 

36.74 

40.23 

TOTAL 

349. 80 

506.70 

STATOR 

643.45 792.34 

TOOTH IS SATJRATEO 

937.63 


100.00 

115.00 

130.00 

L35.30 

3. 

120 .00 

I 38.00 

156.00 

162.00 

0. 

207.85 

239.02 

270.20 

233.59 

3. 

1.45 

1.67 

1.95 

2.13 

0. 

97 .82 

112.49 

127.15 

132.05 

3. 


60.25 

69.28 

78. 20 

81 . .33 

3. 

75.04 

86.29 

97.43 

133.99 

3. 

53.72 

61 . 94 

70.73 

74.43 

3. 

49.19 

56. 72 

54. 77 

5 3.15 

9. 

62.69 

72.26 

82.27 

8 5.1 9 

3. 

59.18 

68.23 

77.75 

8! . 57 

3. 

63.76 

73.51 

83.77 

87.37 

3. 

223.04 

256.48 

269.52 

333.17 

3. 

154.49 

178.13 

203. 01 

>12.93 

). 

5.79 

9.42 

22.48 

2 8.55 

3. 

2.09 

5.53 

15.52 

38. 1 1 

3. 

15.00 

18.13 

22.25 

24.33 

3. 

29.13 

35. 60 

44. L? 

43.29 

3. 

65 . 64 

70.52 

77.21 

91.55 

3. 

38.07 

40.66 

43. 61 

45. 3? 

3. 

43.54 

47. 71 

55.03 

52.21 

3. 

1089.11 

1255.27 

1470.22 

1503.25 

3. 



APPENDIX C 


COMPLETE FORTRAN LISTINGS AND FLOW CHARTS OF 
LUNDELL ALTERNATOR ANALYSIS PROGRAM 

The complete FORTRAN listings of the main program and the four subroutines, 
which together constitute the Lundell alternator analysis program, are contained here- 
in. The main program is LUNDEL and the four subroutines are, in the order given 
SUBLUN, OUTPUT, MAGNET, and WINDGE. Each program listing, is followed by its 
flow chart. 


c 

c 

c 


c 

c 

c 


c 


c 

c 


1 

2 

3 


MAIN PROGRAM LUNDEL 

FOR use WITH STATIONARY. OUTSIDE COIL LUNDELL ALTERNATORS 

COMMON /COM 2/ AC.AS,a3.6K.aG.3S,C.Cl,Dl.DI.DlSH,DISHl,0J,DYC,EE.Ef> 
l. F. GC. HC. HM.NFC ,PA l.PE .PF .Pi . P L, P F , RPM. RR . RS , SC , SO ,S H. S HL2 . SN, SS.T 
2S.TT.TY.TYL, WF, WL.XA.XB.XD.XL 

COMMON /CUM 3/ ACOR E . A I . A SHF T . A TOO TH , A YKl *BAG2 • BC L, BPL . B ROT , BS HFT , 8 
IT L . B YCL X. BYCL2,QCQRE.DSHMIQ,DYKI,0 YK2 ,FC L.FF L, FGL1 . FGL 2 • F 3XL » FPL • F 
2R OT. FS HFT. FTL. F YCL l .FYCL2.KSA T.PAA. I TRTN 

COMMON /COM4/ AB, 4LPHAR , ALPHA S . BG . FGML .FH ,F 0 . G . G A . PA ,PFC . QQ, RY , T 1, 
XT 11. T2, T2 2. TG. VMIN.ZG.ZZ 

COMMON /CGM5 / A A . AG2 .CL. CP .D Y2 .DY3 .HS.LG2.PX .P10.P5tP6.P7.P8tP9.PA 
12. PPL. PPL 1.PPL2.PPL3.PPL4.PX.TY2 . TY3 .TYPY 

INTEGER ZZ 

REAL LG2 

DIMENSION POLL ( X 0) . UPERVUO), UVLL(iO). OVLN(IO). FGLLillOJ, BPLL 
LI 10), FFLLI10J. BCLL(XO). BTLL(IO). BYCLLl(iO). BYCLL2<1J). FYCLL1 
21 10), FYCLL2UG), FPLLUO), FTLLUO), FCLL(iO). FGLL2(10), FK10), 
3 BSHFTL < 10 ) « FSHFTLI LO) • BROTL(iQ), FROTLUO). TTB ( 5 ) • TFA15). RRA 
4(5). R RBI 5 ) • E Z ( 5 ) • WMIS(5), CUD ( 5 ) , EF ( 5) . PR I 5 ) , S T ( 5 ) * A OL ( 5 ) . 
5P P( 5 ) • PSI 5) , £ X( 5 ) . SP I 5) . AKVA(5), *A(5>, EX(5), E(5). AIOO), G 
615). YMAT ( 6 ) , RMAT ( 6 ) * SMAT(6). PFC(5), PAiil). PAAC 10) . CMPNTI28) 

DATA CMPNTI 1I/24H ROTOR END SEC TI ON/ .CMPNT ( 5 ) /24H TAPERED 

XROTOK SECTION/. CMPNT(9)/24H STATOR TOOT H/ . CMPNT ( 13 ) /24H 

2 STATOR CORE / « CMPNT ( 17 ) /24H POLE/.C 

3MPNTI21)/24H YOKE END S EC T I UN/ , CMPNT 1 25 ) /24HC YU NDR 1 CAL YG< 

4E SECTION/ 


READ I 5. 1 ) SMAT 

R FAD 15,2) (All I). 1 = 1* 29) 

READ < 5.1 ) RMAT 

READ 15.2) I AH 11*1=31.59) 

READ I 5. XI YMAT 

READ (5.2) ( AI I I ) . i =61 . 891 

FORMAT I6A6) 

FORMAT ( 8F 10-1 1 
CALL SUBLUN 
CALL OUTPUT 


A 1 
A 2 
A 3 
A 4 
A 5 
A 6 
A 7 
A 8 
A 9 
A 10 
A 11 
A 12 
A 13 
A 14 
A 15 
A 16 
A 17 
A 18 
A 19 
A 20 
A 21 
A 22 
A 23 
A 24 
A 25 
A 26 
A 27 
A 28 
A 29 
A 30 
A 31 
A 32 
A 33 
A 34 
A 35 
A 36 
A 37 
A 38 
A 39 
A 40 
A 41 
A 42 
A 43 


1 

3 

10 

12 

19 

21 


29 
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COMPUTE TOOH WIDTH AT L/3 DISTANCE FkOM NARROWEST SECTION 


A 

A 

A 


44 

45 

46 


C 

C 

C 


IF ill-3) 4,5,6 

4 SM=TT-BS 
GO TO 8 

5 SM = ( 3 * 141 6* ( 3 U2.*HS) /OD)~B3 
GO TO 8 

6 IF (ZZ-4) 5, 7,4 

7 SM=TT- ,94*b S 

8 CONTINUE 
C 

C AREAS AN 3 LENGTHS FOR MAGNETIC CALCULATiUNS 

C 

AYK 1=3.L416*(UYC-TY)*TY 
DYK 1= ( TYL +TY31/2. 

D YK 2= ( QYG-2.* TY-DY2) /6. 

AT3UTH=Q3*SS* SM*Pt /PX 
ACORE^SS* 1 ( I DU- { D l +2 .*M S) ) /2 • ) 

DC3RE=< /3, ) «3. 1U5/PX 

AG2=L G 2* 3 ♦ 141 6* D l SHI 

ASHFT = 0.7854- { D I SH i ** 2 -0 1 SH**2) 

D SHM I D =( SHL2-CL ) /2. 
r=i 

9 1 = 1 + 1 

IF (?A(I) .LE.O.) GO TO 10 
IF (I * LT ♦ 9 ) GO TO 9 
GO TO 14 

10 J = I 

11 J=J + l 

IF { P A ( J 1 .GT.O. ) GO TO 12 
GO TO 11 

12 JU=J-1 

DO 13 K = I , JJ 

13 PAnO-K’AIJI-PMI-l)) /FLOAT l J- I +1 ) ) *F LOAT ( K- I + l ) * PAl I -l ) 
r=j 

GO TO 5 

14 DO 15 1=1,10 

15 paa( i )=ipa(i ) +pai inn / 2 . 

C 

C INITIALIZE SUBSCRIPTED VARIABLES USED IN LOAD CHARACT. CA_ CS 

C 

DO 16 J=1 , 275 

16 FGLLi<J) = 0* 

C 

C 

C NO-LOAD, RATED VOLTAGE MAGNE TI Z ATI UN CHARACT c RI ST I CS 

c 

C 

FGXL = 0 . 

P PL =F J 
ITR TN = 10 
CALL MAGNET 
J = 1 

FYCLL 1C J) =FY:L1 
FYCLL 2(J ) -FYCL2 
FPL L ( J )=F PL 
FTLL( J )=F TL 
FCLL ( J )=FCL 


A 47 
A 48 
A 49 
A SO 
A 51 
A 52 
A 53 
A 54 
A 55 
A 56 
A 57 
A 58 
A 59 
A 60 
A 61 
A 62 
A 63 
A 64 
A 65 
A 66 
A 67 
A 68 
A 69 
A 70 
A 71 
A 72 
A 73 
A 74 
A 75 
A 76 
A 77 
A 7b 
A 79 
A 80 
A 81 
A 82 
A 83 
A 84 
A 85 
A 86 
A 87 
A 88 
A 89 
A 90 
A 91 
A 92 
A 93 
A 94 
A 95 
A 96 
A 97 
A 98 
A 99 
A 100 
A 1 0 L 
A 102 
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FFLL(J)=FFL A 103 

FSHFTL ( J)=FSHFT A 104 

BYCLLK J)=BYCL1 A 105 

BYCLL2I J)=BYCL2 A 106 

B PLL ( J ) =BPL A 107 

8 ILL ( J ) =B TL A 108 

B CLL ( J ) =BCL A 109 

BSHFTLI J)=BSHFT A 110 

FGLL l( J ) = FGL 1 A 111 

BROTLI J )=BROT A 112 

F ROTL ( J ) = FRO 7 A 113 

F GLL2( J ) = FGL 2 A 114 

C A 115 

C SHORT CIRCUIT RATIO AND SHORT CIRCUIT AMPERE-TURNS CALCS A 116 

C A 117 

FSC=XA*(FGL1+FGL2)*0.02 A 118 

SCR=FFL/FSC A 119 

WRITE 16,17) F SC , SCR A 120 

17 FORMAT <IHK,9X,27H SHORT CIRCUIT A MPE RE-TURNS , FI S . 3/ 13 K , 2 3 H SHORT A 121 

1CIRCUIT RATIO, F 23, 3 ) A 122 

WRITE ( 6, 18) SMAT A 123 

18 FORMAT (1H1.18H STATOR MATERIAL -,1H ,6A6) A 124 

WRITE I 6, 19) RMAT A 125 

19 FORMAT ( 1 HK , 18H ROTOR MATERIAL --,1H ,6A6) A 126 

WRITE 16.20) YMAT A 127 

20 FORMAT I1HK.18H YOKE MATERIAL ,IH ,6A6) A 128 

WRITE 16,21) TG.F0.PPL1,PPL2.PPL3,PPL4 A 129 

21 FORMAT (1HL.30H MAGNETIZATION CHARACTERI STI CS, 25H (NO LOAD, RATED A 130 

1V0LTAGE)//10X,22H TG , TOTAL USEFUL F LUX ,F20. 2 , 1 0 H KI LOL I NE S / 10X , 27 A 131 

2H FO. MAIN AIR-GAP FLUX/POLE ,F 1 5. 2/1 0 X , 32H PPL1, FLUX LEAVING BASE A 132 

3 OF POLE, F10.2/10X.34H PPL2 , FLUX IN TAPERED ROTOR SECT . , F 8 . 2/ l OX, A 133 

43 2H PPL 3. FLUX INTO ROTOR END SEC T. , F 1 0 . 2/ 13 X , 32H PPLA, ( LEAKAGE T A 134 

5HR0UGH P5) +PPL 3 . F 1 0. 2/ / 1 OX, 33 M FLUX DENSITIES ( K IL OL I NES / S Q IN)) A 135 

WRITE (6,22) BG,BAG2,BCL,BTL,3YCLl »B YCL2 .BPL.BSHFT , BROT A 136 

2? FORMAT I13X.13H MAIN AIR GAP , = 27. 3/13X . I 8H AUXILIARY AIR 3AP.F22.3 A 137 

1/13X.12H STATOR CORE ,F 28. 3/ 13 X , 1 3H STATOR TEETH, F27 . 3/ 1 3X . 2 7H YOKE A 138 

2 - CYLINDRICAL SEC T I ON . F 1 3. 3 / 1 8X , 1 4H END SEC T I ON , F21 . 3/ 1 3X , 1 3H R A 139 

30 TOR - P0LE.F27.3/21X.12H END SECT I ON, F20. 3/2 l X, 1 6H TAPERED SECTIO A 140 

4N.F16.3) A 141 

WRITE (6,23) A 142 

23 FORMAT (1HK.9X.13H AMPERE-TURNS) A 143 

WRITE (6,22) FGL1.FGL2.FCL.FTL.FYCL1 .FYCL2 .FPL ,FSHFT , FROT A 144 

WRITE (6.24) FFL A 145 

24 FORMAT (1H . 12X.6H TOTAL ,F34. 3 ///) A 146 

IF (ITRTN.EQ.O) WRITE (6.25) A 147 

25 FORMAT (1H , 20X , 54H ( ABOVE CALCS ASSUME UNIFORM AIR-GAP FLJX DISTR1 A 148 

1BUT ION ) ) A 149 

IF (KSAT.E0.10) GO TO 27 A 150 

J=4*KSAT A 151 

I =J— 3 A 152 

WRITE (6,26) (CMPNT(K) »K=I , J) A 153 

26 FORMAT (20X.4A6.36H SATURATED AT NO-LOAD. RATED VOLTAGE) A 154 

GO TO 3 A 155 
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A 156 
A 157 
A 158 


C 

c WO- LOAD POLE-FACE LOSS CALCULATION 

C 


27 

GT= 80/ GC 

A 

159 


AA = 1.7 5/< GF**l . 351*0.8 

A 

160 


GF = AA*P I* SC/ (C*FH) 

A 

161 


D 2-= BG**2. 5* 0*00 006 1 

A 

162 


0 3=( 0. 0 16 7*0G*RPM) **1.65*1. 5I4 7E-5 

A 

163 


IF ITS-0.9) 28,28,29 

A 

164 

28 

D4=FS**L. 285*0. 81 

A 

165 


GO TO 32 

A 

166 

29 

IF I FS-1. 51 30,30,31 

A 

167 

30 

f)4=TS**i. 145*0.79 

A 

168 


GO TO 32 

A 

169 

31 

D 4=T S** 0. 79*0.92 

A 

170 

32 

IF (GT-L. 7) 33.33,34 

A 

171 

33 

05=GT**2. 31*0.3 

A 

172 


GO TO 39 

A 

173 

34 

IF (GT-3.0) 35,35,36 

A 

174 

35 

D 5=GT **2.0*0.35 

A 

175 


GO TO 39 

A 

176 

36 

IF < GT-5.0) 37,37,38 

A 

177 

37 

D 5=GT** 1 . 4* 0.67 6 

A 

178 


GO TO 39 

A 

179 

38 

05=GT**0. 965*1. 38 

A 

180 

39 

06=10.0** i 0.93 2*01-1. 606) 

A 

181 


N=0 1*02* 03* 04*0 5* 06* G A 

A 

182 

C 


A 

183 

C 

CALCULATE NO-L 0 AO, R ATEO VOLTAGE TOOTH AND CORE LOSS 

A 

184 

c 


A 

185 


WT=( SM)*OQ*SS*HS*0.849*IBTLLI 1 ) /BK) **2. 0*rfL 

A 

186 


M 0= ( DU- EC > * 2 .67*HC * S S* I B C L L (11 /B K ) * *2 . 0* WL 

A 

187 

c 


A 

188 

c 

ARRANGING LOAO POINTS IN INCREASING ORDER 

A 

189 

c 


A 

190 


DO 41 J = 1 , 4 

A 

191 


I A=5-J 

A 

192 


DO 41 1*1* IA 

A 

193 


IF (G< n.GT.Gf I + Dl GO TO 40 

A 

194 


GO TO 41 

A 

195 

40 

POL=G< I ) 

A 

196 


G( I )=G< 1*1) 

A 

197 


G I 1*1) =POL 

A 

198 

41 

CONTINUE 

A 

199 


G( 11=0. 

A 

200 

C 


A 

201 


M M= 5 

A 

202 


DO 42 1=2, 5 

A 

203 

42 

IF (G( I ).GE. 1.0. AND.S( 1-1) .LT. 0.999) MM=1 

A 

204 

C 


A 

205 

c 

DEFINING NPF 

A 

206 

c 


A 

207 


DO 43 1 = 1, 5 

A 

208 

43 

IF (PPC< I l.LT.0.001) GO TO 44 

A 

209 


1=1*1 

A 

210 

44 

NPF=l-l 

A 

211 
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C A 212 

c A 213 

C CALCULATE ALTERNATOR LOAD CHARACTERISTICS A 2 14 

c A 215 

C A 216 

1=0 A 217 

45 1=1+1 A 218 

K =0 A 219 

l TR TN = 10 A 220 

I A= 10 A 221 

PF=PFC ( I ) A 222 

CK= 1.0 A 223 

IF IPF.GE.0.95) CK-l.lO A 224 

AM = ARCOS ( PF ) A 225 

EZ( i) = l.O A 226 

J =0 A 227 

J A= 5 A 228 

C A 2 29 

C ARMATURE TEMPERATURE AND RESISTANCE CALCULATIUN A 230 

C A 231 

46 J=J+1 A 232 

TTAiJ ) =( I 1-T11)*3(J)*G(J)+T11 A 233 

RB=( 1.0E-o)*RS*( 1.0+ALPHAS*! TTA( J)-20. ) ) A 234 

R R A ( J ) =R 3 * R Y A 235 

IF ( J •EQ. 1) GO TO 46 A 236 

C A 237 

C EDDY FACTOR CALCULATIONS A 238 

C A 239 

IF (SH) 47,47,48 A 240 

47 EZ(J) = 1. A 241 

GO TO 49 A 242 

48 AA=0.584+(SM*SN-1.0)*0.0625*(SD*CL/($H*HM/2.>)**2 A 243 

AB=(SH*SC*F*AC/(dS*R3*1.0E6))**2.0 A 244 

ET=AA*AB*0. 00335+1. 0 A 245 

EB=ET-Q.00168*A6 A 246 

EZ( J ) =(ET+tB ) *0 . d A 247 

C A 248 

49 X I = ATANM ( ( J )) / luU. ) tSI N<4N)) / ( PF+RKAI J) * ( P( / £P)*G( J I ) ) A 249 

8B=XI-AN A 250 

GXX = ( C 0 S( BB )+(RRA( J) * l Pi/EP)*G( J)*COS (XI )) +< ( X L + 0. 0 7*X 0 ) *S IN( XI )*G A 251 

1< J )*0 .01) )*;X a 252 

FGXL=FGML*G( J)*SIN(XI )*2. A 253 

PPL=FO*GXX A 254 

CALL MAGNET A 255 

IF (ITRTM.EJ.Q) I TR TM = J A 256 

FYCLL 1( J) =F YC L 1 A 257 

FYCLL 2(J) =FYCL2 A 258 

FPL L ( J )=F PL A 259 

F TL L ( J )=F TL A 260 

FCL L ( J )=FCL A 261 

FFLLt J ) =F FL A 262 

FSHFTL ( J) =F SHF T A 263 

BYCLL 1( J ) =6YCL 1 A 264 

BYCLL 2(J ) =BYCL2 A 265 

8PL L ( J )=B PL A 266 

BTLLt J )=3 TL A 267 

BCL L ( J ) =B CL A 268 

BSHFTL (J ) =BSHF T A 269 

FGLLK J)=FGL1 . A 270 
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n o o 


FGL L 2 ( J ) = FGL 2 A 271 

8R0TL ( J ) = BR 0 T A 27? 

FRUTL ( J ) =FRJ T A 273 

IF (KSAT.NE. 10) 30 TO 50 A 274 

IF (K.NE.O) 30 TO 51 A 275 

IF (J.LT.5) 30 TO 46 A 276 

GO TO 55 A 277 

50 I A=K $A T A 278 

IF ((G(J)-G(J-l)!.LE. 0.0*9! GO TO 52 A 279 

G( J ) = G ( J )-0. 10 A 280 

0=0-1 A 281 

K=K*1 A 282 

GO TO 46 A 283 

51 IF (J.E0.5) 30 TO 55 A 284 

0=0 + 1 A 285 

52 0A = J- 1 A 286 

53 '•1 =2 70 + J A 287 

DO 54 K= J , M , 5 A 288 

54 FGL L 1(KJ = 0. A 289 

0=0+1 A 290 

IF ( J .LT. 6) 30 TO 53 A 291 

IF (MM.GT.JA! A A 292 

55 IF (I.EQ.l) FIMM=FFLL (MM) /{ PT*FLOAT (NFC) ) A 293 

VV=3.*PI* EP*PF A 294 

A 295 

LOSSES AND EFFICIENCY UNDER LOAD A 296 

A 297 

M=0 A 298 

56 M = M + 1 A 299 

UA = G( M ) A 300 

F I< M I =FFL L < M ) / ( P T*F LOA T( NFC ) ) A 301 

CDOIM ) =F I <M)/AS A 302 

IF (MM.NE.i) GO TO 57 A 303 

T TB ( M ) =T2 2 A 304 

GO TO 58 A 305 

5 7 TTB ( M ) =( { ( T2- T 22) ^F I { M ) ** 2 ) + { T2 2 *F I M-M**2-T2*F I ( 1 )**2I ) / ( FI MM**2-FI A 306 

1( 1 )** 2 ) A 307 

53 RRB(M)=(i.OE-6)*RR*(l.O+ALPHAR*(TT8<M)-20.))*ZG A 308 

PRIM) = FI( M)*FI (M)*RR3 < M) A 309 

EHM!=FUM!*RKBCN! A 310 

?S(M)=(3.*(PI*UA)**2)*RRA{M) A 311 

W3L ('I ) =W3* (8CLL(MJ/3CLL<1) )**2 A 312 

ST<M)=WT* (BTLL(M)/BTLL(l) ) * *2 A 313 

WA( M) =VV*UA/1000. A 314 

AK V A( M I =W A ( M ) /PF A 315 

WMISIM )=AKVAI M)*10.0 A 3 16 

PP(M)=((GF*UA)**2.0+i.0)*WN A 317 

EX(M)-(EZ(M)-L.OI*PS(M)*(CL/HM) A 318 

SP I M)=PP( M)+PR( MI + PS( M) +EX( M) +ST ) + WF + WQL ( M) +WMIS { M) A 3 19 

E< v| )=( WA( M|/(WA (M)MSP(M) -WF) /1000. !!*130. A 320 

Ei(MI=(WA IM)/(WA(M)+SP(M) /1000. ) )*100. A 321 

IF <M.LT,JA) 30 TO 56 A 322 

WRITE (6,59) PF f (Gl I) 1 1 = 1 *5) A 323 

IF UiAT.'iv*. 1C) 3 (DA) =3 ( J A ) + FLJAT(K) *0.13 A 324 

IF U>.\T. lc.101 3(JA+l) =3(JA+i)+FLOAf(K)*0*13 A 325 

WRITE (o,60) (BILL (I) ,1=1,5) t (3CLL< I ) ,1=1,5), (BYCLLH I ) , 1= 1, b ) , ( A 326 


60 
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ICLL2< l ). 1 = 1, 5), (BPLLI I ) , I =1 . 51 . ( 8 SHF TL (I ) . I = 1 , 5) , ( BROT L( A 327 
WRITE (6,61) (FGLLll I) ,1=1 ,5) , IFGLLZll) ,1 = 1,5) A 328 

WRITE (6.60) (FTLL< I ).I = 1 .5) , (ECLLU I .1 = 1.5) .(EYCLLl (I ), 1 = 1. 5), ( FT A 329 
1C LL2( I ) • I = l , 5) , (FPU. 11). I *1.5) . (F SHF TL ( I ) , I = l , 5) , ( FROTL ( I) .1=1.5) A 330 
WRITE (6.62) ( FFLL ( I ) , I = 1 . 5) ,(FI (It ,1=1 .5) • ( COO ( l),( = l«5) ,(EF( t ) , I A 331 
1=1.5) A 332 

59 FORMAT ( 1H1, 26X47HALTERNAT0R LOAD CHARACTERISTICS (RATED VOLT AGE, F A 333 

15.2. 14H POWER FAC TQR ) /27X. 1 1 1 SH 1 . //7XI IHPERCENT KV A, l 2X , 2PFI A 334 

27.0, 4F 19. 0// 7X25HFLUX DENSITIES (KL/SQ-IN)) A 335 

60 FORMAI < l0X5FTEETHl5X5F19.2/10X4HCOKEl6X5F19.2/10X26HyOKE - CVLIND A 336 

1RICAL SECTION, F13.2.4F19.2/14X.14H END SECTI ON, 2X.5F19 . 2/ 1DX . 12H A 337 
2R0T0R - P0LE.8X.5F 1 9. 2 /18X. 11HEND SEC TI ON, l X .5F19.2/ 13X , 12HT APEREO A 338 
3 SEC..5F19.2) A 339 

61 FORMAT l IH /7X. 12HAMPERE— TURNS/10X, l 2HMAI N AIR GAP ,8X, 5F19 . 2/ 10X , 1 A 340 

12HAUX. AIR GAP , 8X, 5F 19. 2 ) A 341 

62 FORMAT (1H « 6X l 8HT0TAL AMPERE TURNS5 X5F 1 9. 2 //7 X20HF I ELD CJRRENT (A A 342 

IMPSI3X5F 19.2 /7X21HC UR RENT DENS. (F I E LD) 2 X5F 19 . 2/ 7X 1 1HF I ELD VDLTS12 A 343 
2X5F19.21 A 344 

WRITE (6.63) ITT8( 1) ,1=1, 5) .( TTA< I) ,1=1,5) . (RRBt I) ,1=1,5). IRRA(I), A 345 
1 I =1,5). (PR( I ). 1 = 1, 5) .WF.WF.UF, WF.WF, ( ST( I) .1 = 1 ,5) . (WOLI I ) , I = 1, 5 ) , ( A 346 
2PP( I ). 1 = 1,5) ,(PS(I). 1 = 1.5) . (EX1I t . 1 = 1.5) ,( WMIS(I) ,1 = 1,5) ,(SP( I ), 1= A 347 
31. 5), ( AKVAI I ). 1 = 1.5) ,( WA< I ) .1=1 ,5) , (E( I) ,1 = 1 ,5) A 348 


63 FORMAT ( 1HK, 6X20HTEMPERATURES (0E3.C) /10X5HFI ELD1 5X5F19. 2/ 10X8HARW A 349 
1ATURE12X5F19. 2/7X1 8HRESISTANCES OHMS) / 10X5HFI ELD1 5X5F 19. 2/ IDX8HAR A 350 
2NATURE12X5F19.4//7X25H ALTER NAT OR LOSSES ( WATTS)/ 10X5HFIELDI 5X5F19. A 351 
32/10X7 HW I NO AGE 1 3X 5F19. 2/1 OX 12H STATOR T00TH8 X5F19 .2/ 13X l US T ATOR C3 A 352 
4R E9X 5F1 9. 2/1 OX 9HP0LE FACE1 1X5F 19. 2/10X1 3HSTAT0R COPPER7X5F19 . 2/ 10X A 353 
54 HE DDT 16X5F19.2/10X10HMI SC. L3ADI 0X5F1 9.2/10X5HT0T AL15X5FI 9 .2// 7X2 A 354 
63 HAL TFRNA TOR DUTPUT ( K VA ) 5F 19. 2/7 X22H ALTERNATOR OUTPUT (KW)1X5F19. A 355 


72//7X. 21HEFFIC IENCY (PER CENT) /10X .1 6HELECTR0-MAGNET I C . 4X . 5F 19 . 2) A 356 

IF (WF.GT.O.) WRITE (6,64) ( El ( I ) , l = 1 .5 ) A 357 

64 FORMAT 1 10X, 8HQVER-ALL . 12X, 5F1 9.2) A 358 

IF I ITRTN.NE.10) WRITE (6,65) ITRTN A 359 

65 FORMAT ( 1H . 20X , 2 1 H ( S TAR T I NG WITH COLUMN , I 3 .52H ALL CALCS ASSJME J A 360 

INIFORM AIR-GAP FLUX DISTRIBUTION)) A 361 

KSAT=I A A 362 

IF (KSAT.EQ.10) G3 TO 67 A 363 

J=4*KSAT A 364 

M =J— 3 A 365 

WRITE (6.66J (CMPNT(K).K=M,J) A 366 

66 FORMAT (IH .20X.4A6.13H IS SATURATED) A 367 

67 IF (I.GE.NPF) GO TO 68 A 368 

GO TO 45 A 369 

C A 37 0 

C INITIALIZE VARIABLES USED IN NO-LOAD MAGNETIC CALCS A 371 

C A 372 

68 P F=PFC( 1) A 373 

DO 69 J =1 « 220 A 374 

69 FALL ( J i =0. A 375 

C A 376 

c A 377 

C CALCULATE NO-LOAD SATURATION DATA A 378 

C A 379 

C A 380 

F GXL =0. A 381 

I DELR=15 A 382 
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EDD=VM IN A 383 

ITRTN = 10 A 3 84 

J=1 A 385 

70 PPL=F0*EDD A 386 

CALL MAGNET A 387 

IF (ITRTN.EW.O) I TR TN = J A 388 

IF (<SAT.NE.10J 30 TO 71 A 389 

Ftt. L ( J )=P PL A 390 

QPERV { J)=100.*ED0 A 391 

QVLLI J )=EE*EOD A 392 

OVLNIJ )=.IVLL<J>/SQRT{3.) A 393 

F I( J >=FFL/(PT*FLC)AT(NFC) ) A 394 

FYCLL 1(J)=FYCL1 A 395 

FYCLL 21J)=FY0L2 A 396 

FPlLI J )=F PL A 397 

FTLL ( J ) =F TL A 398 

FCLLI J )=FCL A 399 

FFL L ( J )=F FL A 400 

FSHFTL(JJ=FSHFT A 401 

BYCLL l(J)=BYCLl A 402 

BYCLL 2 ( J ) =BY3L2 A 403 

BPLL ( J J =B PL A 404 

BTLLI J )=3 TL A 405 

BCLL( J )=BCL A 406 

BSHFTL(J)=BSHFT a 407 

FG.L 1< J)=FGL1 A 408 

FGLL21 J)=FGL2 A 409 

BROTL ( J ) =BRO T A 410 

FROTL I J)=FROT A 411 

IF (J.EQ. 10) GU TO 74 A 412 

J=J+1 A 413 

EDD=EDU+FL0AT(lDcLR)/100. A 414 

GO TO 70 A 415 

71 EDD=EDL)-F LOA T( 1DELR) /100. A 416 

IF (IDELR.GT.5) 30 TO 72 A 417 

IF (IDELR.GT.2) 30 TO 73 A 418 

GO TO 74 A 419 

72 IDELR=(I0ELR/6)*5 A 420 

EDD=EDO+FLOAT( IDELR) /100. A 421 

GO TO 70 A 422 

73 IDELR=(I0ELR/3)*2 A 423 

ED0=£DD+FLOAT ( IDtLRJ/100. A 424 

GO TO 70 A 425 

74 WRITE (6* 75) ( QPE RV I K. ) .K = l , l 0 ) , ( QVLN ( K ) ,K= 1 . 10 ) , ( QV LL ( K ) , < = 1, 1 0 ) , I A 426 

1FI(K),K=1,10),(FQLL(K) ,K=1,10) A 427 

75 FORMAT ( 1 HI , 5 OX23HNO-LOAD SATURATION DAT A/51 X23H A 428 

1 //2X 7HV0LTA3E/5X7HPERCENT6X10F11.27/5X12HLINE-NEJTRA. 1X10F11. A 4 29 

22/5X9HLINE-LINE4X10F11.2//2X13HFIELD CURRENT 3 X 13 FI 1 .2/ / 2< , 1 3HF. UX A 430 

3PER POLE, 3X, 10F 1 l . 2//2 X , 1 4NF LUX DENSITIES) A 431 

WRITE (6, 76) (BCLL(K) ,K=1,10) ,(BTLL(K) ,K=1 ,10) ,(BYCLL2(K),K=1, 10), A 432 

l( BYCLL UK ),K = 1, 10) , (3 PL UK) ,K=1 ,10) , (BSRFTLU) ,K=1 , ID) .(BROTLI < ), < A 433 

2=1. 10) A 434 

76 FORMAT ( 1H , 4X , 4HC0RE , 9X, 10F1 1 . 2 /5 X ,5 HTEETH , 3 X , 10F 1 1 . 2/ 5X , 4HY0 < E/ 7 A 435 

IX , 1 1HEND SECTION, 10F 1 1 . 2 /7X, 1 OHC YL. SECT. , 1 X , 10 FI 1 . 2/ 5X , 5H R0T3R/7R A 436 

2»4HP0LE,7X, 10F 11.2/ 7X, 11REND SEC TI ON , 10F1 1 . 2/ 7X , 12HT APEREO SEC., FI A 437 

30.2.9FU.2) A 438 
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WR ITE (6, 77) A 439 

77 FORMAT ( IHK, l X , 12HAMPE Rt- T URN S) A 440 

WRITE (6,78) (FGLLUK) ,K=1,10I ,(FGLL2(KI ,K=l,l3) A 441 

78 FORMAT (1H ,4X,12HMAIN AIR-GAP*F12.2,9F11.2/5X*11HAUX AIR-GAP, 2X, I A 442 

10F11.2) A 443 

WRITE ( t>, 76) (FCLLI K) ,K=l , 10) , (FTLL( K) ,K= l ,10 ) , (F7CLL2(K), K=l, 10), A 4AA 

1( FYCLL 1( < ) » K. = l , 10) , ( : PLL(K) ,K=1 ,10) , (FSHFT L( K ) , K= l , 10 ) , ( FR OT L I < ) , < A 445 

2= 1, 10 ) A 446 

WRITE (6,79) (FFLL(K) ,K = l ,10) A 447 

79 FORMAT ( IHK , 4X , 5HT0TA L , 8X , 1 OF 1 1 . 2 ) A A A 8 

IF ( ITKTN.NE. 10) WRITE (6,65) I TRTN A A A9 

IF (KSAT.EQ.10) GO TO 3 A A50 

J=4*KSAT A A5 1 

I =J- 3 A 452 

WRITE (6,66) (CMPNT(K) ,K = I , J ) A 453 

GO TO 3 A 454 

END A 455- 
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Input; output 

Branch point or decision 

Arithmetic statement 

Call statement for 
subroutine 


Numbers in upper left are external 
formula number in FORTRAN program 


No-load rated voltage 
magnetization characteristics 
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A 


Calculate short-circuit 
ratio (SCR) and short- 
circuit ampere-turns (FSC) 


/ Write name\ 
I of alternator | 
Icomponent thaty 
\ saturated J 


WRITE: 1 

FSC, SCR, 
materials, 
no-load 
rated-voltage 
magnetization 
character- 
istics 


(See sample 
output, 
p. 50) 


* 10 \ KSAT , 


No-load pole-face 
loss calculations 


No-load rated-voltage, 
tooth, and core-loss 
calculations 


Arrange load points G 
in increasing order 


Define MM (index for 
rated load point) 


Define NPF (number of 
power factors for which 
load characteristics will 
be calculated) 
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19 


0 


1 
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"No-load saturation I 
data" calculations J 
(see sample output, 
p. 54) 
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c 

c 

c 

c 


c 


c 
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SUBROUTINE SUBLUN 

COMMON /COM1/ A,3L,B2,BCQIL,BV,CC,CE,CF,CM,CQ»DCOIL»DF, DR* DWt DWl,E 
IC « E L • E W,F E »FK.l,FStS2,SAMMAR,GE,HO,HT ,HV,HW,HX,HY, I PN» IPX* I QQ , III?. 
2TS,P2»P3,P4,Pi3AfPC,PCOIL,PTl,PT2 ,PW1 ,PW2 , QN, RC , RD, RF, RGl, RT, S, SF, S 
3HLUSI ,SIGMA,SK ,SNI ,SPS, STATE T, TC , TF ,TST , T YE , T YR W A ,/ R, WC, W 1 , 4 ROT 3 
4R, WTJTAL, WYOKE , XF ,XQ, XR,XU,YY 

COMMON /CUM 2/ AC,AS,33,BK,B0,BS,C,C1,D1,DI,DISH,DISH1,DU,DYC,EE,E? 
i f Ft GCt HC,HM,NFC f PA 1 , PE t PF , P I , PL , PT , RPM , RR , RS , S C , SO , SH , SHL2, SN, $S, T 
2S f TTt TYfTYLfWFfWLtXA, XB.XDfXL 

COMMON /C 0M4/ ABfALPHARt ALPHAS,BG,FGML ,FH»FOtG,GA, PA, PFC,3Q,RY t T 1 , 
1TU,T2,T22, T3 ,VMIN,Z3 ,ZZ 

COMMON /COM 5/ A A , AS2 , C L , C P *0 Y2 , 0 Y3 , MS , LG2 , PI , P 10 , P5 , P6 , P7, P8,P9,PA 
L2*PPL,PPL l,PPL2,PPL3t PPL^fPXt TY2 ,TY3 ,TYPY 

INTEGER TYPY,ZZ 

REAL L TSt LG2 

DIMENSION DA ( 8 ) t DX<6), 0YC8), DZ<3), G(5>, PFC<5>, PA(ll) 

NAM EL 1ST /RATING/ VA,EE,EP,F ,RPM f I PX , PFC , G , VM I N/ S T ATOR/ DI , DU,CL,HY 
1,BV,SF,LTS, WL,BK/SLOTS/ZZ,BO,63 , BS . HO .HX » HY, HS , HT , I QU/ M I NONG/R F, SC 
2» EL.YYtCf DW, SN,SN1 ,0W1 ,CE ,SD , PBA ,SK ,T1 , RS t ALPHAS t T 1 It T S T/ A IRGAP/ GC 
3, G2» WF tLG2/R0 TOR/C I ,C P tC M ,C Q , PW I ,PW2 ,PTL , PT2 , PL,WROTOR, DI, PAL* PA2, 
4PA,PE,GAMMAR,DI SH,31 SHWSHL1 , SH L2/ YORE / TY , TYL , DYC , W YOK E, T/ E, TY R, T/ 
52, T Y3/FIELD/PC0IL, OC 01 L,PT,RD ,RT ,T2 , BCOI L , TF t T 2 2 , RR , ALPHAR,TYP//PE 
6RMCE/P It? 2,P3,P4,P5,P6,P7,P10 

DATA DA,DX,3Y,DZ/0.05,0. 0 72, 0.125 , 0.165 , 0.22:>, 0.438,0 .688, 1.5, 0.00 
10124, 0.00 021, 0.00 021 , 0. 00084 , 2*0 . 0 01 8 9 , 2 *0 . 0 0 0 1 24 , 2 *0 . 000 3 4, 0.3013 
29, 0.00335, 0.007 54, 0.0302,3*0. 000 12 4 ,2*0.30335 ,0.00754,0.0134,0.033 
32/ 

WRITE (6,1) 

FORMAT ( 1H L , I 4X , 5 OH* * * STATIONARY UUTSIDE-COIL LJNDELL ALTERNATOR 
i*** ) 

DO 2 1=1,5 
PFC ( I ) =0. 

C 1 = 0 

R S = 0 • 694 
RR=0. 694 
Al ; >HAS=0. 00393 
ALP HAR =0. 00393 
P 1 = 0. 

P 2=0 . 

P 3=0. 

? 4= 0 • 

P 5 = 0. 

P 6 = 0. 

P 7= 0. 


B l 
B 2 
B 3 
B 4 
B 5 
B 6 
B 7 
B 8 
B 9 
B 10 
B 11 
B 12 
B 13 
B 14 
8 15 

B 16 
B 17 
B 18 
B 19 
B 20 
B 21 
B 22 
8 23 

8 24 

B 25 
B 26 
B 27 
B 28 
B 29 
B 30 
B 31 
8 32 

B 33 
B 34 
B 35 
B 36 
8 37 

B 38 
B 39 
B 40 
8 4 L 

B 42 
8 43 

B 44 
S 45 
B 46 
8 47 

8 48 

8 49 

B 50 
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P 8 = 0 • 

B 

5 l 

P 10=0. 

B 

52 

SPS=0. 

B 

53 

R T = 0 • 

B 

54 

TF=25. 

B 

55 

T ST= 2 5 . 

B 

56 

PE = 0. 

B 

57 

PL = 0 . 

B 

58 

SHL 1=0. 

B 

59 

P BA =60 • 

B 

60 

SN= i . 0 

B 

61 

0Wi=0 

B 

62 

c w=o 

B 

63 

CP = 0 

B 

64 

EL = 0 

B 

65 

CM = 0 

B 

66 

G( 1 )= 0 . 

B 

67 

G ( 2 ) = 0 . 75 

B 

68 

G { 3 ) = 1 .00 

B 

69 

G( 4 )= 1 . 25 

B 

70 

G< 5) = 1.50 

B 

71 

VM I N= 0 . 7 

B 

72 

co=o 

B 

73 

W F= 0 

B 

74 

T Y = 0 

B 

75 

XF =0 . 

B 

76 

EP = 0. 

B 

77 

EE = 0. 

8 

78 

IPM = 3 

B 

79 

P N = 3 . 

B 

80 

IPX=0 

R 

81 

F =0 • 

B 

32 

RPM =0 . 

B 

83 

SF = 0. 

B 

84 

L T S=0 . 

B 

85 

WR3T3R=0. 

B 

86 

HV = 0 • 

R 

87 

BV = 0. 

B 

88 

BCD IL = 0 . 

B 

89 

SK = 0 

B 

90 

P A 1 = 0 • 

B 

91 

P A2=0 . 

B 

92 

GAM MAR =0.2 83 

B 

93 

DY2=0 . 

B 

94 

T Y 2=0 . 

B 

95 

NFC = 2 

B 

96 

01=7.0 

B 

97 

WY3KE =0. 

B 

98 

03 3 1 = 1, 11 

B 

99 

PA( I ) =0. 

B 

100 

READ ( 5, RATING) 

B 

10 l 

READ ( 5, STATOR) 

B 

102 

READ [ 5, SLOTS) 

B 

103 

READ (5,WINL)NG) 

B 

104 

READ ( 5,A IRGAP ) 

B 

105 

READ ( 5, ROTOR) 

B 

106 
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R EAD ( 5. YOKE ) B 107 

READ I 5* F IELD) B 108 

READ (5.PERMCE) B 109 

IF (TYPY.E0.6) NFC = l B 110 

IF (NFC.EO.l) WRITE (6,4) B 111 

IF (NFC.EQ.2I WRITE (6,5) B 112 

FORMAT (IH , 31 X , 16H (LINE FIELD COIL)) B 113 

FORMAT (1H • 31 X • 17H ( TWO FIELD COILS)) B 114 

P F=PFC( 1) 8 115 

If (EP.EO.O.) EP=EE/1. 732051 B 116 

IF (EE.EO.O.) EE=EP*l. 732051 B 117 

IE (DWi.NE.O.) SH=DW 1 B 118 

If ( IPX.EQ.O.AND.RPM.NE.O.) I P X= ( E * 12 0. ) /RPM B 119 

PX=IP.X B 120 

IF (RPM.EO.O..AND.PX.NE.O.) RP M= ( F *12 0. ) /PX B 121 

IF (F.EO.O.) F=PX*RPM/120. B 122 

HW=HY- HD- HT B 123 

00=100 B 124 

IF IZ Z.NE.3) GO TO 6 B 125 

B1 = (H3FHT-HS)*( 6.283185/0O+B3 B 126 

B?=B1M 6.283185*HW/00) B 127 

BS=(82«-B3)/2. B 128 

6 CONTINUE B 129 

P I=( VA*1000. )/(EE*SQRT(3. ) ) B 130 

IF (Z7. EO.l. OR. ZZ.EO. 5) BO=BS B 131 

I ZZ=Z7. B 132 

D B=. 25 B 133 

IF (DU.GE.8. ) 0B=0. 5 B 134 

FE=3.l416*(PC0lL+DC0IL)/2. B 135 

OR=D 1-2.* GC B 136 

IF (SHL1.E0.0.) SHL1-CL B 137 

PHW=(PWl+PW2)/2. B 138 

IF (PE.EO.O.) PE=(PX/3. 1415927)*(ARSI NIPHW/DR) ) B 139 

HC=( Dli-D I - 2. C*HS)*0.5 B 140 

IF (PA1.E0.0.) PAI=P A I 1 ) B 141 

IF (PA2.EQ.0.) PA2*P A ( 11) B 142 

IF (PA( I) .LE.O. ) PA I l ) =P A1 B 143 

IF (PA( 11 ).LE.O. ) PA (11) =PA2 B 144 

DY?=DISH1+2.*G2+2.*TYE B 145 

DY3=DYC-2.*TY B 146 

IF (TY2.E0.0.) TY2=TY3 B 147 

Z Y=0.7*HS B 148 

DO 7 1 = 1, 5 B 149 

7 IF ( G( I )• GT. 9. ) G(I )=G( I ) /100. B 150 

QN=QO/( PX*PN ) B 151 

CS=YY/(PN*ON) B 152 

C B 153 

C CHECK FOR ERROR CONDITIONS B 154 

C B 155 

IF (CS.GT.i.O. OR.CS.LT. 0.5) WRITE (6,9) CS B 156 

IF (EP*EE.EO.O.. OR. ABS(EE/EP-1. 732051). GT.O. 01) WRITE (5,13) B 157 

IF (PX*F*RPM.E0.O..0R.ABS(F-PX*RPM/l2O.) .GT.0.1) WRITE (5,11) B 158 

IF (HC.LT.ZYI WRITE (6,12) HC.HS B 159 

IF (DC0IL.GT.DY3) WRITE (6,13) B 160 

IF 1RT.LT. 1.0E-10) GO TO 8 B 161 

IF <<((DCO IL-PCOIL )*BCOI L ) / (RT*RD) ).LE.2.*PT) WRITE (6,14) B 162 
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8 

9 

10 
11 
12 

13 

14 

C 

C 

C 

15 


16 

C 

C 

17 


18 

19 
C 
C 

c 

20 


21 

C 

C 

C 

22 

23 

24 

C 

C 

C 

25 


26 


G3 TO 15 

IF ( ( DCOIL-PCOIL)*3COIL/RD**2.LE.1.7146*PT) WRITE (6,14) 

FDRMAT ( 5 X , 2 7H CS (PER UNIT POLE PITCH) = *F7 • 3/ 13 < » 31 H CS MUST BE 
1BETWEEN 0-5 AND 1-0) 

FORMAT (1H , 3 8H EITHER PHASE OR LINE VOLTAGE IS WRONG) 

FORMAT { 1H , 44H FREQUENCY, RPM, OR NO- OF POLES IS IN ERROR) 

FORMAT ( 1H /5X54HDEPTH BELOW SLOT IS LESS THAN 70 PERCENT OF SLOT 
1DEPTH/ 10 X , 4HQ8 S=F 8.4/10X,4H SD=F8.4 > 

FORMAT ( IH , 34H FIELO SOIL 0-0- EXCEEDS YOKE I-D.) 

FORMAT (1H » 8 1H FIELO SOIL DIMENSIONS ARE TOO SMALL FOR THE SPECIF 
1 1 ED NO. OF TURNS AND WIRE SIZE) 

DETERMINE STATOR STACKING FACTOR 

IF ( SF -NE -0. ) GO TO 1 7 
IF (LTS.EQ.O. ) GO TO 16 
SF = 1 • 0- ( 12.5E-4/LTS) 

GO TO 17 
SF=1.0 


SS=SF* (CL-HV* B V ) 

S IGMA= (54.£3/0I**2)*(PF/SS)*( VA/RPM) 

VR=0.262*DR*RPM 

TP = 3. 142*DI/PX 

TS = 3 • 142*DI/3Q 

IF IZZ-4) 18,19,18 

TT = ( • 6 67* H S+D I )*3.142/QQ 

GO TO 20 

TT={DI+2. 0*H0+1 -333*BS)*3. 1416/QQ 

CALCULATE CARTER COEFFICIENT 

IF (ZZ.GT - 1-AND.Z.Z.LT. 5) GO TO 21 
CC=(5.0*3C+BS)*TS/( (5.0*3C+BS)*TS-8S*BSI 
GO TO 22 

3C = ( 4 • 44* GC + 0. 75*B0)*T5 
CC=QC /(QC-BD*BQ| 

PITCH FACTOR AND SKEW FACTOR CALCULATIONS 

CF=SIN(YY* 1. 571/1 PN*QN> ) 

IF (SK) 23,23,24 
FS= 1-0 
GO TO 25 

FS=( SK/TP ) * 1 • 5707 
F S= ( i./FS)*(SIN(FS )) 

CHECK IF WINDING HAS INTEGRAL NO. OF SLOTS PER POLE PER PHASE 


0 = 1.0 

IF 1P8A.GT.61. 0) D=2. 0 
IZY = IPX*IPN 
I DM=0 

IDM= I OM+I Z Y 
IF IIQQ-IDM) 28,27,26 


B 163 
B 164 
B 165 
B 166 
B 167 
B 168 
B 169 
B 170 
B 171 
B 172 
B 173 
B 174 
B 175 
B 176 
B 177 
B 178 
B 179 
B 180 
B 181 
B 182 
B 183 
B 184 
B 185 
B 186 
B 187 
B 188 
B 189 
B 190 
B 191 
B 192 
B 193 
B 194 
B 195 
8 196 
B 197 
B 198 
B 199 
B 200 
B 201 
8 202 
B 203 
B 204 
B 205 
B 206 
B 207 
B 208 
B 209 
8 210 
B 211 
B 212 
B 213 
B 214 
B 215 
B 216 
B 217 
B 218 
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C B 219 

C CALCULATE DISTRIBUTION FACTOR FOR INTEGRAL SLOT WINDING B 220 

C B 221 

27 DF = SIN{ 1.57L*D/PN) / ( QN*D* SI N < 1 . 57 L / < PN*QN))> B 222 

GO TO 32 B 223 

C B 224 

C CALCULATE DISTRIBUTION FACTOR FOR FRACTIONAL SLOT WINDING B 225 

C B 226 

28 II3U=IQQ B 227 

I _2 B 2 2 8 

29 IF U IZY/I)*I.EQ.IZY.AND.UIQQ/n*I.EQ.IIQQ> GO TO 33 B 229 

IF I l .GT. IZY) GO TO 31 B 230 

1= I + L B 231 

GO TO 29 B 232 

30 IZY=IZY/I B 233 

I I3Q= I IW/I B 234 

GO TO 29 B 235 

31 FN3 = 1 I QQ B 236 

DF = S IN ( 1. 57L*D/PN ) /( FNQ*D* SI N( I. 571/ (F NQ*PN) I ) B 237 

32 EC=33*SCt t CF' c FS/G B 238 

C B 239 

C COMPUTE ARMATURE CONDUCTOR AREA B 240 

C B 241 

IF (DW1) 33,33,34 B 242 

33 AC=0.785*DW«OW*SNl B 243 

GO TO 46 B 244 

34 Z Y = 0 . 0 B 245 

DT = AM INK 0W,3W1) B 246 

DG=AMAX1< DW.DWl ) B 247 

35 IF (0T-.05) 33,38,36 B 248 

36 JA=0 B 249 

37 J A = J A + 1 B 250 

IF IDT-DAIJAJ) 39,39,37 B 251 

38 0=0 8 252 

IF IZY) 45,45,58 B 253 

39 IF IDG-O. 188) 40,40,41 B 254 

40 C Y =DX I JA- I) B 255 

CZ=DX I JA ) B 256 

GO TO 44 B 257 

41 IF (DG-U.75) 42,42,43 B 258 

42 CY=DYI JA- 1 ) B 259 

CZ=OY ( JA ) B 260 

GO TO 44 B 261 

43 CY=DZ l JA- 1 ) B 262 

C Z = DZ ( JA ) B 263 

44 D=CY+(CZ-C Y)* (DT-DAI JA-1 ) ) / (DAI JA) -OA ( JA-I ) ) B264 

IF IZY) 45,45,58 B 265 

45 AC=(DT*UG-D)*SNl B 266 

C B 267 

C CALCULATE ENO EXTENSION LENGTH B 268 

C B 269 

46 IF (EL ) 47,47, 55 B 270 

47 IF (RF) 48,48,54 B 271 

48 IF IPX-2. 0) 49,49,50 B 272 

49 U= L . 3 B 273 

GO TO 53 B 274 
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50 

IF (PX-4.0) 51,51,52 

B 

275 

51 

0=1.5 

B 

276 


03 TO 53 

B 

277 

52 

U = 1.7 

B 

278 

53 

EL = 3. 1416* U* Y Y* (OI + HS) /QQ+O. 5 

B 

279 


GD TO 55 

B 

280 

54 

EL = 2.0*CE+< 3. 1416* (0. 5*HX+DB) >+< YY* TS*TS/ ( SORT ( TS* T S-BS *B$ ) ) ) 

B 

281 

55 

HM=CL+EL 

B 

282 

C 


B 

283 

C 

CALCULATE ST A TDK RESISTANCE 

B 

284 

c 


B 

285 


A =P I*SC*CF/(C*TS) 

B 

2 86 


RY=SC*QQ*HM/( PN*AC*C*CJ 

B 

287 


RGi=( l.E-6)*RS*< 1.0+ ALPHA S*(T ST-20. ) ) *RY 

B 

288 


S=PI/(C*AC> 

B 

289 

c 


B 

290 

c 

COMPUTE FIELO CONDUCTOR AREA 

B 

291 

c 


B 

292 


IF <RT) 56,56,57 

B 

293 

56 

AS=.7854*RD*RD 

B 

294 


GO TO 59 

B 

295 

57 

Z Y= 1 • 0 

B 

296 


DT = AM INK RT, RD ) 

B 

297 


DG= AM A Xl( RT, RO ) 

B 

298 


GO TO 35 

B 

299 

58 

AS=DT*DG-D 

B 

300 

C 


B 

301 

c 

COMPUTE FIELD RESISTANCE 

B 

302 

c 


B 

303 

59 

ZG=(PT*FE/AS)*FLDAT(NFC > 

B 

304 


F< 1 = ( 1 .E- 6)*RR* ( 1 . 0+ AL PH AR* ( TF-20. ) ) *ZG 

B 

305 

C 


B 

306 

c 

NO LOAD MAGNETIC CALCULATIONS 

B 

307 

c 


B 

308 


GA= 3. 1416*0 I* (CL-HV*3 V) 

B 

309 


AG2 = 3 • 141 6* ( DI SH 1 + G 2 ) *LG2 

B 

310 


GE = CC * GC 

B 

311 


IF (C 1) 61,60,61 

B 

312 

60 

Cl^I .649* A LOO IPE 1+1.359) 

B 

313 

61 

CW = ( 0.707 / 1.7 32 )*C 1*0 F 

B 

314 


TG=EE/(CW*£C*RPM)*6.0E6 

B 

315 


BG= TG/GA 

B 

316 


FH=BG*GE/ 0.0031 9 

B 

317 


IF (CP) 62,62,63 

B 

318 

62 

CP=PE* (ALQGOC / TP) *.0378+1. 191) 

B 

319 

63 

FO=TG*CP/PX 

B 

320 

C 


B 

321 

C 

DETERMINE DEMAGNETIZING AMPERE TURNS (RATED LOAD) 

B 

322 

c 


B 

323 


IF (CM) 64,64,65 

B 

324 

64 

A A = S I N ( 3. 142* PE ) 

B 

325 


AB = S IN C 1. 571*PE )*4. 0 

B 

326 


CM = ( 3. 142*PE + AA) /AB 

B 

327 

65 

CONTINUE 

B 

328 


FGML= •45*EC*PI*CM*l)F/PX 

B 

329 

C 


B 

330 


SUBLUN 


75 



c 

REACTANCE FACTOR CALCULATION 


B 

331 

c 




8 

332 


XR=.07C7*A*DF/IC1*63) 



B 

333 

c 




B 

334 

c 

SPECIFIC ARMATURE SLOT AND 

END-TURN 

LEAKAGE PERMEANCES 

B 

335 

c 




B 

336 


FAC TOR = YY /( PN*UN) 



B 

337 


IF (P8A.LT.6i.) SO T3 66 



B 

338 


FF=.05*(24.*FACTOR-i. ) 



B 

339 


IF (FACTOR. GE. 0.667) FF=.75 


B 

340 


IF UZ. E^.5) FF =1 . 



B 

341 


G3 TO 67 



B 

342 

66 

PF=.25*(6.*FACTGK-1.) 



B 

343 


IF (FACTOR. GE. 0.667) FF=. 25^(3. *F AC TOR+1. ) 

B 

344 


IF {ZZ.E3.3) FF =1 ♦ 



B 

345 

67 

CX-FF / (CF*CF*3F*3P ) 



B 

346 


Z=CX* 20.0/(P\l*QN) 



B 

347 


BT = 3 . 1 42* DI/Q3-B3 



B 

34B 


ZA=BT*BT/( 16. 0*TS*GC) 



B 

349 


ZB = 0. 35*8 T/TS 



B 

350 


ZC=H0/3U 



B 

351 


ZD^HX* . 33 3/3 S 



B 

352 


ZE=HY/BS 



B 

353 


IF ( Z Z-2 ) 68,69,70 



B 

354 

68 

PC = Z*< ZE+ZD+ZA+ZB) 



B 

355 


GJ TO 74 



B 

356 

69 

?C = Z* ( ZC + ( 2. 0*HT/(SO+BS) ) ■MHW/BS1 + ZD+ZA+ZB) 

B 

357 


GO TO 74 



B 

358 

70 

IF (ZZ-4) 71,72,73 



B 

359 

71 

PC=Z*( ZCM 2.0*HT/(3U+tH) ) *■ (2. 0*HW/ (B1 + B2) ) + ( HXZ ( 3.*B2) )+ZA+ZB) 

8 

360 


GO TO 74 



B 

361 

72 

PC = Z*( ZC+0.62) 



B 

362 


GO TO 74 



B 

363 

73 

PC=Z*I ZE + ZD+I 0.5*GC/TS)M 0. 

25*TS/3C)+0.6) 

B 

364 

74 

E<=£l /{ 10. 0** ( 0.103*YY* TS+O 

.402) ) 


B 

365 


IF (DI-8.0) 75,75,76 



B 

366 

75 

EK = SOR T( EK ) 



8 

367 

76 

Z F= .61 2*A LOG ( 10.0*CS) 



B 

368 


EW = 6.2 8*EK*ZF*(TP*M0. 62- ( . 

2 2 8* ALOG 

(ZF) ) ) )/(CL*DF*DF) 

B 

369 


IF (SPS.NE.O. ) GO TO 78 



B 

370 


AA = TP-PW2 



8 

371 


AB = T S- 83 



B 

372 


IF (AA.LT.AB) GO TO 77 



B 

373 


SPS=0. 



B 

374 


GO TO 78 



B 

375 

77 

SP S = ( 6.67*CX/QN) /(2.*CL*3C) 



B 

376 


A A = ( ( A8- ^ A )*0. 5 ) * *2 



B 

377 


AA={ A A*CL )/(PW2-PWL) 



B 

378 


SPS=SP A 



B 

379 

C 




B 

380 

c 

LEAKAGE PERMEANCES USED IN 

MAGNETIC 

C ALCDLAT I ONS 

B 

381 

c 




B 

382 

78 

IF (P i.NE.O.) GO TO 79 



B 

383 


P 1=< 3 . 19* P A 1 1 / ( PL-C L ) 



B 

384 

79 

IP (P10.NE.0. ) GO TO 84 



B 

385 


IF (P2.NE.0. ) GO TO 80 



B 

386 

76 


SUBLUN 





80 


81 


82 

83 

84 


85 

86 

87 

88 

89 

90 


91 

92 

C 

C 

c 

93 

94 


C 

C 

c 


95 


96 

C 

C 


P 2=< 3. 19*PL*(PTH-PT2)*0. 5) /(TP-PHW) 


B 

387 

IF (P3.NE.0.) GO TO 81 


B 

388 

P 3=1 I 3. 19*PHM*0.64)/I (4./CL)*IPL-CLM-1.0>)*2. 


B 

389 

IF IP4.NE.0.) GO TO 83 


B 

390 

IF I IPX .L T.4) GO TO 82 


B 

391 

P4=l 0.64*3. 19*CL)/U 1 TP-PHW) /I PH W/2. ) 1 *1 . 0 ) 


B 

392 

P4=P4*( 1.0/( 1.0-2.0/FLQATI IPX) > I 


B 

393 

GO TO 83 


B 

394 

P 4=1 3. l9*CL*PHW/2.)/IDR-PTl-PT2> 


B 

395 

P 10=2.* ( P 3+P 2*P 4) 


B 

396 

IF (P5.NE.0.I GO TO 91 


B 

397 

A A=D ISH1/ 2. +G2+TYE 


B 

398 

A8=DYC/2.-TY 


B 

399 

IF (TYPY-1) 85,85,86 


B 

400 

P 5=10.53* (AA *0.425* C LG2-TY2) 1 


B 

401 

GO TO 91 


B 

402 

IF (TYPY-5) 87,88,89 


B 

403 

A A=( AB-AA+I TYL-CL ) /2. ) /2. 


B 

404 

GO TO 90 


B 

405 

A A= I DU- DI+TYL-CLI/4. 


B 

406 

GO TO 90 


B 

407 

A A=I TYL-CL )/3. 


B 

408 

A B= AB-AA 


B 

409 

P5=10.53*(AB*0.425*AA) 


B 

410 

IF ( TYPY.E0.6) P 5=0. 5*P5 


B 

411 

IF (P6.NE.0. J GO TO 92 


B 

412 

P 6=3. 1 9*1 I 0. 7854*1 DR-2.*PT2)**2) /I2.*PL-CU ) 


B 

413 

P 8=1 3. 1 9* ( GA 1 1 / GE 


B 

414 

P 9=1 3. I9*AG2 )/G2 


B 

415 



B 

416 

STATOR MINDING LEAKAGE AND ARMATURE REACTION 

REACTANCES 

B 

417 



B 

418 

IF (CO) 93,93,94 


B 

419 

A B=3. I416*PE 


B 

420 

C0=( 4.*PE+l.)/5.-SIN(AB)/3. 1416 


B 

421 

XL=XR*( PC+EW+SPS) 


8 

422 

F G 1 = BG * GE * 3 1 3. 


B 

423 

FG2 = ( ( FO*PX ) /P 91 *500. 


B 

424 

X D=( 0.45* EC*PI*CM*DF* 100. ) / 1 < PG1+FG2 > *PXI 


B 

425 

X0=< IC0)/ICM*C1))*XD 


B 

426 



B 

427 

FIELD LEAKAGE REACTANCE. SELF INDUCTANCE ANO 

TIME CONS T ANT 

B 

428 



9 

429 

S TATET=00*SC*DF*CF/<2.*PN*CI 


B 

430 

AA=P7*P1Q*PX 


B 

431 

A B=? .*( (P9*AA+P9*P5*P5*AA) /(P9«-AA> ) 


B 

432 

IF ( TYPY.E0.6) AB = P5+0.5*( I P9* AA) / I P9* AA ) > 


B 

433 

A B= AB* ( P T**2 )* 1 . OE-6 


B 

434 

IF (XP.GT.O.l GO TO 96 


B 

435 

XF=9.4247 8*F*AB*USTATET/PT)**2)*IPI/EP) 


B 

436 

A A=AA+CP*P 8 


B 

437 

GO TO 95 


B 

438 

S I = AB*l.0E— 2 


B 

439 

T C=SI/FK1 


B 

440 



B 

441 

SYNCHRONOUS ANO TRANSIENT REACTANCES CALCULATIONS 

B 

442 


SUBLUN 


77 



O NO wD vO o o o 


c B 443 

XA = XL +X0 B 444 

XB = XL+X3 B 445 

XU = XL ■*■( XF* XD ) / ( XF *XD ) B 446 

B 447 

WEIGHT CALCULATIONS B 448 

B 449 

IF <Z Z-3) 97,98,97 B 450 

7 WI=C ( OU+3 I )* (DU-31 )*3.1416> /4. B451 

IF IZZ.NE.4) WI =WI-U3*<BS*HS-(<HO+0.5*HT)*<BS-Ban ) B452 

IF (ZZ.EQ.4) WI=WI-Q3*(BS*BS*3. 141 6 /4 . +H0*B0) B 453 

GO TO 99 B 454 

8 WI=(DU-HC )*3.1416*HC 8455 

W I = W I *HS* ( I0I+2.*HS)*3. 1416-QQ*B3). B 456 

WI=WI+Q4* ( (H0 + 0.5*-IT) *<BS-B0) ) 8 457 

9 WI=WI*0.283*SS 8 458 

B 459 

RC = ( . 321*PT*FE*AS)*FLOATI NFC ) 8 460 

C 8 461 

WC=.321*SC*QQ*AC*HM B 462 

C B 463 

IF I W YOKE .NE » 0. I GO TO 100 8 464 

WY3KE=DYC-2.*TY 8 465 

W YOKE = ( l , / 12 . )*< TY2-TY3) * l D Y2 **2+0 Y2* WYOKE+W Y 0KE**2 ) +0 . 25* < TY3*WY3 B 466 

1<E**2) 8 467 

WY0KE=WY0KE-0.25*IDY2**2)*TY2 8 468 

WY3<E = WY3KE+(2.M31SHH-2.*G2+-TYE)*TYE*LG2) 8 469 

WY3<E=WY3KE-M IDYC-TY) *TY* (TYL+2.*TY3) ) 8 470 

WY0KE=WY3KE*3. 1416*0.283 8 471 

C B 472 

100 IF (WR0T3R.NE.0.) GO TO 101 8 473 

WR3T0R=(DISH1**2-DI SH**2)*< TYR-SHL2) B 474 

WR3 TOR =WRD TO R + I 0. 33 3* (31 SH 1 **2 + D I SH 1 *DR+ DR**2 I -DIS H**2 ) *( S HL 2-S HL l 8 475 

1) 8 4 76 

WR3TOR=WRUT3R+(UR**2-DISH**2)*(SHCn*{GAMMAR/0.283) B 477 

WR3 TOR =WR0T3R*0. 7854*0. 283 8 478 

WR0TUR=WR0T0R-<PX/2.)*{PA1+PA2>*<SHL1)*(0.283-GAMMAR) 8 479 

101 WTOTAL=W:+WH-R: + WYOKE«-WROTOR B 480 

RETURN 8 481 

END B 482- 


78 SUBLUN 


J 



Start 


SUBLUN 



Input/Output 



Branch point or decision 
Arithmetic statement 
Call statement for subroutine 


Numbers at upper left are external formula 
numbers in FORTRAN program 



SUBLUN 


79 






80 


SUBLUN 











Detail A 



SUBLUN 


81 












O U f\J fn >fr 


c 


SUBROUTINE OUTPUT 


C 

C 

C 

C 


5 

6 


COMMON /COMi/ A,B1,92,BC0IL,BV,CC,CE ,CF , C M, CQ, DCOI L , OF , OR, DW,DWl, E 
IC, EL, EW,FE ,FKlfFS»S2?SAMMARfGEfHO,3T *HV»HW,HX ,HY , I PN, IPX , l QQ, 177,. 
2TS,P2,P3,P4,PBA,PC,PC0I L,PTi , PT2 ,PWI , PW2 , QN, RC , RD, RF, RG l , RT , S, SF,5 
3HL i,S I , SIGMA, SK,SNl , SPS, STATE T, TC , TF , TST , T YE , T YR ,V A, VR, WC, W I , W ROTO 
4R, WTOTAL, WYQKE , XF ,XQ,XR, XU,YY 

COMMON /:0M2/ AC , A S ,3 3 , BK , BO , BS ,C ,CL , 01 ,D I , 0 I S H , D I S HI , OJ , DY C, E E, E* 
1,F,GC ,HC,HM,NFC,PAL ,PE . P.= ,PI ,PL, PT,RPM ,RR ,RS , S C , SO , S H , S HL2 , SN, S S , T 
2S, T T, T Y, T YL , WF , WL , XA , XB , XO , XL 

COMMON /COM 5/ AA,AG2,CL,CP ,0 Y2 , D Y3 , HS ,LG2 , Pi , PLO , P5 , PS , P7, P8,P9,PA 
12, PPL ,PPL 1,PPL2,PPL3,PPL4,PX,TY2 ,TY3 , T YPY 

DIMENSION STAR ( 5) , DASHC5) 

INTEGER TYPY 

REAL L TS, LG2 

DATA STAR ( 1) / 30H************* *** **************/ f DASH( 1 ) /30H 

l / 


WF1— 10. 

WRITE 16,1) VA,EE ,EP,PI ,PF,IPN,F ,1 PX * RPM 

FORMAT ( 1HL , 1 8H ALTERNATOR R A TI NG/ / 10 X , 1 5 H ALTERNATOR KV A, F 16. I / 10 
IX , I8H LINE-LINE VOLTAGE ,F12.0/10X,19H L l NE-NEUT • »/ OLT AGE, F 1 I . 0 / l 0< 
2, 14H PHASE CURREN T,F18.2/10X, 13H POWER F ACTOR , FI9 . 2/ 10X , 7H PHASES, 
31 22/ L OX, I OH FREQUENCY, F20.0/iOX,6H POLES , l 23/ L 0 X , 4H RPM,F27.l) 

IF (12Z-2) 3,5,2 

IF (IZZ-4) 7, 9, L i 

WRITE (6,4) 3S,HX,3Y,3S,IQQ,TS,TT 

FORMAT ( IHL , 1 3H STATOR SLOTS//5X10H T YPE-OPEN/ 54X , 9H *, 12XS 

1H* /6 2X1H* , 12X13*/55X2HHY,5X1H*,12XLH*/IDX3H BS ,F26 .3, IX 6H INCH 

2ES, 16X,1H*,12XIH*/10X3H HX,F26.3 ,L3X,9H *, 2 X8 H***** ***, 2X LH 

3*/ 1 OX 3 F HY,F26.3,23 XU* ,2 XI H* ,6X1H*,2XIH*/10X3H HS , F26 . 3, 2 3X 1H*, 2< 
41H* ,6X 1H* ,2X1H* /62X, L3* , 2 X8H ******** ,2 X 1H*2 X2HHS / 5 5X2HHX , 5X, 1H*, 12 
5X LH* / IOX1 3H NO. OF SLOTSI l6,23X,lH*,2X3H********,2XlH*/62< l H*, 2X1H 
6*,6X1H*,2X1H* / 10X1 13 SLOT PI TCH ,F18.3 ,1X6HINCHES , 16X1H*,2< 1H*, 6X1H 


7*, 2XIH*/54X9H * , 2 X83 ** ** **** ,2 X l H*/ ID X 1 1 H SLOT PITCH, 4i<IH* 

8, 12XIH*/10X15H AT 1/3 0 1 ST. ,F 1 4. 3 , 1 X6H I NCHES , 16 X 1 9H***** ******** 


9* /62X1H It 12XLH 1/62X1 4H1 BS 1/62X1H1 , 12X1 Hi) 

GO TO 13 

WRITE (6,6) BG,BS,30,HX,3T,HW,HS ,IQQ,TS,TT 

FORMAT ( IHL , 1 3H STATOR SLOTS//5X22H T YPE-P AR T I ALLY CL OS ED/ 67X4H- BD 

1-/5 7X1 OH * ,4X103* / 5 8 X2 HHO ,b X 1 H* , X IH*/ 57< 10H 

2 *, 4X1H*/10X3H BO , F 2 6. 3 , 1 X6H I NC3E S , 19 XI H*, 6X 1 H*/ 13 X 3H BS, F2 6.3, 

3L9X2HHT,4X1H* , 8X1H*/10X33 HO , F2 6 . 3 , 24 XI H* , 1 0 X 1 H*/ 1 0 X3 H HX , F26.3, 18 

4X6H * , L2X1H*/10X33 H T , F 2 6 . 3 , 23 X 1 H* , 1 2 X IH*/ 10X 3H HN , F2S . 3 , l 9X 2H 

5HW , 2X IH*, 12X IH* /I 0X33 HS,F26.3, 18X5H *2X3H********,2XIH*,2X2HH 


C 1 
C 2 
C 3 
C 4 
C 5 
C 6 
C 7 
C 8 
C 9 
C 10 
C 11 
C 12 
C 13 
C 14 
C 15 
C 16 
C 17 
C 18 
C 19 
C 20 
C 21 
C 22 
C 23 
C 24 
C 25 
C 26 
C 27 
C 28 
C 29 
C 30 
C 31 
C 32 
C 33 
C 34 
C 35 
C 36 
C 37 
C 38 
C 39 
C 40 
C 41 
C 42 
C 43 
C 44 
C 45 
C 46 
C 47 
C 48 
C 49 
C 50 


82 


OUTPUT 



7 

8 


9 

10 


1 1 
12 


13 

14 


15 

16 


6S/62X IH* * 2X H* » 6X H * ,2X1 H*/10X13H NO. UF S LOT S I 16 , 2 3< 1 H*, 2X IH*, 6X 1 
7H* i 2X 1H*/62X1H*,2X8H********«2XIH*/10XIIH SLOT PITCH,F18.3,1X6HINC 
8HES,12X2HHX,2XIH*,12XH*/62XLH*,2X8H********,2XIH*/ 10 X 11H SLOT PIT 
9CH, 41X IH* ,2XIH* ,6X1H* ,2XH*/10X15H AT 1/3 0 l ST . , F 14 . 3 , IX 6HI N CHES 

S. 16X1H*2X 1H*6XLH*2X1H*/57X6H * ,2 X8 H** ******, 2X l H*/ 62X1 H*, 12X H 

i*/62X 19H* ************* /62X1H1 , 12X1HI/62X14H1 8S 1/62X 

$1H1, 12XIH1) 

GO TO 13 

WRITE (6,8) B0,Bl,62,B3,BS,H0,HX,HT,HW,HS,lUQ,TS,TT 

FORMAT (1HL.13H STATOR SLU TS/ /5 X25 H T YPE-CONST A NT TOOTH WI0TH/51X1H 

11. 14X IH1/61X16H l 31 1 /10X3H BO , F2 5 . 3 , 1X5 H INCHES , 15X1 HI, l 

24X IHl / 10X3H 31»F26.3»22XH1,5X4H— B0-»5XIH1/IDX3H B2,F26.3, 1 IX 1 7H— 

3 l * , 4 XI 7H* 1 / 1 0 X3H B3 , F2 6 .3 , 2 2X l HI, 4< IH* 

4,4X1H*,4X1HI, 8X2HH0/10X15H 8S = (B2 + B3 ) / 2 , F 1 4. 3 , 22 X IHl , 4X l H*. 4X 17H 

5* 1 / 1 0X3H H0.F26.3 .2 2XIH1 ,2XIH*,8X1H*,2X IHl, 8X2HHT/1 

60X3H HX,F 26. 3,22X1H*, 14X, 1 2 H* /10X3H HT , F26 . 3, 12X2HHS , 7X 

71H* , 16X1H* , TX2HHW/10X3H HW.F26. 3 ,2 0X1 H* ,3 XI 2 H****** ******* 3X 10 H*— 

8 / 1 0X3H HS,F26. 3,l9X2H*l ,3 X1H*»10X1H*,3X2H1*/57X LH*, IX1H1,3< 

91H*. 10X1H*, 3X1H 1, IX1H* .4X2HHX/10X1 3H NO. OF SLOTS, 116, 17X l H*, 2X IHl 

$, 3X12H************»3X1H1,2XIH*»6H /5 5 XI H* , 3X1 HI , 18X 1H 1 , 3X l H*/ 

JIOXUH SLOT PITCH. F18.3.1 X6HI NCHES ,4X3 4H ************** ****** ** 

$********/ 54X IH l, 4X HI, 18X1 HI ,4 XI HI / 10X11 H SLOT P I T CH, 33X 1H l, 4X 20HI 

$ 32 1.4X1H1 /10X15H AT 1/3 01 S T . , F 1 4. 3 , IX 6H I NCH ES, 3 

tX IHl, 4X1H 1, 18X1H1 ,4X1H 1 / 54X30H1 B3 1/54X 1H 

28X IHl ) 

GO TO 13 

WRITE (6,10) BO,HO,BS»HS,IQQ»TS,TT 

FORMAT (1HL.13H STATOR SLOTS/ /5 X ,1 1 H T YPE-ROJ NO// l OX , 1 3H S LOT OPEN 
IING.F 16.3, 1X6HINCHES/10X.19H SLOT OPENING DE PT H , F l 0 . 3/ 10X , 14H S_3T 
2 DIAMETER .F15.3/10X11H SLOT DEPTH , FI 8 . 3/ / 10X , l 3 H NO. OF SLOTS, 116/ 
3/10X, 11H SLOT P ITCH ,F18. 3.1 X6HI NCHES//I0X.11H SLOT PI T CH/ LOX , l 5H 

4 AT 1/3 D I ST. , F 1 4. 3 , 1 X6HI NCHE S) 

GO TO 13 

WRITE (6,12) BS,HX»HY,HS,IQQ»TS»TT 

FORMAT (1HL.13H STATOR SLOTS//5X25H TYPE-OPEN (I CONO . / SLOT ) / 5 7X , 6 

1H *12X6H* /62X,1H*,12XIH*/53X5HHY *,12X1H*/62X IH*, 12X IH*/ 

210X.3H BS,F26.3,lX6HINCHES,ll X.6H *,2 X8H********, 2X IH*/ 10X , 3H 

3HX.F26.3, 23X, IH*, 2X1H * ,6X1H* , 2 XI H*/ 1 0 X , 3H HY , F26 . 3 , 23X , 1 H* , 2X l H*, 6 
4X1H*, 2XH*/10X, 3H H S, F 26. 3 ,23 X , 1 H* ,2 XIH* , 6 XI H* , 2X 1 H*, 2X 2HHS /58X 2HH 
5X. 2XIH*, 2X1H* , 6XIH*,2XIH*/10X,13H NO. OF S LOT S , II 5 . 23X 1 H*, 2X IH*, 6X 
61H*,2X1H*/62X1H*,2X1H*,6XH*,2X1H*/10X,11H SLOT P I T CH , F 1,8 . 3, IX 6H IN 

7CHES, 16X1H*, 2XIH* , 6XH* , 2 XIH* /5 7 X6H *,2X3 H********» 2XIH*/ 10X 11 

8H SLOT PITCH, 41 X1H* , 1 2X1H*/ 1 0 XI 5H AT 1/3 01 ST . , F 14. 3 , IX 6 HINCHES, 

916X 19H*** *********** /62X1HI .12X1H1/62X14H1 BS 1/6 2X IH 

$1, 12XIHI) 

CONTINUE 

WR £ TE ( 6» 14) GC , G2 , LG 2 » GE » CC 

FORMAT { IHL » 8H AIR GAP//lOX,13H MAIN AIR GAP , F23 .3 , IX , 6HINCHES/ 10X 
It 1 8H AUXILIARY AIR GA P f F l 8. 3/ 10X ,2 8H LENGTH OF AJXILIARY AIR GAP* F 
28 • 3/ i OXf 1 8H EFFECTIVE AIR GAP ,F 1 8- 3/ / 10 X. 19H CARTER COEFFI C IENT * FI 
37.3) 

IF (RF.LT..5) WRITE (6,15) 

IF (RF.GE..5) WRITE (6,16) 

FJRMAT ( IH I , 45H ARMATURE WINDING ( Y-CONNECTED , RANDOM WOJND)//) 
FORMAT ( IH 1 v 43H ARMATURE WINDING ( Y— CONNECTED, FORM WOJNDI//) 

IF (0W1.EQ.0.) WRITE (6,17) DW 


C 51 
C 52 
C 53 
C 54 
C 55 
C 56 
C 57 
C 58 
C 59 
C 60 
C 61 
C 62 
C 63 
C 64 
C 65 
C 66 
C 67 
C 68 
C 69 
C 70 
C 71 
C 72 
C 73 
C 74 
C 75 
C 76 
C 77 
C 78 
C 79 
C 80 
C 81 
C 82 
C 83 
C 84 
C 85 
C 86 
C 87 
C 88 
C 89 
C 90 
C 91 
C 92 
C 93 
C 94 
C 95 
C 96 
C 97 
C 98 
C 99 
C 100 
C 101 
C 102 
C 103 
C 104 
C 105 
C 106 


OUTPUT 
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I 



17 

18 


19 


20 


21 

22 


23 


24 

25 


26 


21 


28 


29 


FORMAT (IH * 9X* 16H STRAND DIANE TER*F32.4*1X6 HI NCHES) 

IF (DWi.GT.O.) WRITE (6*18) DW.DWl 

FORMAT <1H • 9X * 1 8H STRANO Dl ME NSl ONS • F 30. 4 , 2H X, IX * F6. 4 , 1 < * 6H I NCHE 
IS ) 

WRITE 16* 19) SO * SN * SNl • AC *S*CE *HM * EL * SK * RS, TST *RGL *STAT ET . EC,Y Y * GN 
FORMAT (IH * 1QX36HDI STANCE BTWN CL OF STRANDS ( RADI ALJ , F 1 1 . 4/ / 1 OX * 
13 8H STRANDS/CONDUCTOR IN RADIAL D I R. * F 1 1 . 0/10X ,24H TOfAL 5 X RANDS /C 
20 NDUCTOR* F20.0/ 10X* 1 5H CONDUCTOR ARE A* F33. 4 ♦ IX6HSQ-I N. / ID X * 29H CUR 
3RENT DENSITY AT FULL LOAD ,F 1 7. 2 *3X10H AMP/SQ-I N. / / 10X , 2 7H COIL EX T E 
4NSI0N BEYOND C0RE.F20.3, 2X6HI NCHE S/10X.24H MEAN LENGTH OF 1/2 TJRM 
5, F23.3/ 10X* 16H END TURN LENGTH ,F3 1* 3/ 10X . 1 7H STATOR SLOT SKEW * F30. 
63/ / 10X * 25H RESISTIVITY AT 20 DEG. C*F23. 4*1X1 6HMICR0 OHM INCHES/ID 
7X * 2 1 H STATOR RESISTANCE AT*F6.0,7H DEG. C *F 14. 4 * 1 X4H0HMS/ / i OX * 30H 
8N0. OF EFFECTIVE SERIES TURNS* F16. 2 /10X* 27H TOTAL EFFECTIVE CONDJC 
9T OR S * F 1 9. 2 / 1 OX • 1 4H SLOTS SPANNED *F 30. 0/10X.25H SLOTS PER POLE PER 
$PHASE*F21.2) 

WRITE (6.20) SC *C *PBA *FS*DF*CF 

FORMAT <1H *9X*16H CONOUC TORS/ SLOT *F28. 0 /10X * 2 5H NO. DF PARALLEL C 
IIRCUITS.F 19.0/10X* 1 7H PHASE BELT ANGLE , F2 7 .D . 5 X7HDEGREES/ / 10X , 1 2H 
2SKEW FACTOR * F35 .3/ 10X* 2 OH DISTRIBUTION F ACTOR ,F27. 3/ 10X , l 3 H PITCH 
3F ACTOR* F34. 3 I 

IF (RT.EU.O.) WRITE (6*21) RD 
IF CRT. GT. 0.1 WRITE (6*22) RD* RT 

FORMAT (1HL.14H FIELD WI NO I NGV /10X19H CONDUCTOR DI AMETE R, F 29 . 4 , 1X5 
1HINCHFS/) 

FORMAT (1HL*14H FIELD WI NO I NG/ /10X*2 l H CONDUCTOR D I MENS I ONS *F27. 4* 
12 H X* 1XF6.4* 1X6HINCHES/I 
WRITE (6*23) AS*NFC*PT*FE,RR*TF,FK1 

FORMAT (IH * 9X 1 5H CONDUCTOR ARE A *F33 . 4 * 1 X5HS0-I N. / / 1 OX , 1 9 H NO. OF 
1FIEL0 COILS* I24/10X*24H NO. OF TURNS (PER COIL) *F20.0/1DX20H MEAN 
2LENGTH OF TURN ,F27. 3.2X6HI NCHE S//10X2 5H RESISTIVITY AT 2D DEG. C*F 
323.4, IX16HMICR0 OHM I NCHE $/ 10X20H FIELD RESISTANCE AT*F5.D*7H DEG. 
4 C*F16.4* 1X4H0HMS) 

IF (NFC. ED. 2) WRITE (6,24) 

FORMAT (IH • 15X* 23H (BOTH COILS IN SERIES)) 

WRITE (6,25) PCO IL * DCOI L *BCOI L * TYP Y 

FORMAT ( 10X2 IH COIL INSIDE DI AME TE R*F26. 3.2X6HI NCHES/ 10X22H COIL 3 
10 TS IDE DIAMETER, F25.3/10X11H COIL WI D TH , F36 . 3 / / 10X , 2 1 H LOCATION CD 
2DE ( TYPY) • 123) 

WRITE (6*26) C 1 *CP , CM*CG *D1 

FORMAT ( 1HL * 10H C0NSTANTS//10X,35H Cl, FUNDAMENT AL/MAK. OF FIELD F 
ILUX.F8.3/ 10X.18H CP* POLE CONSTANT * F25. 3 / 10X • 2 7H CM, DEMAGNET 1 1 AT I 
20 N FACTOR, F 16. 3/10X, 31H CO* CROSS MAGNETIZATION FACT OR, F 1 2 . 3/ 10X , 2 
36 H 01, POLE FACE LOSS FAC TOR, F 17. 3 ) 

WRITE (6*27) D ( ,DU* CL * SS »HC *SF *HV • BV * BK • WL 

FORMAT (1H1.7H STATQR//10X *234 STATOR INSIDE DI AMETER* F 2 1. 2 • IX 6H I N 
1CHES/10X* 24H STATOR OUTSIDE DI AMETER *F 20 .2/10X ,20H OVERALL CORE IE 
2NGTH*F24.2/I0X*22H EFFECTIVE CORE LENGTH ,F22. 2/ LOX , 1 7H DEPTH BELOW 

3 SLOT,F27.2//10X* 16H STACKING FACTOR ,F28. 2/ /10X.21H NO. OF COOLING 

4 DUCTS. F2 U0/10X, 15H WIDTH OF DUCTS ,F29. 2 *1 X5HI NCHES// 1D< , 1 3H CORE 

5 LOSS AT* F6. 1* 1 7H KI LOLI NE S/SQ. I N. ,F 7. 1 * 2 X9HWATTS/ LB. ) 

IF (LTS.NE.O.) WRITE (6*28) LTS 

FORMAT ( 1 OX * 21 H LAMINATION THI CKNESS * F24. 3 *4H IN.) 

WRITE ( 6. 29) DASH( 1 ) *DASH ( 1 ) * DASH ( 1 ) , ( STAR( I ) * I * 1 . 5 ) * DASH( 1 ) * ( ST AR 
1( I). 1 = 1*4) *DASH( 1 1 * STAR! l)*ST4RIl> ,DASH ( 1 ) * ( DASH ( I ) *1 = 1*5) 

FORMAT (IHL.5H Y0KE///16X , 1H1 *22 X , 9HI l T Y3/ 16 X , 1H1 * A6 , 3 X , 3HT YL* 


C 107 
C 108 
C 109 
C 110 
C ill 
C 112 
C 113 
C 114 
C 115 
C 116 
C 117 
C 118 
C 119 
C 120 
C 121 
C 122 
C 123 
C 124 
C 125 
C 126 
C 127 
C 128 
C 129 
C 130 
C 131 
C 132 
C 133 
C 134 
C 135 
C 136 
C 137 
C 138 
C 139 
C 140 
C 141 
C 142 
C 143 
C 144 
C 145 
C 146 
C 147 
C 148 
C 149 
C 150 
C 151 
C 152 
C 153 
C 154 
C 155 
C 156 
C 157 
C 158 
C 159 
C 160 
C 161 
C 162 
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OUTPUT 



o o o 


30 


31 


32 


33 


34 


35 


36 


37 


14X * A 6 » 4H1 1, A3/16X,1H1 ,22X,4Hi 1 , 5X , 1H l / 43 X , 1H1/ 13X , 5 A6, IX, A6/13 

2X, LH* , 28 X y 1H* # 4X,2HTY/13X,1H*,28X,IH*/13X,1H*,2X,4A6,4H — AS/2( 1 
33X.4H* * « 22X f 4H* * , 5X , 1 H 1 / ) ,3 ( 12 X5 H* *,22X5H* */ ) , 2( 1 IX , 6H* 

4 *, 17X, 1H1, 4X,6H* */) , 1 1 X, 8H* ***,HX,12H *** VllK, 

58H* *, 14X,5HTYE * ,6 X , IH* /i 1 X, A6 , 2H** , 14X , 3H , IX, A6, 2H**,/ 11 

6X, lHi, 6X f 1H1 , 15X, H1,4X*1HI ,4X,1H1/5X»A5*8HI LG2 1 , AS , 10< , 18H 

71 1 TY2/11X.1H1 , 6X , 1 H 1 ,25X,1HL/1LX V 1H1* 2H — , 2 A6 , 4HT Y R , 2A6, 

82H — , LHI/ 1 IX, 1H1 .32X.1HL) 

WRITE ( 6, 30) TY , T Y2 , TY3 , TYL , T YE , TYK ,LG2,DYC 

FORMAT ( LHL,9X,2HTY»F30.2,7H £ NCHE S/ 1 OX ,3 HT Y2 «F29.2/ 13X, 3HTY3, F29. 
12/LOX, 3HTYL,F29.2/10X,3HTYE ,F29 . 2 / 1 OX , 3HT YR, F 29 . 2/ l OX , 3 HL G2, F29.2/ 
2/ 10X, 2 1HYGKE OUTSIDE D I AME TE R ,F i L. 2 ) 

WRITE <6, 31) DA SH< L) , SHL2 , D ASH ( 1 ) .STAR ( 1 ) . ST ARC 2 ) . ( DASHCI) . 1= L , 4), 
1ST ARC l)tSTAR{ 1 ) t STAR ( 1 ) ,STAR(1) ,DASH(1) , DR, 01 SHI 
FORMAT ( 1H1,6H RO TOR //i OX ,1 HI ,20X,lHi/10X,lHl,A6,F7.2,lX,A6, lHl/3( 
110X, 1H 1, 20X, 1 H 1 / ) IOXiHI ,4X,2 A6 ,2 X A5 / 1 3 X , 1 HI 3 X , A2 , L OX , A2 , 3X , L H 1, 
219X, 1HL/10X, 6H1 * *,10X,6H* * 1 , 1 9 X , H 1 / 2 X , 2 ( 10 X , 4H* *),2K,lHi 

3/ LX , 2 ( L0X,5H* *) ,20X,iHL/5X,A6,4X,lH*,10X,lH*,4X, A6 , LX,A6, 7X, 1H1 

4/2C 10X , LH*,4X, LH*) , 10X.H1 ,8X.1H1 /2 ( L0X,1H*,4X , IH*I , LOX, 1H1/2C LOX, 
51H*,4X, IH* ) , 10X,LH 1, = 1 L. 2 /2 ( 1 0X,1H* ,4X ,iH*)/2 < LOX, IH*, 4X , IH*), FI 3. 
62, 6X, 1H1/2C 10X,1H*,4X, H*) , l 9 X , 1 HL / 2 ( 1 OX , l H* , 4X , 1 H* ) , 2X , 2< 8X,H1)I 
WR I TE (6, 32) STAR <l) , STAR(l) ,DASHU> ,STAR< 1) .STARCH, < DAS-tCII, 1= 1, 
14), SHL 1 

FORMAT <5X,A6,4X, IH* , 1 OX , 1 H* , 4X , A6 , 1 X , A6 , 7 X , 1 H 1 / 1 1 X , 1 H*, 3* , IH*, LOK 
1, IH* * 3X,lH*,2OX»Hl/12X,lH*,2X»lH*,L0XtlH*,2X »1H*,21X» LHI / 13X, 1 H*, 
21X, IH* ,10X, IH*, IX, H* ,22X,1H1 /I X , 2 H** , 10 X ,2 H**,23X,LHi/15X, 2 AS, 2X 

3, 4A6//2( 15X,1H1»10X,1H1/) 1L X,5H L , F8 . 2 , 2 X , 5 H l /5X,2( 10K, LHI) 

4) 

WRITE <6,33) PW1,PU,PA1 ,PW2 ,PT2 ,PA2 , PL , QI SH , P E , VR , $ IGMA,SAMMAfR 
FORMAT ( 1HL , 9X , 25H POLE DIMENSIONS ( I NCHES ) / / 12X • L LH NARROW ENIO/LS 
LX , 6H WIDTH, F 1 8. 2 / 16 X, L OH TH l C KNESS , F 14 . 2 / 1 SX , 5 H AR E A, F 1 9. 2 / / l 2X , 9H 
2 WIDE END / 16X , 6H Wl D TH , F L 8. 2 / L 6 X ,L OH THICK NESS , F 1<* . 2/ L6X, 5H AREA, F 
31 9 • 2/ / 12X , 1 2H POLE LE NG TH , F 1 6 . 2 / 12 X ,2 2 H INSIDE SHAFT D I AME T ER/ l 4X , 
419H (FOR HOLLOW SH AF T ) , F 8 . 3 / / 1 0 X , L 3 H POLE EMBR ACE, F 17 . 2/ / i OX , L 7H 3 
5ERIPHERAL SP E E D , = 1 3 . 2 , 1 OH F E E T/ MI N ./ L 2 X , L 8H (AT MAIN AIR GAPJ/LOX, 
623H SPEC. TANGENTIAL F URGE ,F 7 . 2 , 1 1 H L BS/ S Q. I N. / / 1 OX ,2 4 H DENSITY OF 
7 NON-MAGNETIC /12X, 15H ROTOR M ATE RI A L , F 14 . 3 , 1 1 H LBS/CJ.IN.) 

WRITE (6, 34) WC , RC , Wl , WkOTOR,WYOKE , WTOT AL 

FORMAT (1HI,8H WE I 3 H T S/ /I 0X1 3 H STATOR CUNO . , F 17 . 3 , IX 6HP0JN DS / I OK 1 2 
1H FIELD OQND . ,F 18. 3/LOXI2H STATOR I RON ,F I 3 . 3 / l 0 X , 6 H ROT OR, F24. 3/ 10 
2X.5H YOKE ,F25.3//10X,6H TOTAL/L 1 XI 8H ( ELECTRUMAGNE T IC ) , F L 1 . 3 ) 

COMPUTE WINDAGE LOSS (MAIN AND AUXILIARY AIR GAPS) 

IF (WF— 1.0) 37,35,37 
WF 1 = 0 • 

C A - L WINDGt (WFl,RPM,:L,DR/2. ,GC,RE1) 

WF 2 = 1 .0 

CALL WINOGE ( WF 2 , RPM , LG2 , D I SH l /2 . , G 2 , RE2 ) 

WF = WF l+WF 2*2 . 

WRITE (6,36) RE1,RE2,WF1,WF2, WF 

FORMAT { 1 HK , 9X , 1 6H REYNOLDS N UM BER / l 3 X , 1 3 H MAIN AIR GAP, F17.2/ L 3 < , 
I13H AUX. AIR GAP, FI 7.2/lOX»21H WINDAGE LOSS ( W AT T 3 ) / 1 3X , 1 3 H MAIN A 
21 R GAP ,FI 7.2/ 13X, 13H AUX. AIR G AP , F I 7 . 2/ 1 3 X , 5 H T0TAL,F24.2) 
CONTINUE 

WRITE (6,38) PC,EW,SPS,Pl 


C 163 
C 164 
C 165 
C 166 
C 167 
C 168 
C 169 
C 170 
C 171 
C 172 
C 173 
C 174 
C 175 
C 176 
C 177 
C 178 
C 179 
C 180 
C 181 
C 182 
C 183 
C 184 
C 185 
C 186 
C 187 
C 188 
C 189 
C 190 
C 191 
C 192 
C 193 
C 194 
C 195 
C 196 
C 197 
C 198 
C 199 
C 200 
C 201 
C 202 
C 203 
C 204 
C 205 
C 206 
C 207 
C 208 
C 209 
C 210 
C 211 
C 212 
C 213 
C 214 
C 215 
C 216 
C 217 
C 218 
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IF ( (P2+P3+P4) .GT.1.0E-4) WRITE (6,39) P2,P3,P4 C 219 

WRITE 16,40) P5,P6,P7.P8.P9,P10 C 220 

38 FORMAT UHUllH PER MEANCE S//10 X, 1 6H WINDING LE AKAGE/ 13X , 2 6 H (PER I C 221 

INCH OF CORE LENGTH )/ 12X , 12H STATOR SLOT ,F26. 3 , I8H LI NES/ AMPERE T JR C 222 

2N,/13X, 10HSTAT0R END ,F27. 3/13X ,13H SPEC I AL ( SPS ) , F24. 3// 10K , 14H PI C 223 
3( POLE T IP ) , F 26. 3 ) C 224 

39 FORMAT U0X,3H P2. F 37. 3/ 1 OX , 3H P3 , F37. 3 / 10 X , 3H P4.F37.3) C 225 

40 FORMAT (10X.3H P5. F 37. 3/1 0X.3H P6 , F37. 3/ 10 X, 3H P7 , F37. 3/ ID X , 18 H P8 C 226 

1 (MAIM AIR GAP ) .F22.3/10X, 18H P9 (AUX. AIR G AP ) , F2 2. 3/ 10 K , 4 H P10.F C 227 
236.3) C 228 

IF IWF1.GT. -0.0001) WRITE (6,41) C 229 

41 FORMAT ( 1HI ) C 230 

WRITE (6.42) A.XR,XL.XD,XQ.XA,XB,XF,XU.SI ,TF,TC C 231 

42 FORMAT (IHL.UH RE AC TANCE S//1 0X23H AMPERE CONDUCTORS/ I NCH, F 20 .3/ 10 C 232 

IX 17H REACTANCE FAC TOR ,F26. 3//1 0X23H STATOR WINDING LEAKAGE, F20 .3, 1 C 233 
2X.7HPERCENT/10X23H ARM. REACTION ( 01 RECT ) ,F2 0. 3/ 10 X22H ARM. REACTI C 234 
30 N (OUA0.1.F21.3/10X21H SYNCHRONOUS (DIRECT) ,F22. 3 / 10X20H SYNCHROS C 235 
40 US OUAD.),F23.3/10X14H FIELO LEAKAGE ,F29. 3/ 10X10H TRANS I ENT , F33. C 236 
53//10X22H FIELD SELF I NDUCTANC E ,F2 1. 3 , 1 X7 HHENRI ES/ / 10X 27H OPEN CIR C 237 
6CUIT TIME CONSTANT / 13X, 16H (FIELO ONLY, AT.F4.0.8H DEG. C).F15.5,1 C 238 
7X7HSECUNDS) C 239 

RFTURN C 240 

END C 241- 
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SUBROUTINE MAGNET 

COMMON /COM 3/ ACCH E * A I * A SHE T f A TOOTH « A YKl * B AG2 * BC L, BPL* BROT * B$ HFT * B 
IT Lt BYCL 1, B YCL2 *DC3RE*DSHMID*DYK1*DYK2 *FCL*FFL .FGLi ,FGL2 * FSXL , FPL* F 
2ROT* FSHFT *FTL* F YCL1 *FYCL2 *KSAT *PA A » I TRTN 
COMMON /COM 5/ A A • A G2 *CL . C P * OY2 • DY3 ,HS,LG2,P1 * P10 * P5 * PS * P7 * P8* P9 * PA 
1 2 * PPL * P PL 1*PPL2*PPL3*PPL4*PX*TY2*TY3*TYPY 

INTEGER TYPY 

REAL LG2 

DIMENSION A I 190 ) * PAA(IO)* ATKIO), AT2I10)* AT3U01. FL JK 1(13)* F 
11 UX 2 < 101* FLUX3U0)* ERRORi (10) , ERRQR3 I 13 ) 

CGUNT=0. 

BYCL 1=0- 
BYCL 2=0- 
B PL =0- 
BTL = 0- 
ft CL = 0- 
F FL = 0- 
F YCL 1 = 0- 
F YCL 2=0 • 

FPL = 0- 
F f.L =0- 
BROT = 0 
FROT = 0 
BSHFT=0. 

FSHF T=0. 

F TL = 0. 

00 2 1 = 1 * 10 
ERR0R3< U = l .0E10 
Ft UX i< U = PPL/10- 


1=0 

COUNT=COUNT+1-0 

1=1 + 1 

PPM = FL UX1 ( I I 

CALCULATE MAIN AIR-GAP AMPERE-TURNS 

F GL 1 = 1 ( IPPM/CP )*PX)/P8)*i-0E4 

FLUX OENSITY AND AMPERE-TURNS FOR TEETH 

B TL= ( PPM/ A TOOTH > *10* 

X =8TL 
N A= 1 
X =1 


D l 
D 2 
D 3 
D 4 
D 5 
0 6 
D 7 
D 8 
0 9 

D 10 
0 11 
0 12 
0 13 

0 14 

0 15 

D 16 
0 17 

0 18 
D 19 
D 20 
D 21 
D 22 
D 23 
0 24 

0 25 

D 26 
D 27 
0 28 
D 29 
0 30 

0 31 

D 32 
0 33 

D 34 
D 35 
D 36 
0 37 

D 38 
D 39 
D 40 
D 41 
D 42 
D 43 
0 44 

D 45 
D 46 
0 47 

D 48 
D 49 
0 50 
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KSA T= 3 0 51 

G3 T3 32 D 52 

FTL=AT*HS 0 53 

D 54 

FLUX DENSITY AND AMPERE-TURNS FOR CORE D 55 

D 56 

BCL=(PPM/(AL0RE*2.) )*10. D 57 

X =BCL D 58 

N A= I D 59 

KSA T= 4 D 60 

K = 2 D 61 

GD TU 32 D 62 

FCL=AT*DC0RE D 63 

D 64 

AT3( 1 )=(FGL l + FTL+FSL+FSXL/2. I *2. 0 65 

IF ( I .LT. 10) GO TO 4 D 66 

D 67 

POLE-TIP LEAKAGE FLUX D 68 

D 69 

PPL 1 1 = AT 3 ( 1)*0.001*P1 D 70 

D3 7 1 = 1, 10 0 71 

FLUX2I I ) =A T3 ( I ) *P 1 0*0. 000 1 D 72 

IF (I.EQ.l) FLUX3I 1) =PPL1 1 + ( F LUX1 ( 1 ) +FLJX2 ( 1) ) / 2 . 0 73 

IF (I.GT.l) FLUX3II) =FLUX3 ( I — 1 > + IFLUX1 1 1 -1 )+FLJX2( 1-1) +FLJK U I )+F. D 74 

lUX2(I>)/2. D 75 

0 76 

FLUX DENSITIES AND AMPERE TURNS FOR POLE D 77 

D 78 

K= 3 D 79 

KSA T= 5 D 80 

BPL =0 . 0 81 

1=0 D 82 

8 1=1+1 D 83 

N A = 31 D 84 

X = FLUX 31 I )/PAA(I) 0 85 

BPL=AMAX1IBPL,X) D 86 

G3 TO 32 D 87 

9 AT 1 ( I )=AT*CL/10. D 88 

IF (I.GT.l) AT1I I ) =AT1( I >+ATl (I -1) D 89 

IF < I .LT. 10) GO TO 8 D 90 

FPL =A T II 1 0) D 91 

FBB=FP L+A T3( 1 ) -0. 5* (AT1(1 )+ATl (9)-ATl (10) ) D 92 

IF ( ITRTN.NE. 10) GO TO 18 D 93 

D3 10 1=1,10 0 94 

11=11-1 D 95 

AT2( I ) =FB B-2.*FPL+0.5* (ATI (I ) +AT1 ( I I ) I 0 96 

IF (I.NE.IO) AT2(I) =A T2 ( I ) +( ATI ( II -1 ) ) *0.5 0 97 

10 IF (I.NE.l) AT2< I) =AT2 (I ) + ( ATI ( 1-1 ) ) *0.5 0 98 

C D 99 

£RR0R=0. D 100 

D3 11 1=1,10 D 101 

ERROR 11 I)=AT2( I)-AT3(I) D 102 

11 ERROR = ERR 0R+A3S(ERR0R1(I)) 0 103 

DO 14 1=2,9 D 104 

AA=ERR0R3 ( I ) -ERROR! ( I ) D 105 

IF (ABS(AA)-1.0E-8) 12,12,13 0 106 
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I 



12 

AA = 1.0 


D 

107 


GO TO 14 


D 

108 

13 

AA=ERR0R3 ( I ) /A A 


D 

109 

14 

FLUX L I I}=((AT2(I)/AT3(I))* 5 MAA) ) *FLUX1 U ) 


D 

110 


DO 15 1=2 , 9 


D 

111 

15 

ERROR 3(11 =ERROR 1(1 ) 


D 

112 


FL U X= 0 . 


D 

113 


DO 16 1=1,10 


D 

114 

16 

FLUX=FLUX+FLUX1 (I ) 


D 

115 


DO 17 1=1,10 


D 

116 

17 

FL UX L { i)=(FLUXl(l) )*(PPL/FLUX> 


D 

l 17 


IF ( ERROR .LT. 0* 56-3 ) GO TO 18 


D 

118 


IF {COUNT .LE. 24. 9) GO TO 3 


D 

119 


ITR TN = 0 


D 

L20 


GO TO 1 


D 

121 

18 

PPL INFLUX 3 { 10) + (-LUX1 ( 10) + FLUX2 (10) ) / 2. 


D 

122 

C 



D 

123 

C 

TOTAL FLUX THROUGH TAPERED ROTOR SECTION 


D 

124 

C 



D 

125 


PPL 2=(PPL l*PX/2. ) + (AT3tl) /2. ) *(O.OOL»* <P7)+PPL11*< PX/2. ) 

D 

126 

C 



D 

127 

C 

FLUX DENSITY AND AMPERE-TURNS FOR TAPERED SECTION OF ROTOR 

D 

128 

c 



D 

129 


K= 7 


D 

130 


tn 

o 

• 

o 

II 

z 

Q 


D 

13! 


KS A T= 2 


0 

132 


ARM ID=PA2^ (PX/2. ) 


D 

133 

19 

DM = 0M + 0.1 


0 

134 


NJA=31 


D 

135 


AREA=ARM1 D-( { ARM I D- A SHF T ) *DN ) 


D 

136 


X=(PPL2)/(ARtA) 


D 

137 


BROT=AMAX 1 (X, 3R0 1 ) 


n 

138 


GO TO 32 


D 

139 

20 

FR JT=FfOT+AT*L>SHMID /I 0. 


D 

140 


IF 10N.LT.U.9) GJ TO 19 


D 

141 

C 



D 

L 4? 

C 

TOTAL FLUX THRuUGH ROTOR END SECTION 


D 

143 

C 



D 

144 


PPL 3=PPL2 + (F8ti + 2.*FR0T)*P6*0. 001 


D 

145 


3 Pl X = P PL 3 


D 

146 

c 



D 

147 

c 

FLUX DENS ITT AND AMPERE-1 URNS FOR AUXILIARY 

AIR-GAP 

D 

148 

c 



D 

149 


8AG2=PPL3/A32 


D 

150 


FGL 2=P PL 3 / ( p 9* 0 • 001 ) 


D 

151 

c 



0 

152 

c 

FLUX DENSITY AND AMPERE-TURNS FOR ROTOR END 

SECT ION 

D 

153 

c 



D 

154 


K = 6 


D 

155 


DM =-0 • 05 


D 

156 


KS AT= 1 


D 

157 

21 

DM =DM +0. 1 


D 

158 


M A = 3 1 


D 

159 


X ~P PL 3*DN 


D 

160 


X=X/ASHFT 


D 

161 


6SHFT = AMA Xl(X,bSHFT) 


D 

162 
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GO TO 32 


D 

163 

22 

FSHFT=FSHFT+AT*LG2/10. 


D 

164 


IF IDN.LT.O. 90) GU TO 21 


D 

165 

C 



0 

166 

C 

FLUX DENSITIES AND AM PE Kb- TUR NS FUR END SECTION OF 

YOKE 

D 

167 

c 



D 

168 


<K=1 


0 

169 

23 

IF ( T YPY • LE • 4) GO TO (27,24), KK 


D 

170 

24 

DM =-0.05 


D 

171 


K SA T= 6 


D 

172 


K = 4 


D 

173 


DY< 2= ( DY3-DY2 ) /20. 


D 

174 

25 

ON =DN + 0.1 


D 

175 


AR E A= 3 . 14 1 bM D Y 2+ ( 2 0. * D YK2 *DN H 


n 

176 


AR b A= AREA* ( TY3+( ( TY2-TY3) *< 1. O-DN) ) ) 


D 

177 


X=?PL X/AREA 


D 

178 


BYCL2 = AMAXKX,8Y;L2) 


D 

179 


N A = 6 1 


D 

L 80 


GO TO 32 


D 

181 

26 

FYCL 2= AT* (DYK2J+F YCL2 


D 

182 


IF (DN .Lf.0.9) GO TO 25 


0 

183 

C 



D 

184 

C 

Flux DtMSITY ANO A rtPe RE - T UR NS FOR CYLINDRICAL POKF [ON OF YOKE 

□ 

185 

c 



D 

186 

27 

IF (TYPY.NE.l) GJ TO (30,28) ,KK 


D 

187 

28 

BYCL 1=PPL X/AYKI 


D 

188 


X = 3 YC L i 


D 

189 


.M A = 6 1 


D 

190 


KSAT= 1 


0 

191 


.< - 5 


0 

192 


GO TO 32 


D 

193 

29 

FYCL 1 = AT* D YK 1 


0 

194 


IF (KK.E0.2) GO TO 31 


n 

195 

C 



n 

196 

C 

FLUX THROUGH SECTION OF YOKE AFFECTED BY PERMEANCE 

P5 

D 

197 

c 



0 

198 

30 

AT = FBB +FGL2*FSHFT + FR0 1 


D 

199 


IF (TYPY.GT.4) A T=A T+ F YC L2 


D 

200 


IF (TYPY.EQ.6) AT=2.*AT 


D 

20 l 


IF (TYPY.EQ.i) AT=AT«-FYCLi 


D 

202 


PPL 4= A T*P 5*0. 00L+PPL3 


D 

203 


PPL X=P PL 4 


0 

204 


K K = 2 


0 

205 


GO TO 23 


D 

206 

C 



D 

207 

c 

TOTAL AMPERE-TURNS 


D 

208 

c 



0 

209 

31 

FFL =( FYCL 2+FBB/2.+FGL24-F YCLU-F SHFT+FROT) *2. 


D 

210 


< S A T= 1 0 


D 

211 


RETURN 


D 

212 

C 



D 

213 

c 

INTERPOLATION PROCEDURE FOR MATERIAL CURVES 


D 

214 

c 



D 

215 

32 

I r UI(NA>.LT. X) GO TO 36 


D 

216 


NA=NA+3 


D 

217 

33 

IF 1 A I (NA i — X J 34*35 *35 


D 

218 

34 

NA=NA+2 


D 

219 


GO TO 33 


D 

220 

35 

XX= I A I (NA ) — A I i NA— 2 > 1 /{ ALOG I All NA+1 ) / { AI( NA-1 ) +0. 3 33 1 ) ) J 

D 

221 


Y =A l < N A ) - XX* A LOG I A I I NA+1) I 


D 

222 


AT=EXPUX-Y)/XX) 


D 

223 


GO TO 15*6*9*26*29*22 *20) *K 


D 

224 

36 

IF (ITRTN-EO.O) RETURN 


D 

225 


IF I KSAT . LT.3.0R. KSAT.GT .5) RETURN 


D 

226 


ITRTN^O 


D 

227 


GO TO 1 


D 

228 


END 


D 

229 


MAGNET 
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Calculate PPL2, total 
flux through tapered 
rotor section 


19-20 


Calculate flux density 
BROT and ampere-turn 
drop FROT for tapered 
rotor section 


Calculate PPL3, total 
flux through rotor 
end section 


PPLX « PPL3 


Calculate BAG2 and FGL2 


21-22 


{See detail C) 


Calculate flux density BSHFT 
and ampere-turn drop FSFIFT 
for rotor end section 


(See Detail D) 



24-26 

Calculate flux density BYCL2 
and ampere-turn drop FYCL2 
for yoke end section 


(Similar to 
detail C) 



Calculate flux density BYCL1 
and ampere-turn drop FYCL1 
for cylindrical section of yoke 



(Similar to 
detail A) 


30 


Calculate leakage flux through 
permeance P5 and total flux 
PPL4 through section of yoke 
affected by P5 


Calculate total ampere- 
turn drop FFL 


i 

KSA' 

l 

r - 10 


i 

( RETURN ) 


i 

PPLX = 

i 

- PPL4 

j 


| KK 

1 

= 2 

i 



Detail a 



Detail B 


MAGNET 


93 




















Detail E 


MAGNET 











c 


SUBROUTINE WINDGE ( P .RPM f L #R*T .RE I 


C 

C 


l 


2 

3 

4 


5 


REAL L • NU • M 

NAMELIST /WIND/ RO , VI S . PRE S S. M ♦ TEMP 

IE (P.NE.O.) GO TO 3 
Rn=o.d 
PRESS- 14. 7 
TEMP=25. 

V IS=1.24E-5 

N =0. 

READ <5. WIND) 

IF (R3.EQ.0.0) RQ*(PRESS*M*5. 16E-02) /(TEM P+273.0) 

WRITE ( 6. U VIS. R3. PRESS. TEMP 
IE (M.NE.O.) WRITE (6,21 M 

FORMAT I1HL.8H WINDAGE //10X.17H FLUID PROPERT I ES/ 1 3X , 1 3 H VISCOSITY 
1. 1PE25.3. 11H L8S/SEC-FT/13 X.83 DE NS I T V .F24. 4 ,1 1 H L BS/CU. Ff . / 1 3X , 9H 
? PRESSURE. 0PF21. 2. 11H LBS/SO. IN/I3X.I2H TEMPE RATU RE . F 1 9 . 2 . 7H DEG. 
3 C) 

FORMAT (IH .12X.17H MOLECULAR WE I 3HT . E 1 2 . I ) 

NU=< VI S/ROI* 144. 0 
R E=3. 14 159*R*T*RPM/{ 30. 0*NU) 

C DN = I . 0 
C CO = CON 

CONa 1.0/1 2- 04+ 1 . 768* A LOG < RE*SQRTI CDOi ) 1**2 
IF < ABSICDN-CDCU/CON.LE. 1.0E-04) GO TO 5 
GO TO 4 

P = < 5.85E- 10)*CDN*RQ*(RPM**3)*I R**4)*L 

R ETURN 

END 


E 1 
E 2 
E 3 
E 4 
E 5 
E 6 
E 7 
E 8 
E 9 
E 10 
E 11 
E 12 
E 13 
E 14 
E 15 
E 16 
E 17 
E 18 
E 19 
E 20 
E 21 
E 22 
E 23 
E 24 
E 25 
E 26 
E 27 
E 28 
E 29 
E 30 
E 31- 


WINDGE 
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APPENDIX D 


DEFINITION OF FORTRAN VARIABLES 

Following are two alphabetic lists of nearly all FORTRAN variables used in the pro- 
gram. The first list gives the variables used in LUNDEL, SUBLUN, OUTPUT, and 
MAGNET. The second list gives the variables used in WENDGE. Each variable is de- 
fined and the units used in the program for each variable are given. Variables that 
are subscripted (dimensioned) or that are permissible program input are identified by 
a + or *, respectively. 


Typical straight iine segment 
between input data points on 



Log of magnetizing force, ampere-tur ns/in. 

Figure 27. - Variables used in interpolation 
between points on magnetization curves. 


FORTRAN SYMBOLS USED IN LUNDEL, SUBLUN, OUTPUT, AND MAGNET 


A 

AA 

AB 

AC 

A CORE 
AG2 

*+ AI 

AIRGAP 


ampere conductors per inch of stator periphery, A/in. 

dummy variable used for variety of calculations 

dummy variable used for variety of calculations 

2 

armature conductor area, in. 

2 

core cross-sectional area used in magnetic calculations, in. 
auxiliary air-gap area, in. ^ 

coordinates of points on material magnetization curves (see figs. 
20 and 27) 

NAMELIST name 
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+ A KVA 

* ALPHAR 

* ALPHAS 
AN 
AREA 
ARMID 
AS 

ASHFT 

AT 

+ ATI 

+ AT2 

+ AT3 

ATOOTH 

AYK1 

B1 

B2 

* B3 
BAG 2 
BB 
BB 
BCL 

+ BCLL 

* BCOIL 
BG 

* BK 

* BO 
BPL 

+ BPLL 


alternator output at load point G, kVA 
temperature coefficient of RR, °C * 
temperature coefficient of RS, °C - ^ 
power factor angle (see fig. 7), rad 
dummy variable (generally an area) 

=PA2*(PX/2.), in. 2 

o 

field conductor area, in. 

2 

cross-sectional area of rotor under auxiliary air gap, in. 

ampere-turns per inch at flux density X for magnetic materials 
(see fig. 27), A-turn/in. 

ampere-turn drop from pole tip to i pole segment 
A-turns 

ampere-turn drop across teeth, core, and main air 

ampere-turn drop across teeth, core, and main air 

2 

tooth area used in magnetic calculations, in. 

2 

yoke area used in magnetic calculations, in. 
stator slot dimension (see fig. 21), in. 
stator slot dimension (see fig. 21), in. 
stator slot dimension (see fig. 21), in. 

2 

auxiliary air-gap flux density, kilolines/in. 
displacement angle (see fig. 7), rad 
dummy variable 

2 

core flux density, kilolines/in. 

2 

core flux density, kilolines/in. 
field coil width, in. 

2 

main air-gap flux density (no load rated voltage), kilolines/in. 

2 

flux density at which core loss WL is specified, kilolines/in. 

stator slot dimension (see fig. 21), in. 

2 

pole flux density, kilo lines/in. 

2 

pole flux density, kilolines/in. 


(see fig. 9), 

gap, A-turn 
gap, A-turn 
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BROT 

+ BROTL 

* BS 
BSHFT 

+ BSHFTL 

BT 

BTL 

+ BTLL 

* BV 
BYCL1 
BYCL2 

+ BYCLL1 

+ BYCLL2 

* C 

* Cl 
CC 

+ CDD 

* CE 
CF 
CK 

* CL 

* CM 

+ CMPNT 

COM1 
COM2 
COM3 
COM4 
COM 5 
COUNT 


p 

flux density in tapered section of rotor, kilolines/in. 

p 

flux density in tapered section of rotor, kilolines/in. 
stator slot dimension (see fig. 21), in. 

2 

shaft flux density (under auxiliary air gap), kilolines/in. 

2 

shaft flux density (under auxiliary air gap), kilolines/in. 

dummy variable used in slot and end turn leakage permeance 
calculations 

2 

tooth flux density, kilolines/in. 

2 

tooth flux density, kilolines/in. 
width of cooling ducts, in. 

2 

yoke flux density (cylindrical section), kilolines/in. 

2 

yoke flux density (end section), kilolines/in. 

2 

yoke flux density (cylindrical section), kilolines/in. 

2 

yoke flux density (end section), kilolines/in. 
number of parallel armature circuits per phase 
ratio of fundamental to actual value of field form 
Carter coefficient 

2 

current density in field winding, A/in. 
straight portion of coil extension, in. 
pitch factor 

power-factor adjustment factor 

length of stator stack (axial direction), in. 

demagnetizing factor (direct axis) 

array of names used for saturation message print-out 

COMMON block name 

COMMON block name 

COMMON block name 

COMMON block name 

COMMON block name 

counts number of iteration cycles in MAGNET 
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* 

CP 

ratio of average to maximum value of field form 

* 

CQ 

cross-magnetizing factor (quadrature axis) 


CS 

per unit pole pitch 


CW 

winding constant 


CX 

dummy variable used in slot leakage permeance calculation 


CY 

2 

used to calculate area correction D, in. 


CZ 

2 

used to calculate area correction D, in. 


D 

used in distribution factor calculation 


D 

area correction for corner radii in rectangular conductor, in. 

* 

. D1 

pole-face loss factor 


D2 

factor used in pole-face loss calculation 


D3 

factor used in pole-face loss calculation 


D4 

factor used in pole-face loss calculation 


D5 

factor used in pole-face loss calculation 


D6 

factor used in pole-face loss calculation 

+ 

DA 

2 

used to calculate area correction, in. 

+ 

DASH 

used in OUTPUT to draw diagrams 


DB 

diameter of bender pin for forming armature coils, in. 

* 

DCOIL 

field coil outside diameter, in. 


DCORE 

core length used in magnetic calculations, in. 


DF 

distribution factor 


DG 

maximum dimension of rectangular conductor, in. 

* 

DI 

stator inside diameter, in. 

* 

DISH 

shaft inside diameter (of hollow shaft) (see fig. 23), in. 

* 

DISH1 

rotor diameter under auxiliary air gap (see fig. 23), in. 


DN 

used in magnetic calculation for rotor and yoke 


DR 

rotor diameter at main air gap (see fig. 23), in. 


DSHMID 

(SHL2-CL)/2. (see fig. 23), in. 


DT 

minimum dimension of rectangular conductor, in. 
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* 

DU 

stator outside diameter, in. 

* 

DW 

armature strand diameter or width (see fig. 24), in. 

* 

DW1 

armature strand thickness (uninsulated) (see fig. 24), in. 

+ 

DX 

2 

used in rectangular conductor area calculation, in. 

+ 

DY 

2 

used in rectangular conductor area calculation, in. 


DY2 

yoke dimension (see fig. 22), in. 


DY3 

yoke dimension (see fig. 22), in. 

* 

DYC 

yoke outside diameter (see fig. 22), in. 


DYK1 

yoke dimension used in magnetic calculations, in. 


DYK2 

1/10 of yoke end section length (see fig. 22), in. 

+ 

DZ 

2 

used in rectangular conductor area calculation, in. 

+ 

E 

alternator electromagnetic efficiency, percent 

+ 

El 

alternator overall efficiency, percent 


EB 

eddy factor (bottom) 


EC 

number of effective armature conductors 


EDD 

excitation voltage, per unit 

* 

EE 

line-to-line rated alternator voltage, V rms 

+ 

EF 

field voltage, V dc 


EK 

leakage reactive factor for stator end turns 

* 

EL 

end extension length of armature coil, in. 

* 

EP 

line-to-neutral rated alternator voltage, V rms 


ERROR 

ERROR1 (l)+ERRORl (2)+ . . . +ERROR(10), A-turn 

+ 

ERROR1 

error in interation of air-gap flux distribution, A-turn 

+ 

ERROR3 

values of ERRORl from pervious iteration step, A-turn 


ET 

eddy factor (top) 


EW 

specific stator end winding leakage permeance, lines/A-turn 

+ 

EX 

eddy loss, W 

+ 

EZ 

eddy factor (average) 

* 

F 

frequency, Hz 
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FACTOR 

dummy variable in slot leakage permeance calculations 


FBB 

ampere-turn drop from pole base to pole base, A-turn 


FCL 

core ampere turns, A-turn 

+ 

FCLL 

core ampere turns, A-turn 


FE 

mean length of one field coil turn, in. 


FF 

dummy variable in slot leakage permeance calculation 


FFL 

total ampere turns, A-turn 

+ 

FFLL 

total ampere turns, A-turn 


FG1 

main air-gap ampere-turn drop (uniform flux distribution), A-turn 


FG2 

auxiliary air-gap ampere-turn drop (negligible leakage flux), 
A-turn 


FGL1 

main air-gap ampere turns, A-turn 


FGL2 

auxiliary air-gap ampere turns, A-turn 

+ 

FGLL1 

main air-gap ampere turns, A-turn 

+ 

FGLL2 

auxiliary air-gap ampere turns, A-turn 


FGML 

demagnetizing ampere-turns per pole at rated load, A-turn 


FGXL 

demagnetizing ampere-turns at load G (see table m), A-turn 


FH 

air-gap ampere-turns at no-load rated voltage, A-turn 

+ 

FI 

field current, A 


FIELD 

NAMELIST name 


FIMM 

field current at rated load, A 


FK1 

field winding resistance at temperature TF, ohm 


FLUX 

total air-gap flux after calculation new flux distribution, kilolines 

+ 

FLUX1 

useful air-gap flux distribution (see fig. 9), kilolines 

4 - 

FLUX2 

flux through rotor leakage permeance P10 (see fig. 9), kilolines 

+ 

FLUX3 

flux at center of each of 10 pole segments (see fig. 9), kilolines 


FNQ 

used in distribution factor calculation 


FPL 

pole ampere turns, A-turn 

+ 

FPLL 

pole ampere turns, A-turn 
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FQ 

FQLL 

FROT 

FROTL 

FS 

FSC 

FSHFT 

FSHFTL 

FTL 

FTLL 

FYCL1 

FYCL2 

FYCLL1 

FYCLL2 

G 

G2 

GA 

GAMMAR 

GC 

GE 

GF 

GT 

GXX 

HC 

HM 

HO 

HS 

HT 


flux per pole (no-load rated voltage), kilolines 
flux per pole at no-load voltage QPERV, kilolines 
ampere-turn drop across tapered rotor section, A-turn 
ampere-turn drop across tapered rotor section, A-turn 
skew factor 

short-circuit ampere turns for rated armature current, A-turn 

ampere turn drop across rotor end section, A-turn 

ampere turns across portion of shaft under auxiliary air gap, 
A-turn 

tooth ampere turns, A-turn 
tooth ampere turns, A-turn 

yoke (cylindrical section) ampere turns, A-turn 
yoke (end section) ampere turns, A-turn 
yoke (cylindrical section) ampere turns, A-turn 
yoke (end section) ampere turns, A-turn 

volt-ampere alternator output at which load characteristics are 
calculated, per unit or percent 

auxiliary air gap (radial length) (see fig. 23), in. 

main air-gap area, in. ^ 

3 

density of nonmagnetic rotor material, lb/in. 
main air gap (radial length) (see fig. 23), in. 
effective main air gap, in. 

constant used in load pole-face loss calculation 
ratio of slot opening width to main air gap 
flux per pole multiplying factor (see fig. 7) 
stator depth below slot, in. 

armature conductor length (=1/2 coil length), in. 
stator slot dimension (see fig. 21), in. 
stator slot dimension (see fig. 21), in. 
stator slot dimension (see fig. 21), in. 


HV 

HW 

HX 

HY 

EDELR 

EDM 

ETQQ 

IPN 

IPX 

IQQ 

ITRTN 

IZY 

IZZ 

JA 

KSAT 

LG2 

LTS 

LUNDEL 

MAGNET 

MM 

NFC 

NPF 

OUTPUT 

PI 

P10 

P2 

P3 

P4 

P5 


number of stator cooling ducts 

stator slot dimension (see fig. 21), in. 

stator slot dimension (see fig. 21), in. 

stator slot dimension (see fig. 21), in. 

increment by which QPERV is increased, percent 

(Integer) * (Number of poles) * (Number of phases) 

used to calculate DF for fractional slot winding 

number of phases 

number of poles 

number of stator slots 

if ITRTN * 10 uniform air-gap flux distribution was assumed 
(Number of poles) * (Number of phases) 
type of stator slot (see fig. 21) 

number of load calculations made before alternator saturated 
saturation indicator (if KSAT + 10 alternator is saturated) 
auxiliary air-gap length (axial direction) (see fig. 23), in. 
stator lamination thickness, in. 
main program name 
subroutine name 

subscript of load point G such that G(MM) = 1.0 
number of field coils 

number of power factors for which load calculations are made 
subroutine name 

pole tip leakage permeance (see fig. 6), lines/A-turn 
leakage permeance = 2. * (P3 + P2 + P4), lines/A-turn 
leakage permeance (see fig. 6), lines/A-turn 
leakage permeance (see fig. 6), lines/A-turn 
leakage permeance (see fig. 6), lines/A-turn 
leakage permeance (see figs. 5 and 8), lines/A-turn 



* P6 

* P7 
P8 
P9 

*+ PA 

* PA1 

* PA2 

+ PAA 

* PBA 
PC 

* PCOIL 

* PE 
PERMCE 
PF 

PF 

*+ PFC 

PHW 
PI 

* PL 
PN 
POL 

+ PP 

PPL 
PPL1 
PPL11 

PPL2 

PPL3 

PPL4 


leakage permeance (see fig. 6), lines/A-turn 
leakage permeance (see figs. 5 and 8), lines/A-turn 
main air-gap permeance, lines/A-turn 
auxiliary air-gap permeance, lines/A-turn 

2 

area at ends of each of 10 pole segments (see fig. 23), in. 

2 

pole area at tip of pole (see fig. 23), in. 

2 

pole area at base of pole (see fig. 23), in, 

2 

area at center of each of 10 pole segments, in. 
phase belt angle, deg 

specific stator slot winding leakage permeance, lines/A-turn 

field coil inside diameter, in. 

pole embrace 

NAMELIST name 

design power factor 

power factor at which load characteristics are calculated 

power factors at which load characteristics are calculated 

average pole width, in. 

rated line current, A 

pole length (see fig. 23), in. 

number of phases 

dummy variable 

pole-face losses, W 

flux per pole in main air gap, core, and teeth, kilolines 

total flux leaving base of pole (see figs. 8 to 10), kilolines 

pole-tip leakage flux through permeance Pi (see figs. 9 and 10), 
kilolines 

flux in tapered rotor section (see fig. 10), kilolines 

flux in rotor end section, auxiliary air gap and yoke (see fig. 10), 
kilolines 

flux in yoke including leakage through P5, kilolines 
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PPLX 

dummy variable (PPLX = PPL3 or PPLX = PPL4), kilolines 


PPM 

=FLUX1 (I) 

+ 

PR 

field losses, W 

+ 

PS 

stator copper loss, W 

* 

PT 

number of field turns per field coil 

* 

PT1 

pole thickness at pole tip (see fig. 23), in. 

* 

PT2 

pole thickness at pole base (see fig. 23), in. 

* 

PW1 

pole width at pole tip (see fig. 23), in. 

* 

PW2 

pole width at pole base (see fig. 23), in. 


PX 

number of poles 


QC 

used in Carter coefficient calculation 


QN 

slots per pole per phase 

+ 

QPERV 

voltage points on no-load saturation curve, percent 


QQ 

number of stator slots 

+ 

QVLL 

line-to-line voltage on no-load saturation curve, V rms 

+ 

QVLN 

line-to-neutral voltage on no-load saturation curve, V rms 


RATING 

NAMELIST name 


RB 

resistivity of stator conductor at temperature TTA, /i ohm -in, 


RC 

field coil weight, lb 

* 

RD 

field conductor diameter or width, in. 


RE1 

Reynolds number for main air gap 


RE2 

Reynolds number in auxiliary air gap 

* 

RF 

type of stator winding (random or form wound) 


RG1 

stator winding resistance at temperature TST, ohms 

*+ 

RMAT 

rotor material name 


ROTOR 

NAMELIST name 

* 

RPM 

rotor speed, rpm 

* 

RR 

field coil resistivity at 20° C, ^ohm-in. 

+ 

RRA 

armature resistance, ohm 

+ 

RRB 

field coil resistance of NFC coils in series, ohms 
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* 

RS 

stator conductor resistivity at 20° C, /^ohm-in. 

* 

RT 

field conductor thickness, in. 


RY 

RRA(J)/RB, M in. _1 


S 

2 

armature current density at rated load, A/in. 

* 

SC 

number of conductors per stator slot 


SCR 

short circuit ratio 

* 

SF 

distance between centerlines of stator strands (depth) (see fig. 24), in. 

* 

SF 

stator stacking factor 


SH 

uninsulated strand height, in. 

* 

SHL1 

shaft length (see fig. 23), in. 

* 

SHL2 

shaft length (see fig. 23), in. 


SI 

field self- inductance, H 


SIGMA 

2 

specific tangential force on rotor, lb/in. 

* 

SK 

stator slot skew at stator inside diameter, in. 


SLOTS 

NAMELIST name 


SM 

tooth width at 1/3 distance from narrowest section, in. 

*+ 

SMAT 

stator material name 

* 

SN 

strands per stator conductor in depth 

* 

SN1 

total strands per stator conductor 

+ 

SP 

total losses, W 


SPS 

special permeance (see ref. 2), lines/A-turn 


ss 

solid stator stack length, in. 

+ 

ST 

stator tooth loss, W 

+ 

STAR 

used in subroutine OUTPUT to draw diagram 


STATET 

number of effective aramature winding turns 


STATOR 

NAMELIST name 


SUBLUN 

subroutine name 

* 

T1 

armature winding temperature at rated load, °C 

* 

Til 

armature winding temperature at no-load, °C 

* 

T2 

field winding temperature at rated load, °C 
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* 

T22 


TC 

* 

TF 


TG 


TP 


TS 

* 

TST 


TT 

+ 

TTA 


TTB 

* 

TY 

* 

TY2 

* 

TY3 

* 

TYE 

* 

TYL 

* 

TYPY 

* 

TYR 


U 


UA 

* 

VA 

* 

VMIN 


VR 


vv 

+ 

WA 


wc 

* 

WF 


WFl 


WF2 


WI 


field winding temperature at no-load, °C 

open-circuit time constant (field only), sec 

field coil temperature at which FK1 and TC are calculated, °C 

total air-gap flux at no-load rated voltage, kilolines 

pole pitch, in. 

stator slot pitch at stator inside diameter, in. 

armature winding temperature at which RG1 is calculated, °C 

stator slot pitch at 1/3 distance from narrow section, in. 

armature winding temperature, °C 

field winding temperature, °C 

yoke dimension (see fig. 22), in. 

yoke dimension (see fig. 22), in. 

yoke dimension (see fig. 22), in. 

yoke dimension (see fig. 22), in. 

yoke dimension (see fig. 22), in. 

field coil location code (see fig. 8) 

yoke dimension (see fig. 22), in. 

used to calculate end extension length 

=G(M), per unit 

rating of alternator, kVA 

minimum voltage at which no-load saturation curve is calculated, 
per unit 

rotor peripheral velocity at main air gap, ft/min 

rating of alternator, W 

alternator output power, kW 

armature winding weight (copper only), lb 

total windage loss, W 

windage loss in main air gap, W 

windage loss in auxiliary air gap, W 

stator iron weight, lb 
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WINDGE 

WHSTDNG 

* WL 

+ WMIS 

WN 
WQ 

+ WQL 

* WROTOR 
WT 

WTOTAL 

* WYOKE 

X 

XA 

XB 

XD 

XF 

XI 

XL 

XQ 

XR 

XU 

XX 

Y 

*+ YMAT 

YOKE 

* YY 
Z 


subroutine name 
NAMELIST name 

core loss at flux density BK, W/lb 
miscellaneous load losses, W 
no-load pole-face loss, W 
no-load rated-voltage core loss, W 
core loss, W 
rotor weight, lb 

no-load rated-voltage tooth loss, W 
total electromagnetic weight, lb 
yoke weight, lb 

flux density at which AT is found by interpolation (see fig. 27), 
2 

kilolines/in. 

synchronous reactance (direct), percent 
synchronous reactance (quadrature), percent 
armature reaction reactance (direct), percent 
field winding leakage reactance, percent 

angle between line current phasor and D-axis (internal power 
factor angle) (see fig. 7), rad 

stator winding leakage reactance, percent 

armature reaction reactance (quadrature), percent 

reactance factor 

transient reactance (direct), percent 

slope of magnetization curve at flux density X (see fig. 27) 
used in interpolation procedure of material curves 
yoke material name 
NAMELIST name 

slots spannd per coil (=No. of slots between coil sides + 1) 
dummy variable used in slot leakage permeance calculation 
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ZA 

dummy variable used in slot leakage permeance calculation 


ZB 

dummy variable used in slot leakage permeance calculation 


ZC 

dummy variable used in slot leakage permeance calculation 


ZD 

dummy variable used in slot leakage permeance calculation 


ZE 

dummy variable used in slot leakage permeance calculation 


ZF 

dummy variable used in slot leakage permeance calculation 


ZG 

(length of field winding conductor)/ (field conductor area), in. 


ZY 

0.7* HS, in. 


ZY 

used as an index in armature and field conductor area calculations 

* 

ZZ 

type of slot (see fig. 21) 



FORTRAN VARIABLES USED IN WINDGE 


CDN 

skin friction coefficient (new value in iteration) 


CDD 

skin friction coefficient (old value in iteration) 


L 

air-gap length (axial direction), in. 

* 

M 

molecular weight of fluid in air gap 


NU 

o 

kinematic viscosity, in. /sec 


P 

windage power loss, W 

* 

PRESS 

2 

pressure of gas in air gap, lb/in. 


R 

radius of rotor at air gap, in. 


RE 

Reynolds number 

* 

RO 

density of gas in air gap at temperature TEMP and pressure 
PRESS, lb/ft 3 


RPM 

rotor speed, rpm 


T 

air-gap length (radial direction), in. 

* 

TEMP 

temperature of gas in air gap, °C 
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VIS fluid viscosity, lb/ft-sec 

WIND NAMELIST name 

WINDGE subroutine name 
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APPENDIX E 


ROTOR LEAKAGE PERMEANCE MEASUREMENT 

The permeance & of a section of flux path is defined as the ratio 

&=1L 

Ni 

where <p is the flux through the section of the flux path under consideration and Ni is 
the ampere-turn drop across the same flux path section. 

In a Lundell alternator, a permeance that is of particular interest is the rotor leak- 
age permeance (^ L ). This is the permeance of the path followed by the flux that goes 
from one rotor half to the other without entering the stator. This flux is termed leakage 
flux in that it bypasses the armature and therefore does not generate a voltage in 
the armature winding. The magnitude of the rotor leakage permeance has a large 

effect on the calculation of field excitation requirements. In the evaluation of the com- 


r Winding (190 
\ turns, 13AWG) 



C -69-2113 


Figure 28. - Rotor leakage permeance measuring fixture. 
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puter program accuracy, it was of some interest, therefore, that an accurate value be 
known. 

To obtain an accurate value for a fixture was built which would allow measur- 

ing the leakage permeance of the rotor. This fixture, with the rotor inserted, is shown 
in figure 28. Essentially, the fixture establishes a well defined, low reluctance flux 
path from one end of the rotor to the other. The distributed winding around the yoke of 
the fixture provides a means of exciting the magnetic circuit, which consists of the yoke, 
an air gap, and the rotor in series. The purpose of the air gap is to allow a gaussmeter 
probe to be inserted into the magnetic circuit to measure flux densities within the yoke. 

To obtain using this fixture, the ampere-turn drop across and the total flux 

through the nonmagnetic section of the rotor are measured. These values are then sub- 
stituted into the basic definition of permeance. 

The ampere-turn drop across the nonmagnetic rotor section is measured as follows. 
A dummy rotor made of solid SAE 1010 steel in the shape of a cylinder is inserted in the 
fixture in place of the actual rotor. The winding is then excited and a curve of excitation 
current against air gap flux density is obtained. This is the lower curve shown in fig- 
ure 29. Essentially, this curve gives, at any given value of air gap flux density, the 



Figure 29. - Experimental data used to obtain ampere- 
turn drop across nonmagnetic rotor section. Wind- 
ing used to obtain these curves had 190 turns. 
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ampere-turn drop across the yoke, the air gap, and the magnetic rotor sections. A simi- 
lar curve is then obtained with the actual alternator rotor inserted in the fixture. This 
is the upper curve shown in figure 29. At any value of air-gap flux density, the differ- 
ence between the two curves gives the ampere-turn drop across the nonmagnetic section 
of the alternator rotor. For example, at the air-gap flux densities of 1. 5 and 2. 5 kilo- 
gauss, the ampere-turn drop across the rotor leakage permeance is 351 and 627, re- 
spectively. 

The flux through the nonmagnetic rotor section is measured at the air-gap in the 

yoke. It consists of two components: the flux through the air-gap proper and the flux 

due to fringing around the air-gap. The flux in the air-gap proper is merely the product 

2 

of the air-gap area (3. 31 in. ) and the flux density, which is constant throughout the air 
gap. The flux in the fringe pattern was obtained by measuring the flux density at many, 
regularly spaced grid points in an approximately 14-square-inch area surrounding the 
yoke at the air gap. The results of these measurements are as follows: at an air-gap 
flux density of 1. 5 kilogauss, the total flux through the yoke is 38. 9 kilolines; at 
2. 5 kilogauss air-gap flux density, the total yoke flux is 66. 5 kilolines. 

When these values, and the values for ampere-turn drop are substituted into the 
basic definition for permeance, the result is, for an air-gap flux density of 1. 5 kilogauss, 

& - 3^- 9 X 10 _ in lines/ampere-turn 

RL 351 

or, for an air-gap flux density of 2. 5 kilogauss, 

g, _ 66. 5x10 _ jQg ij nes / am p ere _i urn 
RL, g27 

The difference between the two values is within the experimental error and will be 

taken to be the average, that is, = lines per ampere-turn. 

An assumption made in these calculations is that all the flux in the air gap of the 
fixture, including that due to fringing, passes through the nonmagnetic rotor section. 

This is, of course, not completely correct. Some of the flux leaks across the fixture 
without ever entering the rotor. Some additional flux leaks around the nonmagnetic 
material and into the space normally occupied by the stator. This flux, as regards the 
alternator, is not leakage flux. The import of this is that the measured value, obtained 
as described in this appendix, is higher than the actual rotor leakage permeance. 
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